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0\ ' Abstract 
(N 

\ The present experimental results on neutrino flavour oscillations provide evidence 

^> ' for non-zero neutrino masses, but give no hint on their absolute mass scale, which 

. is the target of beta decay and neutrinoless double-beta decay experiments. Crucial 

\ complementary information on neutrino masses can be obtained from the analysis of 

' data on cosmological observables, such as the anisotropics of the cosmic microwave 

^ ■ background or the distribution of large-scale structure. In this review we describe 

^) \ in detail how free-streaming massive neutrinos affect the evolution of cosmological 



perturbations. We summarize the current bounds on the sum of neutrino masses 



Q^' that can be derived from various combinations of cosmological data, including the 



most recent analysis by the WMAP team. We also discuss how future cosmological 



|H \ experiments are expected to be sensitive to neutrino masses well into the sub-eV 

^ ■ range. 



Key words: Neutrino masses, Cosmology, Dark matter 
PACS: 14.60.Pq, 95.35.-hd, 98.80.Es 



Contents 

Introduction! 3 



1 



2 Neutrino oscillations and absolute neutrino mass searches! 5 



3 The cosmic neutrino background! 9 



3.1 Basics on relic neutrinos, including neutrino decoupling! 9 



3.2 Extra radiation and the effective number of neutrinos! 13 



Email addresses: lesgourg@lapp.in2p3.fr (Julien Lesgourgues), 
pastor@ific.uv.es (Sergio Pastor). 



Preprint submitted to Elsevier Science 



3 February 2008 



3.3 Massive neutrinos as dark mattei 



4 Massive neutrinos and cosmological perturbations! 



4.1 Qbservables targets: definition of the power spectral 



4.2 Baclcground evolutioni 



4.3 Gauge transformations and Einstein eauationsi 



4.4 Linear perturbation theorv in a neutrinoless Universe (pure ACDmH 



4.5 Linear perturbation theory in presence of neutrinos fAMDM) 



4.6 Summary of the neutrino mass efFectsI 



5 Current observations and boundsl 



5.1 Statistical methods for parameter inference from cosmological datal 



5.2 CMB anisotropiesi 



5.3 Galaxy redshift survevd 



5.4 Lvman-g forest 



5.5 Summary of current bound'3 



5.6 Extra parameters! 



5.7 Non-standard relic neutrino^ 



6 Future sensitivities and new experimental techniauesi 



6.1 Fisher matrix! 



6.2 Future CMB experiments! 



6.3 Future galaxy redshift survevi 



6.4 CMB weak lensing| 



6.5 Galaxy weak lensing (cosmic shear! 



6.6 Galaxy cluster survey^ 



7 Conclusion^ 



Referenced 



1 Introduction 



Neutrino cosmology is a fascinating example of the fecund interaction between 
particle physics and astrophysics. At the present time, the researchers working 
on neutrino physics know that the advances in this field need a combined effort 
in the two areas. 

Prom the point of view of cosmologists, the idea that massive neutrinos could 
play a significant role in the history of the Universe and in the formation of 
structures has been discussed for more than thirty years, first as a pure spec- 
ulation. However, nowadays we know from experimental results on flavour 
neutrino oscillations that neutrinos are massive. At least two neutrino states 
have a large enough mass for being non-relativistic today, thus making up a 
small fraction of the dark matter of the Universe. At a stage in which cos- 
mology reaches high precision thanks to the large amount of observational 
data, it is unavoidable to take into account the presence of massive neutrinos. 
In particular, the observable matter density power spectrum is damped on 
small scales by massive neutrinos. This effect can range from a few per cent 
for 0.05 — 0.1 eV masses, the minimal values of the total neutrino mass com- 
patible with oscillation data, up to 10 — 20% in the limit of three degenerate 
masses. 

From the point of view of particle physicists, fixing the absolute neutrino mass 
scale (or, equivalently, the lightest neutrino mass once the data on flavour os- 
cillations are taken into account) is the target of terrestrial experiments such 
as the searches for neutrinoless double beta decay or tritium beta decay ex- 
periments, a difficult task in spite of huge efforts and very promising sched- 
uled experiments. At the moment, the best bounds come from the analysis 
of cosmological data, from the requirement that neutrinos did not wash out 
too much of the small-scale cosmological structures. Recently the cosmolog- 
ical limits on neutrino masses progressed in a spectacular way, in particular 
thanks to the precise observation of the Cosmic Microwave Background (CMB) 
anisotropics by the Wilkinson Microwave Anisotropy Probe (WMAP) satel- 
lite, and to the results of a new generation of very deep galaxy redshift surveys. 
Neutrino physicists look with anxiety at any new development in this field, 
since in the next years many cosmological observations will be available with 
unprecedented precision. 

At this very exciting moment, which might be preceding a breaking discovery 
within a few years, our aim is to present here a clear and comprehensive re- 
view of the role of massive neutrinos in cosmology, including a description of 
the underlying theory of cosmological perturbations, a summary of the cur- 
rent bounds and a review of the sensitivities expected for future cosmological 
observations. Our main goal is to address this manuscript simultaneously to 
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particle physicists and cosmologists: for each aspect, we try to avoid useless 
jargon and to present a self-contained summary. Compared to recent short 
reviews on the subject, such as Refs. [1,2,3,4], we tried to present a more de- 
tailed discussion. At the same time, we will focus on the case of three flavour 
neutrinos with masses in accordance with the current non-cosmological data. 
A recent review on primordial neutrinos can be found in [5], while for many 
other aspects of neutrino cosmology we refer the reader to [6]. Finally, a more 
general review on the connection between particle physics and cosmology can 
be found in [7]. 

We begin in Sec. 2 with an introduction to flavour neutrino oscillations, and 
their implications for neutrino masses in the standard three-neutrino scenario. 
We also briefly review the limits on neutrino masses from laboratory exper- 
iments. Then, in Sec. 3, we describe some basic properties of the Cosmic 
Neutrino Background. According to the Big Bang cosmological model, each of 
the three flavour neutrinos were in thermal equilibrium in the Early universe, 
and then decoupled while still relativistic. We review the consequences of this 
simple assumption for the phase-space distribution and number density of the 
relic neutrinos, and we summarize how this standard picture is confirmed by 
observations from Big Bang Nucleosynthesis (BBN) and cosmological pertur- 
bations. Finally we explain why neutrinos are expected to constitute a fraction 
of the dark matter today. 

In Sec. 4, we describe the impact of massive neutrinos on cosmological pertur- 
bations. This section is crucial for understanding the rest, since current and 
future bounds rely precisely on the observation of cosmological perturbations. 
Although the theory of cosmological perturbations is a very vast and technical 
topic, we try here to be accessible both to particle physicists willing to learn 
the field, and to cosmologists willing to understand better the specific role of 
neutrinos. We introduce the minimal number of concepts, technicalities and 
equations for understanding the main results related to neutrinos. Still this 
section is quite long because we wanted to make it self-contained, and to ex- 
plain many intermediate steps that are usually hidden in most works on the 
subject. 

In Sec. 4.1, we define the quantities that can be observed, i.e. those that we 
want to compute. In Sec. 4.2, we describe the evolution of homogeneous quan- 
tities and in Sec. 4.3 we define the general setup for the theory of cosmological 
perturbations. Then in Sec. 4.4 we described in a rather simplified way what 
the evolution of perturbations would look like in absence of neutrinos. For a 
cosmologist, these four sections are part of common knowledge and can be 
skipped. In Sec. 4.5, we present a detailed description of the impact of neu- 
trinos on cosmological perturbations. This long section is the most technical 
part of this work, so the reader who wants to avoid technicalities and to know 
the results can go directly to Sec. 4.6, for a comprehensive summary of the 
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effects of neutrino masses on cosmological observables. 

Then, in Sec. 5 we review the current existing bounds on neutrino masses, 
starting from those involving only CMB observations, and then adding differ- 
ent sets of data on the distribution of the Large Scale Structure (LSS) of the 
Universe, which include those coming from galaxy redshift surveys and the 
Lyman-a forest. We explain why a unique cosmological bound on neutrino 
masses docs not exist, and why there are significant variations from paper to 
paper, depending on the data included and the assumed cosmological model. 

Finally, in Sec. 6 we describe the prospects for the next ten or fifteen years. 
We summarize the forecasts which can be found in the literature concerning 
the sensitivity to neutrino masses of future CMB experiments, galaxy redshift 
surveys, weak lensing observations (either from the CMB or from galaxy el- 
lipticity) and galaxy cluster surveys. We conclude with some general remarks 
in Sec. 7. 



2 Neutrino oscillations and absolute neutrino meiss seeirches 

Neutrinos have played a fundamental role in the understanding of weak inter- 
actions since they were postulated by W. Pauli in 1930 to safeguard energy 
conservation in beta decay processes. These chargeless leptons are massless 
in the framework of the succesful Standard Model (SM) of particle physics. 
However, this is an accidental prediction of the SM, and there are many well- 
motivated extended models where neutrinos acquire mass and other non-trivial 
properties (see e.g. [8,9,10,11]). Thus the measurement of neutrino masses 
could give us some hints on the new fundamental theory, of which the SM is 
just the low-energy limit, and in particular on new energy scales. 

It was realized by B. Pontecorvo in 1957 that if neutrinos were massive there 
could exist processes where the neutrino fiavour is not conserved, that we call 
neutrino oscillations. For small neutrino masses, these oscillations actually 
take place on macroscopic distances and can be measured if we are able to 
detect neutrinos from distant sources, identify their flavour and compare the 
results with the theoretical predictions for the initial neutrino fluxes. Actually, 
after decades of experimental efforts in underground facilities, neutrinos have 
been detected from various natural or artificial sources. The former include the 
neutrinos produced in the nuclear reactions in the Sun, as secondary particles 
from the interaction of cosmic rays in the atmosphere of the Earth and even 
a few neutrinos were detected from a supernova explosion (SN1987A). We 
have also measured the neutrino fiuxes originated in artificial sources such as 
nuclear reactors or accelerators. 



5 



atm 



m 


m 

3 


2 




Am2 < 




sun 

2 


1 








1 


3 



Am" 



NORMAL 



INVERTED 



Fig. 1. The two neutrino schemes allowed if Aml^^ S> Arn^^^: normal hierarchy 
(NH) and inverted hierarchy (IH). 

Nowadays there exist compelling evidences for flavour neutrino oscillations 
from a variety of experimental data on solar, atmospheric, reactor and accel- 
erator neutrinos. These are very important results, because the existence of 
flavour change implies that neutrinos mix and have non-zero masses, which 
in turn requires particle physics beyond the SM. There are many excellent 
reviews on neutrino oscillations and their implications, to which we refer the 
reader for more details (see e.g. the recent ones [13,14]). 

We know that the number of light neutrinos sensitive to weak interactions 
(flavour or active neutrinos) equals three from the analysis of the invisible Z- 
boson width at LEP, N^, = 2.994 ± 0.012 [12], and the three flavour neutrinos 
(i/g, I'n, Vr) s-re linear combinations of states with deflnite mass i/j, where i is 
the number of massive neutrinos. 

In a three-neutrino scenario flavour and mass eigenstates are related by the 
mixing matrix [/, parametrized as [12,14] 



C12C13 



~Si2C23 — Ci2S23Sl3e 



id 

C12C23 ~ 5125235136 523Cl3 



512523 ~ Ci2C23>Si3e*^ " Ci2523 — 5i2C23<Si3e'^ C23C13 



(1) 



Here Cij = cos 6ij and Sij = sin 6ij for ij = 12, 23 or 13, and 5 is a CP-violating 
phase. Together with the three masses, in total there are seven flavour param- 
eters in the neutrino sector. This would be all if neutrinos were Dirac particles, 
like the charged leptons, and the total lepton number is conserved. Instead, if 
neutrinos are Majorana particles (i.e. a neutrino is its own antiparticle) , the 
matrix U is multiplied by a diagonal matrix of phases that can be taken as 
diag(l, e**^^/^, 6**^'^^+^'^)/^). These two phases do not show up in neutrino oscil- 
lations, but appear in lepton-number-violating processes such as neutrinoless 
double beta decay, as discussed later. 
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Fig. 2. Expected values of neutrino masses according to the values in Eq. (2). Left: 
individual neutrino masses as a function of the total mass for the best-fit values of 
the Am?. Right: ranges of total neutrino mass as a function of the lightest state 
within the 3a regions (thick lines) and for a future determination at the 5% level 
(thin lines). 

Oscillation experiments can measure the differences of squared neutrino masses 
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the relevant ones for solar and atmo- 



spheric neutrinos, respectively. As a reference, we take the following 3a ranges 
of mixing parameters from an update of ref. [13], 



Aml^ = (7.9+^:°) X 10-^ eV^ \Aml^\ = (2.2+^]) x 10"=^ eV^ 

4 = 0.30^°:^° 4 = o.5ol°:i^ 4 < 0.043 (2) 

Unfortunately oscillation experiments are insensitive to the absolute scale of 
neutrino masses, since the knowledge of Am^i > and |Am|]^| leads to the two 
possible schemes shown in Fig. 1, but leaves one neutrino mass unconstrained 
(see e.g. the discussion in the reviews [14,15,16,17,18]). These two schemes 
are known as normal (NH) and inverted (IH) hierarchies, characterized by 
the sign of Arn^^, positive and negative, respectively. For small values of the 
lightest neutrino mass mo, i.e. mi (ms) for NH (IH), the mass states follow 
a hierarchical scenario, while for masses much larger than the differences all 
neutrinos share in practice the same mass and then we say that they are 
degenerate. In general, the relation between the individual masses and the 
total neutrino mass can be found numerically, as shown in Fig. 2. 

It is also possible that the number of massive neutrino states is larger than 
the number of flavor neutrinos. In such a case, in order to not violate the 
LEP results the extra neutrino states must be sterile, i.e. singlets of the SM 
gauge group and thus insensitive to weak interactions. At present, the results 
of the Liquid Scintillator Neutrino Detector (LSND) [19], an experiment that 
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has measured the appearance of electron antineutrinos in a muon antineutrino 
beam, constitute an independent evidence of neutrino conv6rsions dbt db l&,rff6r 
mass difference than those in Eq. 2. In such a case, a fourth sterile neutrino 
is required with mass of C(eV) [13]. The LSND results will be checked by the 
ongoing MiniBoone experiment [20], whose first data are expected for 2006. In 
this review, we will mainly consider the three-neutrino scenario, but we briefly 
comment on the cosmological bounds on the four-neutrino mass schemes that 
include the LSND results in Sec. 5.6. 

As we discuss in the next sections, cosmology is at first order sensitive to 
the total neutrino mass if all states have the same number density, provid- 
ing information on rriQ but blind to neutrino mixing angles or possible CP 
violating phases. Thus cosmological results are complementary to terrestrial 
experiments such as beta decay and neutrinoless double beta decay, which are 
respectively sensitive to the effective masses 



/ N 1/2 

I |7-r|2 21 /22 2,22 2,2 2\l/2 

^f^^yl^l Uei \ rriij = (C12C13 mi + s^^^is ^2 + ^is ^3) ' > 



A A \ A A 109 I A 



(3) 



The search for a signal of non-zero neutrino masses in a beta decay exper- 
iment is, in principle, the best strategy for measuring directly the neutrino 
mass, since it involves only the kinematics of electrons (see [21] for a recent 
review) . The current limits from tritium beta decay apply only to the range of 
degenerate neutrino masses, so that ~ mo (if unitarity of U is assumed). 
The bound at 95% CL is mo < 2.05 - 2.3 eV from the Troitsk [22] and Mainz 
[23] experiments, respectively. This value is expected to be improved by the 
KATRIN project [24] to reach a discovery potential for 0.3 — 0.35 eV masses 
(or a sensitivity of 0.2 eV at 90% CL). 

The neutrinoless double beta decay {Z, A) {Z + 2, A) + 2e~ (in short Oi/2/3) 
is a rare nuclear processes where lepton number is violated and whose obser- 
vation would mean that neutrinos are Majorana particles. If the Oz/2/3 process 
is mediated by a light neutrino, the results from neutrinoless double beta de- 
cay experiments are converted into an upper bound or a measurement of the 
effective mass mpp in Eq. (3). A recent work [25] gives the following upper 
bounds at 99% CL 

\mpp\ < (0.44 - 0.62) Hn eV (4) 



where /i at is a parameter that characterizes the uncertainties on the corre- 
sponding nuclear matrix elements (see [26] for a recent discussion and [15] 
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for a review). The above range corresponds to the resuhs of recent double 
beta decay experiments based on ^^Ge, such as Heidelberg-Moscow [27] and 
IGEX [28], or ^-^"Te such as Cuoricino [29] (see e.g. [30] for results from other 
isotopes). In addition, there exists a claim of a positive 0i/2/? signal which 
would correspond to the approximate range 0.1 < |m/3/3|/eV < 0.9 [31]. Fu- 
ture Oz/2/9 projects will improve the current sensitivities down to values of 
the order jm^g^l ~ 0.01 — 0.05 eV [15]. An experimental detection of Oz/2/9, 
in combination with results from beta decay or cosmology, will also help us 
to discriminate between the two neutrino mass spectra and to pin down the 
values of the neutrino flavour parameters, as discussed for instance in [32,33]. 

There are other ways to obtain information on the absolute scale of neutrino 
masses, such as the measurement of the time-of-flight dispersion of a supernova 
neutrino signal (at most sensitive to masses of order some eV, see [34] and 
references therein), but to conclude this section let us note that the sum of 
neutrino masses is restricted to the approximate range 



where the upper limit comes exclusively from tritium beta decay results and 
the lower limit reflects the minimum values of the total neutrino mass in the 
normal (inverted) hierarchy. 



3 The cosmic neutrino background 

3.1 Basics on relic neutrinos, including neutrino decoupling 

The existence of a relic sea of neutrinos is a generic prediction of the standard 
hot big bang model, in number only slightly below that of relic photons that 
constitute the CMB. Produced at large temperatures by frequent weak interac- 
tions, cosmic neutrinos were kept in equilibrium until these processes became 
ineffective in the course of the expansion. While coupled to the rest of the 
primeval plasma, neutrinos had a momentum spectrum with an equilibrium 
Fermi-Dirac form with temperature T, 



Here we have included a neutrino chemical potential /i^ that would exist in the 
presence of a neutrino-antineutrino asymmetry, but it was shown in [35,36,37] 
that the stringent BBN bounds on ii^,^ apply to all flavours, since neutrino 



0.056 (0.095) eV < ^ 



"^i ^ 6 eV 



(5) 




(6) 
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oscillations lead to approximate flavour equilibrium before BBN. The current 
bounds on the common value of the neutrino degeneracy parameter = fJ'u/T 
are —0.05 < < 0.07 at 2a [38] (relaxed to —0.13 < < 0.3 in presence of 
extra relativistic degrees of freedom [39,40]). Thus the contribution of a rehc 
neutrino asymmetry can be safely ignored. 

As the universe cools, the weak interaction rate falls below the expansion 
rate given by the Hubble parameter H and one says that neutrinos decouple 
from the rest of the plasma. An estimate of the decoupling temperature can 
be found by equating the thermally averaged value of the weak interaction 
rate 



where oc G'j^ is the cross section of the electron-neutrino processes with Gp 
the Fermi constant and Ui, is the neutrino number density, with the expansion 
rate 



where p is the total energy density and Mp is the Planck mass. If we ap- 
proximate the numerical factors to unity, with F,^ ^ G'j^T^ and H ^ T'^/Mp, 
we obtain the rough estimate Tdec ~ 1 MeV. More accurate calculations give 
slightly higher values of T^ec which are flavour dependent since electron neu- 
trinos and antineutrinos are in closer contact with e^, as shown e.g. in [6]. 

Although neutrino decoupling is not described by a unique Tdec, it can be ap- 
proximated as an instantaneous process. The standard picture of instantaneous 
neutrino decoupling is very simple (see e.g. [41,42]) and reasonably accurate. 
In this approximation, the spectrum in Eq. (6) is preserved after decoupling, 
since both neutrino momenta and temperature redshift identically with the 
universe expansion. In other words, the number density of non-interacting 
neutrinos remains constant in a comoving volume since the decoupling epoch. 

We have seen in Sec. 2 that active neutrinos cannot possess masses much 
larger than 1 eV, so they were ultra-relativistic at decoupling. This is the 
reason why the momentum distribution in Eq. (6) does not depend on the 
neutrino masses, even after decoupling, i.e. there is no neutrino energy in the 
exponential of /eq(p). 

When calculating quantities related to relic neutrinos, one must consider the 
various possible degrees of freedom per flavour. If neutrinos are massless or 
Major ana particles, there are two degrees of freedom for each flavor, one for 




(7) 




(8) 
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neutrinos (one negative helicity state) and one for antineutrinos (one posi- 
tive helicity state). Instead, for Dirac neutrinos there are in principle twice 
more degrees of freedom, corresponding to the two hehcity states. However, 
the extra degrees of freedom should be included in the computation only if 
they are populated and brought into equilibrium before the time of neutrino 
decoupling. In practice, the Dirac neutrinos with the "wrong-helicity" states 
do not interact with the plasma at temperatures of the MeV order and have 
a vanishingly small density with respect to the usual left-handed neutrinos 
(unless neutrinos have masses close to the keV range, as explained in Sec. 6.4 
of [6], but such a large mass is excluded for active neutrinos). Thus the relic 
density of active neutrinos does not depend on their nature, either Dirac or 
Majorana particles. 

Shortly after neutrino decoupling the photon temperature drops below the 
electron mass, favouring annihilations that heat the photons. If one as- 
sumes that this entropy transfer did not affect the neutrinos because they 
were already completely decoupled, it is easy to calculate the ratio between 
the temperatures of relic photons and neutrinos T^/Ty = (11/4)^/^ ~ 1.40102. 
Any quantity related to relic neutrinos can be calculated at later times with 
the spectrum in Eq. (6) and T,^. For instance, the number density per flavour 
is fixed by the temperature, 

3 6C(3) „ 



which leads to a present value of 113 neutrinos and antineutrinos of each 
flavour per cm'^. Instead, the energy density for massive neutrinos should in 
principle be calculated numerically, with two well-defined analytical limits. 



Pv{my > Ty)^myny . (10) 

Thus we see that the contribution of massive neutrinos to the energy density 
in the non-relativistic hmit is a function of the mass (or the sum of masses if 
all neutrino states have rrii » T^). 

In a more accurate analysis of neutrino decoupling, the standard picture de- 
scribed above is modified: the processes of neutrino decoupling and an- 
nihilations are sufficiently close in time so that some relic interactions be- 
tween and neutrinos exist. These relic processes are more efficient for 
larger neutrino energies, leading to non-thermal distortions in the neutrino 
spectra and a slightly smaller increase of the comoving photon temperature, 
as noted in a series of works (see the full list given in the review [6]). A 
proper calculation of the process of non-instantaneous neutrino decoupling 
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Fig. 3. The upper (lower) lines are the frozen distortions of the electron (muon or 
tau) neutrino spectra as a function of the comoving momentum y, calculated with 
(solid) and without (dotted) the effect of flavour oscillations. Left: real neutrino 
distribution functions normalized to the equilibrium one. Right: contribution of the 
distortions to the comoving number density. Here the scale factor was normalized 
so that a(t) l/T^ at large temperatures. 

demands solving the Boltzmann equations for the neutrino spectra, a set of 
integro-differential kinetic equations that are difficult to solve numerically. The 
momentum-dependent calculations were carried out in refs. [43,44,45], while 
the inclusion of finite temperature QED corrections to the electromagnetic 
plasma was done in [46]. 



A recent work [47] has considered the effect of fiavour neutrino oscillations 
on the neutrino decoupling process (see also [48]). The frozen values of the 
neutrino distributions are shown in Fig. 3, for the cases calculated with and 
without the effect of fiavour oscillations. One can see that the distortions grow 
with the neutrino momentum, and their contribution to the number density 
of neutrinos is maximal around ?/ = 4. For the best-fit values of the mixing 
parameters in Eq. (2), Ref. [47] found an increase in the neutrino energy den- 
sities of 0.73% and 0.52% for i/g's and t'^,r's, respectively. At the same time, 
after annihilations the comoving photon temperature is a factor 1.3978 
larger, instead of 1.40102 in the approximation of instantaneous decoupling. 
Note, however, that for any cosmological epoch when neutrino masses can be 
relevant one should consider the spectra of the neutrino mass eigenstates 2^1,2,3 
instead of fiavour eigenstates t'e./x.T- In principle, for numerical calculations of 
the cosmological perturbation evolution such as those done by the Boltzmann 
codes CMBFAST [49] or camb [50], one should include the full distribution 
function of each mass eigenstate, including the distortions shown in Fig. 3. 
It is only when the neutrinos are relativistic or when one is interested in an 
observable that does not depend on neutrino masses (such as the number den- 
sity) that the effect of the distortion can be simply integrated over momentum. 
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For instance, the contribution of relativistic relic neutrinos to the total energy 
density is taken into account just by using A^^cfr = 3.046 as defined later in 
Eq. (11). In practice, the distortions calculated in Ref. [47] only have small 
consequences on the evolution of cosmological perturbations, and for many 
purposes they can be safely neglected. 

3.2 Extra radiation and the effective number of neutrinos 

Neutrinos fix the expansion rate during the cosmological era when the Universe 
is dominated by radiation. Their contribution to the total radiation content 
can be parametrized in terms of the effective number of neutrinos A'eff [51,52], 
through the relation 



where is the energy density of photons, whose value today is known from the 
measurement of the CMB temperature. This equation is valid when neutrino 
decoupling is complete and holds as long as all neutrinos are relativistic. 

In Sec. 2 we saw that from accelerator data the number of active neutrinos 
is three, while in the previous subsection we learned from the analysis of 
neutrino decoupling these three active neutrinos contribute as N^s = 3.046. 
Any departure of A'efr from this last value would be due to non-standard 
neutrino features or to the contribution of other relativistic relics. 

A detailed discussion of cosmological scenarios where A^^cfr is not fixed to three 
can be found in the reviews [6,53], while the particular case of active-sterile 
neutrino mixing was recently analyzed in [54]. Since in the present work we 
focus on the standard case of three active neutrinos, here we only give a brief 
review of the most recent bounds on A^eflf from cosmological data. 

The value of N^^ is constrained at the BBN epoch from the comparison of 
theoretical predictions and experimental data on the primordial abundances 
of light elements, which also depend on the baryon-to-photon ratio rjb — n\)/nj 
(or baryon density). The main effect of N^s is to fix the Hubble expansion rate 
through its contribution to the total energy density. This in turn changes the 
freezing temperature of the neutron-to- proton ratio, therefore producing a 
different abundance of ^He. 

The BBN bounds on Neg have been recently reanalyzed taking in input the 
value of the baryon density derived from the WMAP first year data [55] 
VCMB — 6.14 ± 0.25. In Fig. 4 we show the results from [40], where the al- 
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Fig. 4. Bounds on A'efr from BBN including D and ^He (contours at 68% and 
95% CL) and from a combined analysis of BBN (D only) and CMB data. Here 
= ^hh"^ - 10i°(77b/274). This Figure is taken from Ref. [40]. 



lowed raiiGre N,s = 2.5lJi (95% CL) was inferred from data on light element 
abundances (see also [56,57,58] and [59] for a recent review on BBN). This 
range is perfectly compatible with the standard prediction of 3.046. However, 
the reader should be cautious in the interpretation of the BBN allowed range 
for such as that in Fig. 4. It is well-known that the main problem when 
deriving the primordial abundances from observations in astrophysical sources 
is the existence of systematics not accounted for, in particular for "^He (a dis- 
cussion of the observational situation can be found in Ref. [60]). For instance, 
the detailed BBN analysis in [61] gives the allowed range 1.90 < Ncs < 3.77 
for the *^He "box" range in [60]. 

Independent bounds on the radiation content of the universe at a later epoch 
can be extracted from the analysis of the power spectrum of CMB anisotropies. 
We will describe later the neutrino effects on the CMB spectrum. Assum- 
ing a minimal set of cosmological parameters and a flat universe, Ref. [62] 
found the value N^s = 3.512^ (95% CL) (see also [56,57,63]) from the com- 
bination of WMAP first year data with other CMB and LSS data. More re- 
cently, it has been shown that the addition of Supernovae and new CMB data 
leads to a reduction of the allowed range : iVeff = 4.2^1:7 [64] (with the latest 
BOOMERANG data) and N^s = 3.3l^i [65] (with WMAP three year data). 
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3.3 Massive neutrinos as dark matter 

A priori, massive neutrinos are excellent candidates for contributing to the 
dark matter density, in particular because we are certain that they exist, in 
contrast with other candidate particles. Together with CMB photons, relic 
neutrinos can be found anywhere in the Universe with a number density of 
339 neutrinos and antineutrinos per cm^. Depending on the value of the mass, 
this density is enhanced when neutrinos cluster into gravitational potential 
wells, although recent analyses show that the overdensity is limited to small 
factors (see e.g. [66]). 

In addition, it is easy to have a neutrino contribution of order unity to the 
present value of the energy density of the Universe, just by considering eV 
neutrino masses. In such a case, we saw before that all neutrinos should ap- 
proximately share the same mass mo and their energy density in units of the 
critical value of the energy density (see Eq. (21)) is^ 



where h is the present value of the Hubble parameter in units of 100 km s~^ 
Mpc~^ and J2inT-i = 3mo. Even if the three neutrinos are non-degenerate in 
mass, Eq. (12) can be safely applied. Indeed, we know from neutrino oscillation 
data that at least two of the neutrino states are non-relativistic today, since 
both (Am^i)^/^ ~ 0.047 eV and (AmiJ^/^ ~ 0.009 eV are larger than the 
temperature T,^ ~ 1.96 K ~ 1.7 x 10~^ eV. If the third neutrino state is very 
light and still relativistic, its relative contribution to is negligible and Eq. 
(12) remains an excellent approximation of the total density. 

If we demand that neutrinos should not be heavy enough to overdose the 
Universe (Qj, < 1), we obtain an upper bound mo < 15 eV for the absolute 
neutrino mass scale mo = Sj^i/S (fixing h = 0.7). This argument was used 
many years ago by Gershtein and Zeldovich [67] (see also [68]). Since from 
present analysis of cosmological data we know that the approximate contri- 
bution of matter is — 0.3, the neutrino masses should obey the stronger 
bound mo < 5 eV. 

Dark matter particles with a large velocity dispersion such as that of neutrinos 
are called hot dark matter (HDM). The role of neutrinos as HDM particles 
has been widely discussed since the 1970s, and the reader can find a historical 
review in Ref. [69]. It was reahzed in the mid-1980s (see e.g. [70,71,72]) that 

^ For high precision, the relation between (mi, m2, ma) and Q^h^ must be evaluated 
numerically using the distorted distributions described in [47] . 




93.14 /i2 eV ' 
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HDM affects tlic evolution of cosmological perturbations in a particular way: 
it erases the density contrasts on wavelengths smaller than a mass-dependent 
free-streaming scale (we will discuss in more detail the effects of massive neu- 
trinos on the evolution of cosmological perturbations in Sees. 4.5 and 4.6). In a 
universe dominated by HDM, this suppression is in contradiction with various 
observations. For instance, large objects such as superclusters of galaxies form 
first, while smaller structures like clusters and galaxies form via a fragmen- 
tation process. This top-down scenario is at odds with the fact that galaxies 
seem older than clusters. 

Given the failure of HDM-dominated scenarios, the attention then turned to 
cold dark matter (CDM) candidates, i.e. particles which were non-relativistic 
at the epoch when the universe became matter-dominated, which provided a 
better agreement with observations. Still in the mid-1990s it appeared that 
a small mixture of HDM in a universe dominated by CDM fitted better the 
observational data on density fluctuations at small scales than a pure CDM 
model (see e.g. [73]). However, within the presently favoured ACDM model 
dominated at late times by a cosmological constant (or some form of dark 
energy) there is no need for a significant contribution of HDM. Instead, one can 
use the available cosmological data to find how large the neutrino contribution 
can be. In Sec. 4, we will explain the effect of neutrino masses on cosmological 
observables. In Sec. 5, we will review the upper bounds on the sum of all 
neutrino masses which have been derived from current data. 



4 Massive neutrinos and cosmological perturbations 

Let us start with some generalities on the theory of cosmological perturba- 
tions. This field has been thoroughly investigated over the past thirty years, 
and many excellent reviews have been written on the subject (see for instance 
Refs. [42,74,75,76,77]). Many equations in this section are reminiscent of those 
in Ma & Bertschinger [74] and Bertschinger [77] - however, with several sign 
differences, since we choose the metric signature to be (+, — , — , — ). This sig- 
nature tends to be the most popular nowadays in cosmology. 

In our Universe, the metric and the energy-momentum tensor are inhomoge- 
neous. Their perturbations, given by 



are known to be small in the early Universe, typically 10^ times smaller than 
the background quantities, as shown by CMB anisotropies. As we shall see in 



%i/(x, t) = ^^^(x, t) - g^^{t) , 



(13) 
(14) 
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this section, after photon decouphng, the matter pcrtTirbations grow by grav- 
itational collapse and reach the non-linear regime, starting with the smallest 
scales. However, the linear perturbation theory is a good tool both for describ- 
ing the early Universe at any scales, and the recent universe on the largest 
scales. The most reliable observations in cosmology are those involving mainly 
linear (or quasi-linear) perturbations. In particular, current cosmological neu- 
trino mass bounds arc based on such observations. Therefore, our goal in this 
section is to describe the evolution of linear cosmological perturbations. The 
great advantage of linear theory is, as usual, to obtain independent equations 
of evolution for each Fourier mode. Note that the Fourier decomposition must 
be performed with respect to the comoving coordinate system: so, the quantity 
{27T/k) is the comoving wavelength of a perturbation of wavevector k, while 
the physical wavelength is given by 

A-a(t)^, (15) 



where a{t) is the scale factor of the Universe. For each mode k, the perturba- 
tion amplitudes evolve under some equations of motion (which depend only on 
the modulus k, since the background is isotropic), and on top of this evolution, 
the physical wavelength is stretched according to the Universe expansion. 

For concision, we will restrict ourselves to the main-stream standard cosmo- 
logical model: the inflationary ACDM scenario, whose background evolution 
is reviewed in Sec. 4.2. In Sec. 5 we will briefly comment on some results 
based on more exotic scenarios. However, the theory of cosmological pertur- 
bations is so rich that in order to present a short and pedagogical summary, it 
is necessary to start from strong assumptions concerning the underlying cos- 
mological background. Therefore, we will assume that at early times (e.g. at 
BBN) the Universe contains a stochastic background of Gaussian, adiabatic ^ 
and nearly scale-invariant primordial perturbations, as predicted by inflation, 
and as strongly suggested by CMB observations. Since the perturbations are 
stochastic and Gaussian, what we will call "amplitude" should often be un- 
derstood as "variance" : the linear equations of evolution give the evolution of 
each individual mode (5>/k of a perturbation 5f, and also, more interestingly, 
that of the root mean square {\Sf\'^)]/'^, obtained by averaging over all modes 
of flxcd wavenumber k. It is this last quantity which is the target of most 
observations, and which carries information on the cosmological parameters. 

The perturbations deflned in Eqs. (13) contain many degrees of freedom: the 

^ If we relax the assumption of adiabaticity, i.e. if we allow for isocurvature modes 
in the early Universe, the neutrino background can play a very specific role and 
provides some non-trivial initial conditions for the evolution of cosmological per- 
turbations. In this review, we will not cover this possibility and refer the interested 
reader to Refs. [78,79,80,81,82,83,84]. 
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homogeneity and isotropy of the background imphes that cji^^^ and T^j^ are 
diagonal, but in general this is not true at the level of perturbations. However, 
we will see in Sec. 4.3 that some of these degrees of freedom are just artifacts 
of the relativistic perturbation theory set-up; moreover, only a fraction of 
the physical degrees of freedom contribute significantly to the CMB and LSS 
observables. So, wc will see at the end of Sec. 4.3 that the problem can be 
reduced to the integration of a small number of linear equations of evolution. 

In order to perform this integration, one must specify the properties of each 
fluid contributing to the energy density. Apart from the cosmological constant, 
the standard ACDM scenario includes contributions from photons, neutrinos, 
baryons and cold dark matter. In Sec. 4.4, we will review the behavior of the 
most relevant perturbations under the assumption that neutrinos are absent, 
in order to better identify the specific role of neutrinos in Sec. 4.5. In Sec. 4.6, 
we will summarize the effects of neutrino masses on cosmological observables, 
which are important for understanding the present bounds. Thus, the reader 
already familiar with cosmological perturbation theory in ACDM and AMDM 
(A Mixed Dark Matter) models can go directly to the last subsection. 

4-1 Observables targets: definition of the power spectra 

Our goal is to compute observable quantities like 

• the CMB temperature anisotropy power spectrum, defined as the angular 
two-point correlation function of CMB maps ST/T{h) {h being a direction 
in the sky). This function is usually expanded in Legendre multipoles 

\^jr{n)jrin)j = 2^ CiPi{n-n), (16) 

where Pi{x) are the Legendre polynomials. So, for Gaussian fiuctuations, all 
the information is encoded in the multipoles Ci which probe correlations on 
angular scales 9 = n/l. Since for each neutrino family the mass is already 
known to be at most of the order of 1 eV (see Sec. 5), the transition to 
the non-relativistic regime is expected to take place after the time of re- 
combination between electrons and nucleons, i.e. after photon decoupling. 
The shape of the CMB spectrum is related mainly to the physical evolution 
before recombination. Therefore, the CMB will only be marginally affected 
by the neutrino mass (we will see later that there is only an indirect ef- 
fect through the modified background evolution). For this reason, we will 
only mention in this review some essential aspects of the CMB temperature 
spectrum, without entering into details. For the effects on CMB of heavier 
neutrinos (with masses of a few eV) , we refer the reader to Ref . [85] . 
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• the CMB polarization anisotropy power spectra bring interesting comple- 
mentary information to the temperature one, but for the same reasons as 
above, plus the fact that current polarization measurements are still rather 
imprecise, we will not discuss polarization in this section. 

• the matter power spectrum, observed with various techniques described in 
the next section (directly or indirectly, today or in the near past), probes 
the current Large Scale Structure of the Universe. It is defined as the two- 
point correlation function of non-relativistic matter fluctuations in Fourier 
space 



where 5^ — (^Pm/Pm- In the case of several fluids (e.g. CDM, baryons and 
non-relativistic neutrinos), the total matter perturbation can be expanded 
as 



Since the energy density is related to the mass density of non-relativistic 
matter through E = mc^, Sm represents indifferently the energy or mass 
power spectrum. When the redshift z is not explicitly written, we assume 
that P{k) refers to the matter power spectrum evaluated today (at z = 0). 
The shape of the matter power spectrum is the key observable for constrain- 
ing small neutrino masses with cosmological methods. Therefore the main 
purpose of this section is to review the physics responsible for its shape, 
in a pedagogical but sufficiently detailed way, starting with linear structure 
formation in absence of neutrinos (Sec. 4.4), and then including the impact 
of free-streaming massive neutrinos (Sees. 4.5, 4.6). 

• finally, one can build some cross-correlation power spectra (which measure 

the angular or spatial correlation between two different types of observ- 
ables). These correlations are in principle very interesting to extract from 
observations, since they provide some additional information beyond that 
contained in the self-correlation power spectra Ci or P{k). For instance, the 
cross-correlation between CMB temperature and polarization maps is non- 
zero, since the two signals are affected by the same physical inhomogeneities 
close to the last-scattering surface. The CMB temperature maps are also 
correlated to some extent to the matter distribution at low redshift, due to 
the late integrated Sachs- Wolfe effect (that we will briefly introduce in Sec. 
4.4.7). Finally, different measurements of the matter distribution in differ- 
ent overlapping redshift ranges or with different techniques are likely to be 
correlated with each other. The measurement of various cross-correlation 
power spectra is expected to be crucial in the future, but since nowadays it 
only plays a marginal role in neutrino mass determinations [86] , we will not 




(17) 



5, 



m 



Ei Pi 
Ei Pi 
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introduce the corresponding formalism in this review. 



4-2 Background evolution 



In most of this section, we will use the conformal time r instead of the 
proper time t. The two are related through dr — dt/a{t). A dot will de- 
note a derivative with respect to r. For instance, the Hubble parameter reads 
H — d{\a.a)/dt = a/a^.ln absence of spatial curvature, the background Pried- 
mann metric reduces to 

ds'^ = a(r)2[dr2 - 5ijdx'dx^\ , (19) 



and the Priedmann equation reads 

H = (Pt + Pcdm + Pb + Pl. + Pa) , (20) 



where the homogeneous density of photons scales like a~^, that of non- 
relativistic matter (pcdm for CDM and pb for baryons) like a~^, that of neu- 
trinos pi, interpolates between the two behaviors at the time of the non- 
relativistic transition, and finally the cosmological constant density p\ is of 
course time-independent. At any time, the critical density pc is defined as 
pc = 3H^/8ttG, and the current value Hq (or reduced value h) of the Hubble 
parameter gives the critical density today 

p° = 1.8788 X lO-^^/i^ g cm"^ . (21) 



The normalization of the scale factor a is arbitrary. If wc first neglect neu- 
trino masses, the photon and neutrino densities can be combined into a total 
radiation density pr = P7 + Pu scaling like a^^, while the CDM and baryons 
combine into a matter density pm = Pcdm + Pb scaling like a~^. Then, we can 
choose the normalization of a such that pm = 3/(87rGa^) in order to obtain 
a particularly simple Fricdmann equation valid during radiation domination 
(RD) and matter domination (MD) 

(d)' = a + Oeq , (22) 



where Oeq is the scale factor at equality, when pm = Pr- There is a simple 
solution, 

2 

a^^ + ^T^j(T + 2{V2 + l)T^) , (23) 
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Fig. 5. Evolution of the background densities from the time when = 1 MeV (soon 
after neutrino decouphng) until now, for each component of a flat AMDM model 
with h = 0.7 and current density fractions = 0.70, Ob = 0.05, Qi, = 0.0013 and 
^cdm = 1 — ^^A — — ^lu- The three neutrino masses are distributed according 
to the Normal Hierarchy scheme (see Sec. 2) with mi = 0, m2 = 0.009 eV and 
771-3 = 0.05 eV. On the left plot we show the densities to the power 1/4 (in eV units) 
as a function of the scale factor. On the right plot, we display the evolution of the 
density fractions (i.e., the densities in units of the critical density). We also show on 
the top axis the neutrino temperature (on the left in eV, and on the right in Kelvin 
units). The density of the neutrino mass states 1^2 and 1^3 is clearly enhanced once 
they become non-relativistic. On the left plot, we also display the characteristic 
times for the end of BBN and for photon decoupling or recombination. 



where r, 



cq 



2(V2- 



is the value of the conformal time at equality. At 



low redshift (typically z < 0.5), the cosmological constant density takes over, 
causing a departure from the above solution, with an acceleration of the scale 
factor. Finally, if we include the effect of small neutrino masses, the solution 
is also slightly modified, since the non-relativistic transition of each neutrino 
species amounts in converting a fraction of radiation into matter. This can 
be seen in Fig. 5, where we plot the evolution of background densities for a 
AMDM model in which the three neutrino masses follow the Normal Hierarchy 
scheme (see Sec. 2) with mi = 0, 1712 = 0.009 eV and = 0.05 eV. 



4-3 Gauge transformations and Einstein equations 



In the real Universe all physical quantities (densities, curvature...) are func- 
tions of time and space. Thanks to the covariance of general relativity, they 
can be described in principle in any coordinate system, without changing the 
physical predictions. The problem is that in order to obtain simple equa- 
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tions of evolution, wc wish to use a linear perturbation theory, in which the 
true physical quantities are artificially decomposed into a homogeneous back- 
ground and some small perturbations. This is artificial because the homoge- 
neous quantities are defined as spatial averages over hypersurfaces of simul- 
taneity: f{t) — (/(i,x))x. Any change of coordinate system which 

(1) mixes time and space (therefore, redefining hypersurfaces of simultaneity, 
and changing the way to perform spatial averages), and 

(2) remains small everywhere, so that the differences between true quantities 
and spatial averages are still small perturbations, 

gives a new set of perturbations (new equations of evolution, new initial con- 
ditions), although the physical quantities (i.e., the total ones) are the same. 
This ambiguity is called the gauge freedom in the context of relativistic per- 
turbation theory. 

Of course, using a linear perturbation theory is only possible when there exists 
at least one system of coordinates in which the Universe looks approximately 
homogeneous. We know that this is the case at least until the time of photon 
decoupling: in some reference frames, the CMB anisotropics do appear as small 
perturbations. It is a necessary condition for using linear theory to be in such 
a frame; however, this condition is vague and leaves a lot of gauge freedom, i.e. 
many possible ways to slice the spacetime into hypersurfaces of simultaneity. 

We can also notice that the definition of hypersurfaces of simultaneity is not 
ambiguous at small distances, as long as different observers can exchange light 
signals in order to synchronize their clocks. Intuitively, we see that the gauge 
freedom is an infrared problem, since on very large distances (larger than the 
Hubble distance) the word "simultaneous" does not have a clear meaning. The 
fact that the gauge ambiguity is only present on large scales emerges naturally 
from the mathematical framework describing gauge transformations. 

Formally, a gauge transformation is described by a quadrivector field e'^{x,t) 
(see e.g. Ref. [74]). When the latter is infinitesimal, the Lorentz scalars, vectors 
and tensors describing the perturbations are shifted by the Lie derivative along 

5V...(x,i) ^ SA^,,„{^,t) + Ce[SA^,,„{^,t)] . (24) 

Since there are four degrees of freedom (d.o.f.) in this transformation - the four 
components of - we see that among the ten d.o.f. of the perturbed Einstein 
equation 60^1, = SnGST^i,, four represent gauge modes, and six represent 
physical degrees of freedom. 
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In addition, it can be shown that this equation contains three decoupled sec- 
tors. In other words, when the metric and the energy-momentum tensor are 
parametrized in an adequate way, the ten equations can be decomposed into 
three systems independent of each other: 

(1) four equations relate four scalars in the perturbed metric Sgf^u to four 
scalars in ST/j^i,, 

(2) four equations relate two transverse 3-vectors in the perturbed metric (4 
d.o.f. in total) to two transverse 3-vectors in ST^^, 

(3) two equations relate one transverse traceless 3x3-tensor in the perturbed 
metric (2 d.o.f. in total) to a similar tensor in ST^,^. 

Moreover, this decomposition is left invariant by gauge transformations. These 
three types of variables are called scalar, vector and tensor modes; physically, 
they describe respectively the generalization of Newtonian gravity, gravito- 
magnetism, and gravitational waves. It can be shown that two scalar d.o.f. 
and two vector d.o.f. are gauge modes, so each of the three sectors contains 
only two physical d.o.f. In principle, all modes can contribute to the CMB 
anisotropy maps, but CMB observations show that the vector and tensor con- 
tributions are negligible (at least for temperature anisotropics). As for the 
LSS of the Universe, it is related to the mass/energy density distribution of 
non-relativistic matter, i.e., again to the scalar sector. Therefore, in this re- 
port, we will focus only on scalar perturbations, not because the cosmological 
backgrounds of vector /tensor modes are insensitive to neutrino properties, but 
because they are not likely to be observed with great precision in a near future. 

It is possible to build some gauge-invariant combinations, and to reduce the 
Einstein equation into a set of gauge-invariant equations. This is not the most 
economic way to proceed: one can simply choose an arbitrary gauge-fixing con- 
dition, i.e., a prescription that will limit the number of effective degrees of free- 
dom to that of physical modes only, and make all calculations into that gauge. 
When the same problem is studied in two different gauges, the solutions can 
look very different on large wavelength; for instance, the total energy density 
perturbation 5p{t, k) for a given time and wavenumber can appear as grow- 
ing in one gauge, and constant in another gauge, although the two solutions 
describe the same Universe. However, physical observables - like the matter 
density perturbations probed by galaxy surveys or temperature/polarization 
anisotropies probed by CMB experiments - are always limited to small scales, 
at most of the order of the Hubble length. On those scales, the predictions 
arising from different gauge choices always coincide with each other. 

Throughout our discussion, we choose to work in the longitudinal gauge, which 
is probably the most popular for studying cosmological perturbations. In this 
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gauge, one requires that the non-diagonal metric perturbations vanish. This 
ehminates two scalar degrees of freedom, and the remaining ones are defined 
in such a way that Eq. (19) reads 

ds^ = g^.^dx^dx" = a'^(r)[(l + 20)(iT^ - (1 - 2'ijj)5ijdx'dx^] . (25) 

For the energy-momentum tensor, the four scalar degrees of freedom can be 
identified as 



5T^^5p, (26) 
6T^^{p + p)v^, (27) 

5r; = -5p5i + E;ii , (28) 

where 6p is the energy density perturbation, 6p the pressure perturbation, = 
diV the longitudinal component of the velocity field, and E*" = {didj — ^6ij'V^)a 



the traceless and longitudinal-divergence component of the 3x3 tensor ST^. 
For these last two degrees of freedom, we could write the equations in terms 
of the two potentials v and a. However, it is more conventional to deal with 
the two quantities {9, a) which represent respectively the velocity divergence 
and the shear stress (or anisotropic stress) 

e = ^diVi^V^v , (29) 

i 

(p+p)vV^-^(a,a, - ^v%)e;.'i = ~v'a . (so) 

Let us recall that the four scalar components of the energy-momentum tensor 
are not gauge-invariant: we define the above quantities in the longitudinal 
gauge, but in another gauge, for instance, 5Tq could behave differently than 
the above Sp, especially on super- Hubble scales. 

We can now write the perturbed Einstein equation for the scalar sector in the 
longitudinal gauge. 



5G^o = 2a-2 <^ -3 



3-ip + VV } = SttG Sp , 
a 



5G^i = 2a-'^di ^-(j) + i)\ = St^G {p + p) Vi , 



5G) = -2a- 



/ a 




2-- 





\ a 





(/. + -(</. + 2V;) +V; + ^V^(</>-^) 



(0-V')} 



(31) 
(32) 

(33) 
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= 87rG (-5p5i + E}) . (34) 

Using the variables 6 = Sp/p, 6p, 9, a, and switching to comoving Fourier 
space, we obtain 



(• . 2 
-) (P-3-iP-k^i^^A7rG a"" p5 , (35) 
aj a 

(^^(j) + ip^= iTTG {p + p)e , (36) 

{^l - (^) ') <^ + ^l*^ + 2V') + - y (0 - V') = 47rG a' 5p , (37) 

k'^(cf) -ip) = 127iG c? {p + p) a .(38) 

Through the Bianchi identities, the Einstein equation imphes the conservation 
of the total energy- momentum tensor. Actually, the energy-momentum tensor 
of each uncoupled fluid is conserved, and obeys the continuity equation 

5 ^ {1 + w){e + 3ij) , (39) 



and the Euler equation 

9 = ^{3w - 1)0 - -^9 - e<f> - k'^a - ^^k'^S , (40) 



where we assumed an equation of state p[x, t) = wp{x, t) for the fluid, so that 
w = p/p = 6p/6p. 



4-4 Linear perturbation theory in a neutrinoless Universe (pure ACDM) 



4-4-1 Perfect fluids 

Intuitively, in perfect fluids, microscopic interactions enforce a "collective be- 
havior" of the particles. More precisely, they guarantee that the stress tensor 
Tij is isotropic and diagonal, and that the energy-momentum tensor can be 
simply described in terms of functions of time and space: the density p, the 
pressure p and the bulk velocity U^, 

TM- ^ _p gi^- + + p)U''U'' . (41) 



Here, = dx^ / \J ds"^ is the 4- velocity. As long as the bulk 3- velocity ~ 
dx^ /dr remains first-order in perturbations, we can approximate the 4- velocity 
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by = dx^/lail + (t))dT] and obtain: 



C/'^=(a-i[l-0],a-V) , r°-p, T^ = v\ ^ -p , (42) 



where T/ is no^ summed over the repeated index. In a pure ACDM model, the 
Universe contains photons, baryons and cold dark matter. Strictly speaking, 
the perfect fluid approximation does not apply to this situation, since on the 
one hand photons become collisionlcss after decoupling, and on the other hand 
cold dark matter is coUisionless at any time (so it cannot even be called a fluid) . 
However, we will see in the next paragraphs that we can use Eqs. (41, 42) in 
an effective way. 

First, after its own decoupling time, the cold dark matter component is always 
non-relativistic and coUisionless. As long as the bulk motion of CDM particles 
can be described by a single- valued flow, the CDM energy-momentum tensor is 
formally identical to that of a perfect fluid with density pm and zero pressure. 

Second, let us discuss the case of photons and baryons. If our purpose was 
to review the CMB physics, it would be crucial to introduce separately the 
exact form of the photon and baryon energy-momentum tensors. However, 
in this report, we only wish to provide a simplified description of large scale 
structure formation. For this purpose, we can view the baryon and photon 
components before recombination as forming a single tightly-coupled fluid, 
obeying to Eq. (41) with the equation of state of ordinary radiation: pj. = 
(this assumes that the baryon density is negligible with respect to the photon 
density, which is only a crude approximation around the time of equality and 
recombination). The local value of pr is related to the local value of the photon 
blackbody temperature. The bulk velocity of this fluid remains small enough 
to use Eqs. (42). After recombination, the photons are irrelevant for structure 
formation, and we will simply discard this component; the baryons behave like 
a non-rclativistic collisionlcss medium, but for simplicity we will neglect their 
density with respect to that of cold dark matter. 

Finally, since we are only interested in scalar perturbations, we can reduce the 
bulk velocities of radiation and CDM into irrotational vector fields, and define 
two scalars 9^ and 9^ as in Eq. (29). 

In summary, we can solve the Einstein equations with the following perturbed 
energy-momentum tensor 



d\5T^) = (pr + Pr)9r + Pra9^ = ^pA + P^9, 



(43) 
(44) 
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5T^^-5p,^~6p,. (45) 



4.4-2 Jeans length 

In a spatially flat Priedmann Universe, a physical process switched on at a 
time ti and propagating at a velocity v along radial geodesies {v dt — a{t) dx) 
can only affect wavelengths smaller than a causal horizon, defined as 



t 



/r y df' 
dx = a{t) J — - . (46) 



a{t' 



This horizon is simply the maximal physical distance on which the signal can 
propagate between ti and t (sometimes, is it defined with a factor two bigger 
than above). The so-called particle horizon dn obeys this definition in the 
particular case of a signal traveling at the speed of light, v — 1 (we adopted 
units in which c = 1). If both ti and t are chosen during the matter or radiation 
dominated stage, and if t ^ ti, it is straightforward to show that the particle 
horizon can be approximated by the Hubble length Rnit) = 1/H{t), up to 
a numerical factor of order one. Indeed, for a power-law expansion a{t) oc 
with n < 1, one has 

Rnit) = - , dnit > ti) ~ . (47) 



Similarly, the characteristic velocity under which acoustic perturbations 
propagate before photon decoupling defines a characteristic length ds{ti,t), 
called the sound horizon. For a constant sound speed, and under the conditions 
described above for ti and t, the sound horizon is obviously given (up to a 
numerical factor of order one) by the ratio Cs/H{t), which is called the Jeans 
length. The exact definition of the Jeans wavenumber and Jeans length is 
usually 



where the numerical factors -which are unimportant for the purpose of un- 
derstanding the physics- are simply dictated by the particular form of the 
Newtonian equation of evolution of the density contrast of a single uncoupled 
perfect fluid with constant sound speed, 

5 + -s + (e -e^cU^o , (49) 

a 
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which derives from the combination of the continuity, Euler and Poisson equa- 
tions. 

Our physical expectation is that for a fluid of sound speed Cg, modes with 
k > kj will oscillate with a pulsation lo = kCg, due to the competition between 
the gas pressure and the gravitational compression. These modes are said 
to be Jeans stable. On the other hand, for modes k < kj, pressure cannot 
causally resist to gravitational compression, and density perturbations should 
in principle grow monotonically. This Jeans instability explains some basic 
phenomena in the inhomogeneous Universe: before recombination, the photon- 
baryon fluid has a sound speed of order Cg ~ and oscillates on scales 

smaller than Aj; after recombination, Cg and Aj become vanishingly small, kj 
grows to infinity and structures form. 



4.4-3 Radiation domination 

Deep inside the radiation era, we have seen that the radiation component can 
be approximated by a perfect fluid of sound speed c^, due to the tight coupling 
between baryon and photons. This self-gravitating fluid oscillates inside the 
Jeans length (or sound horizon). The cold dark matter evolves as a test fluid 
in this background, until the time (close to equality) at which its gravitational 
back-reaction becomes important. However, for simplicity, we will work at first 
order in the expansion parameter Pm/pr- 

In this limit, it is possible to derive an analytic expression for the evolution 
of (5r and 5m- The absence of anisotropic stress implies that the two metric 
fluctuations are equal, (j) — il^. The full continuity, Euler and Einstein equations 
for the radiation and matter fluid read 



<^m = ^m + 30, + (50) 

L^--Om-k''(f> , er^-'^5,-k^cl>, (51) 
a 4 

-3^0 - (3 ' + j = AnGa^ {p^S^ + pA) , (52) 

-A;' (^0+^0) = inGa^(^Pm0ra+^Pr9?j , (53) 
a ■ ( 'd / d\'^\ , . „ o / 1 



+ 3-0+ (^2- - (^-j J = 47rG«' [■^PA) , (54) 

and only five out of these seven equations are independent. In addition, one 
out of the five independent equations is just a constraint equation, so the 
system admits four independent solutions. In order to find them, let us first 
notice that since we have two fluids, we can have differences in the number 
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density contrasts, usually parametrized by the entropy perturbation 



3 . 



(55) 



The pressure perturbation 5p = 5pr can be uniquely decomposed in the basis 
of the total density perturbation 5p = 5p-c + 5pui and entropy perturbation r;, 



5p 



3 1 + 



3prr 
4pr 



H -1 



{5p + p^vi) . 



(56) 



This defines the effective sound speed in the two-fiuid as 

-1 



-eflf 



3 1 + 



3pm 



P 
P 



(57) 



We can reduce the system (50)-(54) to a pair of second-order coupled equa- 
tions, obtained from the combinations [(54) — (52)] and [A;^(52) — 3d/a(53)], 
as follows 



a.. 2 - 1 



(/) + 3-(l + cfff)0+- 
a a 



3 l + f cl,-^^ 



P} 
Pm 



<l> + k'cl^<P = -clsV, (58) 



3r 



—r] + ——V = — «— -(59) 



eff 



a 



4 Pr 



6 Pr 



Once the four independent solutions have been found, the density contrasts 
can be obtained e.g. from combining Eqs. (54) and (58) 



(60) 
(61) 



and the velocity gradients follow from Eqs. (51). The evolution of the scale 
factor can be taken from equation (23). Deep inside the radiation era, we 
can write a, Cgg and the two equations [r^ x (58)] and (59) at first order in 

Pm/Pr = r/Teq, 



+ 4r0 + (T^a;^)^ = , f) + -?) = ^^^^0 , 

T Z 



(62) 



where uj = kces = k/\/3. In this limit, the first equation decouples: this just 
confirms that the radiation fiuid is self-gravitating, and the CDM behaves as 
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a test fluid. Two independent solutions of the system, both satisfying r){T 
0) = 0, read 



(f) = (a;T)"^ ~ cosa;T^ + D{k) ^-sinu;T 

77 = -3|c(^) (/ 



cos CUT ' 



UJT 



COS X — 1 1 . . 

ax — -{cos CUT — 1) 



j 



smx 1 

dx + - sm u)T 

X 2 



(63) 



The solution proportional to C{k) represents the growing adiabatic (or isen- 
tropic) mode, the other one is the decaying adiabatic mode. There are two 
other obvious solutions (0,7^) = (0, ^(A;)) and (0,7^) = (0, S(A;)lnT) which 

stand respectively for the isocurvature growing and decaying mode. For the 
growing adiabatic mode, we get from Eqs. (60) and (61) that the radiation 
and matter density contrasts evolve like 



, f 2 sin a;r — a;r cos a;r sm ujt — ujt cos ujt\ 
6^^-2C[k)\ 2 — \ , 

1 sm.uJT sin UJT — UT COS UJT 1 /"cosa;— 1 , , 

5^^?,C{k) { + ——: + o + / dx ) . 64 

' UJT [ujtY 2 J a; ' 

The integral is given by 

I COSX ^ ^.^^^^ _ ^ _ ^^^^^^ ^ ^g^^ 



where Ci(a;) is the cosine integral and C ~ 0.5772 is the Euler constant. For 
a; > 1, the function Ci(a;) describes damped oscillations around zero. Let 
us summarize the physics described by these solutions. During the radiation 
era and far outside the sound horizon, all quantities appear as frozen in the 
longitudinal gauge, with 

5^ = = ~(t> = -\c{k) . [RD, UT = kT/^/?, < 1] (66) 

Then, inside the sound horizon, the evolution of the radiation and matter 
perturbations are radically different from each other. The self-gravitating ra- 
diation component oscillates on scales smaller than its Jeans length, and Eq. 
(64) represents a first-order approximation to the acoustic oscillations ob- 
served in the CMB anisotropy spectrum. Instead, the cold dark matter has 
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a considerably smaller Jeans length and does not fluctuate on its own. The 
gravitational driving force only induces damped oscillations on S^^, on top of a 
net gravitational clustering effect. This clustering is not very efficient, because 
the gravitational potential fluctuations are damped by the acoustic oscillations 
and do not feel signiflcant back- reaction from the CDM component. In aver- 
age, the growth of the CDM perturbations is proportional to the logarithm of 
conformal time 

5m ~ 3C(fc)[l/2 + C + ln(a;T)] . [RD, ur ^ kr / Vs I] (67) 



4 -4 -4 Matter domination 

Deep into the matter domination era, when a cx r^, the effective sound speed 
is vanishingly small and Eq. (58) simply reduces to 

+ 3-0=0, (68) 



with solutions 

(j) = C{k)+T-^D{k) . (69) 
The matter density contrast then follows from Eq. (52), which simplifles into 

-(■ 



C4 + kA<p-^^= ^S^ , (70) 



T 



and we flnally obtain 

= - (4+ ^ I C{k) + I 11— - 6 ) T-'D{k) . (71) 



We can identify the solution proportional to C{k) (resp. D{k)) with the adia- 
batic growing (resp. decaying) mode. A detailed and precise matching between 
the solutions (66), (67) -valid deep inside the radiation era- and the solution 
(71) -valid deep inside the matter era- would require some significant work in 
order to take into account various effects close to the time of equality. How- 
ever, a crude matching is sufficient for catching the essential behavior of the 
matter power spectrum. Prom (66), (67) and (71) we obtain that at any time 
T belonging to the MD stage, 

M..x)=(4.^)i^. (72) 
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where 5m(^,Teq) IS equal to Eq. (64) evaluated at r = Tcq. Since we are inter- 
ested only in modes which are inside the Hubble radius in the recent Universe, 
we can take the limit /c^r^ ^ 1. During matter domination, it is also possible 
to replace r by 2/ {aH). Then, we obtain 



4.4-5 Dark energy domination 

After the matter dominated era, the dark energy component becomes the 
dominant contribution to the energy density of the Universe. In this era, the 
evolution of the gravitational potential is still given by Eq. (54), but now the 

effective mass term does not cancel. When the dark energy corresponds to 
the contribution of a cosmological constant with equation of state parameter 
w = -1 (Pa = -Pa), 

+ 3-0 + (SnGa^pA) = 0, (74) 



and the potential decays proportionally to a scale-independent damping factor 
g{T) — (I){t) / (f){TM) , where tm can be any time during matter domination. This 
factor can be found by numerical integration. In the case of a cosmological 
constant, a fairly good analytic approximation of g{To) - where tq is the present 
time - is provided in Ref. [87] 

g{To) ~ n'^yil + 0.003(fiA/^^m)'/'] . (75) 



The modes probed by LSS experiments are all deep inside the Hubble ra- 
dius during the end of matter domination and during A domination. On 
those scales, the Poisson equation of Newtonian gravity (which corresponds 
to Eq. (52) with k ^ aH) gives a simple relation between the gravitational 
fluctuations and the matter perturbations. 



So, is proportional to a(f). which proves that the decay of the gravitational 
potential slows down the growth of matter fluctuations as a function of the 
scale factor. Also, we immediately obtain a relation between S^i during A 
domination and matter domination 



a{TM)H{TM) 



a{T)H{T) 



2 



^^^^ -<5m(TM) . (77) 
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Finally, Eqs. (64), (72) and (77) give us the expression of matter perturbations 
on sub-Hubble scales during A domination. 



9{r)k'C{k) 



X 



(sin kx) I ikx) + (sin kx — kx cos kx) / (kx)^ + i + Ci{kx) 



(78) 



C - la{kx) 



[1 + (kxy/12] 



where X = Tcq/^/S. The matter power spectrum is defined as P{k,T) — 
{\5ra{k,r)\'^). The previous equation allows us to relate it to the primordial 
power spectrum, usually parametrized by a power-law 

k'{\<P{k,T^ 0)1^) = ^{\Cik)\') ^ Ak-'-' , (79) 



where the index Ug is the scalar tilt. We can compute the two asymptotes 
for the modes which were far outside/inside the sound horizon at the time of 
equality. 



kr^/Vs < 1 , (80) 

kr^/Vs > 1. (81) 

Because of the crude approximations performed around the time of equality, 
this solution is not accurate, in particular for the global amplitude and for 
the behavior near k ~ y/S/req. However, the shape of the two asymptotes 
reflects the behavior of the true power spectrum obtained by a full numerical 
simulation, which is close to the BBKS [88] fitting formula 



^^^'^^'^n^(r)2[a(r)i/(r)]4 4 ' 

gjry 9[12YAk-^[lnkxY 



P(k r) - ^(^)' ^ Tia? (82) 

^^'''^^-nUrmrWr)]^ 25 ^^^^ ' ^^^^ 

T{q) ^ + 3.89 q + (16.1 qf + (5.46 qf + (6.71 q)']-'^' , 

k 

q 



n^h'^ exp{-2Qb) Mpc ' 



Note also that the cosmological constant not only damps the perturbations 
5pm but also implies a smaller fraction fim- The net efi'ect is that today and for 
fixed aoHo, the power spectrum normalization is enhanced by the cosmological 
constant (i.e. g{To)/Q.ra{To) > 1 if Qa > 0). 
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Fig. 6. Evolution of the longitudinal metric perturbations (j) = ^ (top) and the mat- 
ter density perturbations 5m = Spm/pm (bottom) in a neutrinoless ACDM model, 
obtained numerically with a Boltzmann code, as a function of the scale factor a and 
Fourier wavenumber k. In order to get a better view, we chose opposite time axes in 
the two plots: time is evolving from back to front in the upper plot, and from front 
to back in the lower plot. The initial condition was set arbitrarily to k^^'^cf) = — 10~^. 

4.4-6 Numerical results 

We show in Fig. 6 the evolution of = (p and 6^ for a neutrinoless ACDM 
model, obtained numerically by solving the exact Einstein and conservation 
equations. We run the public code cmbfast [49] based on the synchronous 
gauge, and convert the output into the longitudinal gauge variables used 
throughout this section. We start from a flat primordial spectrum with ampli- 
tude k^/'^cp = 10~^ (let us recall that the "amplitudes" calculated by the code 
should be interpreted in the real Universe as the variance of random Gaussian 
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Fig. 7. (Red solid lines) The power spectrum of CMB temperature anisotropies 
in multipole space (left), and the LSS matter power spectrum evaluated today in 
Fourier space (right), for a reference neutrinoless ACDM model with uj^ = 0.13, 
Oa = 0.74, = 0.022, Us = 0.95, r = 0.08. (Dashed and dotted lines) Same 
power spectra for three models in which we keep all parameters fixed but one: 
either Wm = 0.12, or Q\ = 0.80, or = 0.026. The changes in the power spectra 
illustrate the effects labeled from (1) to (4) in Sec. 4.4.7, as explained in details at 
the end of this subsection. 

perturbations). We show in Fig. 7 the CMB temperature anisotropy and LSS 
power spectra for the same cosmological model (but now with the primordial 
spectrum amplitude preferred by current data) 

The global evolution is consistent with the analytical solutions found in the 
previous subsections. Indeed, for metric perturbations ip = (p (upper plot, 
time evolving from back to front), we see that during radiation domination, 
super- Hubble modes are constant, and sub- Hubble modes are damped (with 
small oscillations). During matter domination, all modes are constant, so the 
shape imprinted near the time of equality survives until today (at least, as 
long as perturbations remain in the linear regime). During A domination, 
there is a small scale-independent damping of all perturbations, hardly seen 
on the figure. For matter density perturbations 5m (lower plot, time evolving 
from front to back), we see that during radiation domination, super-Hubble 
modes are constant, and sub- Hubble modes grow slowly (asymptotically like 
the logarithm of a). During matter domination, super-Hubble modes are still 
constant, while sub-Hubble modes grow linearly with a. The smallest scales 
are the first ones to reach the condition of non-linearity (5m ~ 1, and large 
scales are still linear today. During A domination, 5m grows a bit more slowly 
as a function of a. 



4-4-7 Parameter dependence 

A minimal flat ACDM models without neutrinos can be described with six 
parameters: the cosmological constant fraction ^a, the total non-relativistic 
matter density tUm = f^m^^, the baryon density ujb = Qhh"^, the primordial 
spectrum amplitude A and tilt n^, the optical depth to reionization r. These 
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parameters control various physical effect which are responsible for the shape 
of the observable power spectra: 

(1) The time of Radiation/Matter equality. In the above parameter 
basis, the time of equality between pm and pr is fixed by cUm only, since 
(aeq/ao) = Pr/Pm i^^^ sciipt ° means "evaluated today"), and p° is fixed 
by the CMB temperature. A late equality implies more metric pertur- 
bation damping for modes entering the Hubble radius during radiation 
domination (see Fig. 6). For the LSS matter power spectrum, the con- 
sequence is the following: since on sub-Hubble scales S^a grows more ef- 
ficiently during MD than during RD (because ip does not decay), it is 
clear that when Um decreases (late equality) the matter power spectrum 
is suppressed on small scales relatively to large scales. Also, the global 
normalization increases, because of the factor uj~^ in Eqs. (80,81). So, 
the net effect of postponing equahty is to amplifies P{k) only for small 
k. For the CMB, the effect is opposite and counter-intuitive: small-scale 
perturbations are boosted. This is related to the dilation effect before 
recombination, and to the early Integrated Sachs- Wolfe effect just after 
recombination. In simple words, the rapid decay of ip{k, r) just after Hub- 
ble crossing during RD tends to boost the beginning of the first acoustic 
oscillation through gravitational redshift effects. So, a later equality in- 
duces higher CMB peaks, especially for the first one. 

(2) The time of Matter/A equality. If the cosmological constant is larger, 
equality between matter and A takes place earlier. We have seen that af- 
ter that time, ijj decays and grows more slowly. Thus the matter power 
spectrum normalization is slightly suppressed by a long A domination. 
The CMB spectrum is affected in a more subtle way: the time-variation 
of the metric perturbations leads to a net redshifting of the CMB photons 
traveling across a gravitational potential well at small z (late Integrated 
Sachs- Wolfe effect). Therefore, the nearby distribution of galaxy cluster 
leaves an imprint in the large-scale CMB spectrum: the small-/ multipoles 
are enhanced (the larger is A, the more significant is the enhancement). 
This effect is difficult to measure experimentally due to cosmic variance 
(on large scales, we see only a few realization of the stochastic mode am- 
plitudes, and we expect a large scattering of the data points around the 
theoretically predicted variance). 

(3) The physical scale of the sound horizon at equality, and the 
angulcir scale of the sound horizon at recombination. From the 
previous sections it is clear that the characteristic scale in the matter 
power spectrum is fixed by the sound horizon at the time of equality 
(i.e., essentially, the Hubble scale /^^(req) divided by VS). Since this 
scale can be computed by extrapolating the present value of the Hub- 
ble parameter back in time, it depends on h, fl\, and on the time of 
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equality. In our parameter basis, it means that changing any of the two 
parameters cjm, Q\ shifts the matter power spectrum horizontally. The 
characteristic scale of the oscillations in the CMB power spectrum is set 
by the sound horizon (is(Trec) at recombination. The time Tree is more or 
less fixed by thermodynamics, but the sound horizon is an integral over 
Csdt/a{t) between and Tree- This integral depends on the time of equal- 
ity, and on the baryon density at late times (through Cg). In addition, the 
actual observable quantity is the angular scale of this sound horizon on 
the last scattering surface, which is given, for fixed cis(Trec), by h and fl\. 
We conclude that the observed angular scale of the peaks constrains a 
combination of the three parameters JIa, cuui and cub. 

(4) The balance between gravity and pressure in the tightly-coupled 
photon-beiryon fluid. In this section, for simplicity, we assumed that 
before recombination the energy density of the tightly-coupled photon- 

baryon fluid was dominated by that of photons. Actually, when the baryon 
density becomes significant (near the time of equality), the fiuid is more 
affected by the gravitational compression -which is more efficient for 
non-relativistic baryons than for relativistic photons- relatively to the 
competing photon pressure. Inside a gravitational potential well of given 
amplitude, the zero-point of oscillations (i.e. the instantaneous value of 
6T/T for which gravity and pressure exactly cancel each other) is dis- 
placed. The main effect in the CMB spectrum is an increase of the first 
and third peaks, but not of the second one. As far as the matter power 
spectrum is concerned, we should remember that P{k) probes the total 
contribution from baryonic and cold dark matter perturbations. In the 
limit of a small CDM contribution {ujh/uj„^ ^1); the small-scale matter 
power spectrum has a low amplitude, because it docs not benefit from 
the slow growth of CDM perturbations before equality on sub-Hubble 
scales (visible in Fig. 6). Also, in this limit, the oscillatory structure of 6b 
at the time of recombination survives until now, and the matter power 
spectrum exhibits visible wiggles. In the opposite limit ujb/uj^-a 0, which 
we adopted in our analytical approximate solution, the small-scale power 
spectrum is higher, and the baryonic oscillations in P{k) are small (still, 
they are visible in Fig. 7 between k — 0.05 Mpc and k — 0.2 
Mpc, and they have actually been observed recently, see Ref. [89]). 

(5) The amplitude of primordial perturbations (A), defined in Eq. (79), 
obviously fixes the global normalization of all spectra. 

(6) The tilt of primordial perturbations (n^), defined in Eq. (79), fixes 

the balance between large and small scale amplitudes. Together with 
the previous parameter it provides interesting constraints on infiationary 
models. 
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(7) The optical depth to reionization (r) is a rather simplified way of 
parameterizing the reionization of the Universe at redshift of order 10 
to 20 (depending on the model). This parameter has no effect on the 
matter power spectrum, but it damps the CMB temperature fluctuations 
for scales which are inside the Hubble horizon at the time of reionization, 
leaving small-Z temperature multipoles unchanged. 

In Fig. 7, we illustrate graphically the consequence of varying either a;m, 
or ujy^, while keeping all the other cosmological parameters fixed. When we 
decrease uj^, the relevant effects arc (1) and (3): equality is postponed, boost- 
ing the first CMB peak as well as P{k) for small k; simultaneously, there is 
a tiny horizontal shift of both spectra. When we increase Q^, we can see the 
effects (2) and (3): the smallest CMB temperature multipoles increase, the 
P{k) normalization decreases and both spectra are shifted to the left. When 
oui, is increased, the effects (3) and (4) are taking place: the first CMB peak 
increases, while the second one decreases slightly; P{k) is slightly suppressed 
on small scales; finally, the scale of CMB peaks is slightly modified. 

We have seen that {A, Ug, r) have very specific effects, while (Qa, cu^n, (^b) have 
intricate effects. Nevertheless, since there arc seven effects for six parameters, 
it is in principle easy to measure all of them from CMB and LSS data, assuming 
a flat ACDM model. However, the neutrinoless model described in this section 
for pedagogical purposes is unable to provide a good flt to the data: it is time 
now to include neutrinos. 

4-5 Linear perturbation theory in presence of neutrinos (AMDM) 

We will now restore the neutrino component, which is expected to behave 
at late times (after the non-relativistic transition) as an extra sub-dominant 
dark matter component: the standard cosmological scenario should actually 
be called A Mixed Dark Matter (AMDM). 

The impact of neutrinos on cosmological perturbations was flrst clearly ex- 
plained in 1980 by Bond, Efstathiou and Silk [70]. Later, during the 1980's, 
a very large number of papers investigated this issue in more details, and 
it would be difficult to provide an exhaustive list of references (see e.g. those 
given in [6,90]). Despite their interest, many of these papers are not of practical 
use nowadays because they generally assumed a neutrino-dominated universe 
with no CDM or A components. For instance, the interested reader can take 
a look at Refs. [71,72] and go through the historical summary of Ref. [69]. 

Let us flrst introduce some notations and formalism -which are very close 
to those of Ma & Bertschinger [74], but have been adapted to our metric 
signature (-I-, — , — , — ). The momentum of individual neutrinos can be labeled 
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by Pi, which is the canonical conjugate of the comoving coordinate x\ It differs 
from the proper momentum pj measured by a comoving observer (with fixed 
a;* coordinates) by 

Pi = a(l - . (83) 



Note that in absence of metric perturbations, Pj would remain constant, while 
Pi would decrease with the expansion like a~^. Usually, the formulas are conve- 
niently expressed in terms of the third variable qi = api, which can be viewed 
as the proper momentum "corrected" from the effect of homogeneous expan- 
sion. The proper energy of a single neutrino measured by a comoving observer 
is equal to (p^ + m^)^/^. We define the related quantity e 



From Pfj^P^ = and Pj = (1 — ■0)gj, we get Pq = (1 + 0)e. Finally, we recall 
that by definition the phase-space distribution / gives the number of particles 
in an infinitesimal phase-space volume: dN = f{x\ Pj, r) 0?%;* d^Pj- 



4-5.1 Collisionless fluids 

In the approximation of a homogeneous Universe, the phase-space distribution 
of neutrinos would be perfectly isotropic, 

f{x\Pj,T) = MP,T) , (85) 



where /o is the Fermi-Dirac distribution already introduced in Eq. (6). Then, 
the energy-momentum tensor would be diagonal with an isotropic pressure 
term 



47r r 

T^-Pu^— q'dqefoiq) , (86) 

Qj J 

Ti = -p. = -^Jq'dq^jfo{q). (87) 



Note that /o as a function of q is time- independent. 



because after neutrino decoupling the product of the neutrino temperature 
by the scale factor a is a constant number. 
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Spatial perturbations in the metric will induce variations in the neutrino 
phase-space distribution depending on time, space and momentum. The energy- 
momentum tensor will be perturbed accordingly. With respect to a perfect 
fluid, the main difference is is that, in absence of microscopic interactions, 
there is no reason for the stress tensor STij to be isotropic at first order in 
perturbations. Therefore, the scalar sector of ST^j^i, contains an extra degree of 
freedom, the anisotropic stress a. 

This can be understood in physical terms. First, strongly-interacting fiuids 
have bulk motions: the coarse-grained fiuid has a unique velocity in a given 
point. In the local rest-frame, an observer can measure the pressure e.g. by 
enclosing some particles in a small rubber balloon. At first order in pertur- 
bation, the balloon will expand or contract isotropically, because anisotropies 
arise from velocity gradients which can be shown to be second order in per- 
turbations. In contrast, in collisionless fiuids, there are no bulk motions: at a 
given point, there can be fiows of particles in all directions. Inhomogeneities 
and anisotropies are described perturbatively by a function ^ <C 1 such that 

f{x\Pj,T)^fo{q)[l + ^x\qj,T)] with Pj^{l-^)qj , (89) 



and anisotropic pressure can now appear at first order in perturbations. In- 
deed, the energy-momentum tensor can be computed from the phase-space 
distribution / and the 4-momentum 

T^, = J dP^dP^dPs i-g)-'/' ^f{x\ Pj, r) , (90) 
where {—g)~^^'^ — a~'^(l — (f) + 3t/j). At first order in perturbations, this gives 

T^{x') = a-^ J q^dq dQ e /o(g) [1 + ^{x\ qn„ r)] , 

T^{x')^a-^J q^dqdnqnifo{q)^{x\qnj,T), (91) 

2 

Tiix') = -a-^ I q^dqdn ^ fnhj fo{q) [1 + *(a;\ qhj, r)] , 

where dD, is the differential of the momentum direction rij — qj/q. It is then 
straightforward to obtain the perturbed components of the energy-momentum 
tensor 

Sp^ = a-^ j q^dq dQ e/o* , SP^ = ^a'^ j q^dq dQ ^/o* , (92) 
STj"^ = a-^ j q^dq dU qfn /o* , E} , = -a"^ j q^dq dU ^ {fnhj - -Sij) /o* . 
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4-5.2 Free- streaming 

We have seen that in perfect fluids, sound waves can propagate at the sound 
speed on scales smaller than the sound horizon. Sound waves cannot prop- 
agate in a coUisionless fluid, but the individual particles free-stream with a 
characteristic velocity -for neutrinos, in average, the thermal velocity Vth- So, 
it is possible to deflne an horizon as the typical distance on which particles 
travel between time tj and t. During MD and RD and for t ^ ti, this horizon 
is, as usual, asymptotically equal to Vth/H, up to a numerical factor of order 
one (see section 4.4.2). Exactly as we deflncd the Jeans length, we can deflne 
the free- streaming length by taking Eq. (48) and replacing by Vth 



As long as neutrinos are relativistic, they travel at the speed of light and their 
free-streaming length is simply equal to the Hubble radius. When they become 
non-relativistic, their thermal velocity decays like 

, M ^ = !Z2 (^^) ^ 160(1 + .) (i^)kms-' , (94) 

m m m \ a J \ m J 



where we used for the present neutrino temperature T° ~ (4/ll)^/^T° and 
T° ~ 2.726 K. This gives for the free-streaming wavelength and wavenumber 
during matter or A domination 



\Fs{t) = 7.7 , ^ + ^ f ^) h-'M^c , (95) 

k,s{i) = 0-82 ^ (^) h Mpc-\ (96) 

where (resp. fi^) is the cosmological constant (resp. matter) density frac- 
tion evaluated today. So, after the non-relativistic transition and during mat- 
ter domination, the free-streaming length continues to increase, but only like 
(aH)"^ oc i^/^, i.e. more slowly than the scale factor a oc t^^^. Therefore, the 
comoving free-streaming length Xps/a actually decreases like {o?H)~^ oc t"^/^. 
As a consequence, for neutrinos becoming non-relativistic during matter dom- 
ination, the comoving free-streaming wavenumber passes through a minimum 
/cnr at the time of the transition, i.e. when m — (p) — ST^, and ao/a— (1-I--2) = 
2.0 X 10^(m/leV). This minimum value is found to be 

1 /2 

k^, ~ 0.018 (^^) h Mpc-i . (97) 
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The physical effect of free-streaming is to damp small-scale neutrino density 
fluctuations: neutrinos cannot be confined into (or kept outside of) regions 
smaller than the free-streaming length, for obvious kinematic reasons. We 
will see that the metric perturbations are also damped on those scales by 
gravitational back-reaction. On scales much larger than the free-streaming 
scale, the neutrino velocity can be effectively considered as vanishing, and after 
the non-relativistic transition the neutrino perturbations behave like cold dark 
matter perturbations. In particular, modes with k < /Cnr are never affected by 
free-streaming and evolve hke in a pure ACDM model. 



4-5.3 Vlasov equation 

The neutrino equation of evolution is the coUisionless Boltzmann equation, 
or Vlasov equation, which just expresses the conservation of the number of 
particles along the phase-space trajectories 

dr dr dr dx^ dr dq dr drti 



We can expand / like in Eq. (89) and keep only first-order perturbations. Using 
the fact that dx^/dr — P'^/P^ — —P^/Pq — —qije at order zero, that dqjdr — 
qip + eni di(f) at order one (from the geodesic equation P^P^ + T^^P'^P'^ = 0), 
and that drii/dr = at order zero (in absence of metric perturbations, the 
momenta keep a fixed orientation), we obtain 

fo^-^fodi^ + iqi^ + erii di<P) ^ = , (99) 



or in Fourier space and after dividing by /o 

if-i^(k-h)-^^-(ij + i-(k-n)A^;^. (100) 
e V q ) omq 



This equation gives the response of the neutrino phase-space distribution to 
metric perturbations. 



4.. 5. 4 Neutrino perturbations during the relativistic regime 

When neutrinos are ultra-relativistic {e — q), the Vlasov equation reduces to 

4f-i(k-n)^^-('iP + i(k-n)(j))—^. (101) 
^ ^ omq 
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The homogeneous equation does not involve g, while the source term has a 
definite g-dependence, fixed once and for all by the shape of the homogeneous 
distribution. This simply reflects the fact that neutrinos at a given point and 
in a given direction are all redshifted or blueshifted by the same amount, so 
that their spectrum always remain that of a blackbody, with a temperature 
shift 5Ti,{x'', fij, t) at a given point and in a given direction 



f{x\q,hj,T)= exp 



+ 1 



dfo 



fo{q) + ^(g) dUx\n„T) . (102) 
Since /o is a function of q/aT^, one has dfo/dTi, — —{q/Ti) dfo/dq, and 



f{x\q,hj,T) = fo{q) 



dfo ST^ i ^ 



dlnq 



(103) 



Identifying this result with the definition of the phase-space perturbation ^, 
we get 



^{x\q,hj,T) 



-^{q) — (a;,n,,r) 



(104) 



We now see explicitly that the presence of the function ^^^^ in the Vlasov 
equation means precisely that the spectrum remains Planckian. Since for rel- 
ativistic neutrinos the local g-dependence is trivial, we only need to compute 
the evolution of the mean momentum, which is a function of position, direction 
and time. One can define the quantity F^{x\ fij, r) as 



Fi,{x\nj,T) 



Jq^dqqfojq) [1 + ^{x\hj,T)] 
!(fdqqfo{q) 



1 , 



(105) 



for which the Vlasov equation reads in Fourier space (after integrating by 
part) 



- i{k ■ h) = 4:(ip + i{k ■ h 



(106) 



We can further reduce the dimensionality of the problem by noticing that the 
evolution depends only on the direction h through the angle k ■ h, which is 
a natural consequence of the isotropy of the homogeneous background. We 
perform a Legendre expansion with respect to this angle, 

oo 

F,{k, h, r) = E(-0'(2^ + l)^w(^, T)PiCk ■ h) . (107) 

1=0 
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By comparing Eqs. (92) and (93) (with e = g) with the above definitions, it 
is straightforward to show that the density contrast (as well as Sp,, = ^Sp,^) 
is given by the monopole: 61, = F^q, the velocity gradient by the dipole: 61, = 
^kF^i, and the shear (or anisotropic stress) by the quadrupole: = — |-F1,2- 
The Vlasov equation now reduces to its most "concentrated" form: an infinite 

— * 

hierarchy of coupled equation for the multipoles Fi,i{k,T) 



4 

5. = -0. + 4z/; , 



k 



(/ + l)F,(,+i) , />2. (108) 



Note that the first two equations consistently reproduce the continuity and 
Euler equations (39)-(40). 

After introducing all this formalism, it is possible to understand the physical 
evolution of relativistic neutrino perturbations. We assume isentropic (adi- 
abatic) initial conditions during radiation domination and on super-Hubble 
scales: so, we start from = S-y and from time-independent metric fluctua- 
tions. On super-Hubble scales, the neutrino distribution is static, with a local 
density (and temperature) perturbation proportional to the metric fluctuation. 
An observer would see locally a quasi-isotropic temperature distribution: the 
dipole is suppressed by a factor k/aH — kr, and the l-th multipole by (/cr)'. 
At leading order in (kr), the Einstein equations (35) and (36) immediately 
give 

s, = s^ = -20 = 4^ , e, = e^ = -\kM , (109) 



while the expressions of ip and require more work (actually, the usual way 
to obtain them is to solve the Einstein equations in the long wavelength limit 
in the synchronous gauge, and then to convert the result into longitudinal 
gauge variables). The result is found to be (see Ref. [74]) 

= fl + (p , (7, = -^{kTf(p , with = _ . (110) 



For three standard non-relativistic neutrinos, on has R„ ~ 0.405 and ijj ~ 
1.162 0. After Hubble crossing, i.e. on scales smaller than the free-streaming 
length, the various multipoles are populated one after each other, at the ex- 
pense of the low multipoles which get reduced. Physically, this corresponds to 
the fact that an observer would see locally the superposition of a growing num- 
ber of neutrino flows, coming from all directions. Since the flows come from 
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Hubble crossing 




Hubble crossing 




Fig. 8. Evolution of (longitudinal gauge) massless neutrino density perturbations 
(5jy = bpyjpy (top), couipared with those of the photon density 5^ = 5p^/p^ (bot- 
tom), in a ACDM model with three massless neutrinos, as a function of the scale 
factor a and Fourier wavenumber A;, obtained numerically with a Boltzmann code. 
Time is evolving from left {a/ao = 10^^) to right {a/ao = 1), and the initial condi- 
tion is set arbitrarily to k^^'^cf) = —10^^. Unlike the photon density contrast, which 
oscillates between Hubble crossing and equality, the neutrino density contrast is 
immediately damped after the first half-oscillation. After equality, 5^ and re- 
main constant: they are preserved from gravitational collapse by free-streaming and 
relativistic pressure. 

a random distribution of over- and under-densities, the local density contrast 
tends to average to zero. 



45 



1.0 10"^ 



0.8 10"^ 



0.6 10"^ 



0.4 10"^ 
0.2 10"^ 
0.0 10"^ 
-0.2 10"^ 

10"^ 10"^ 10"^ 10"^ lO""" 1 

a/ao 

Fig. 9. Evolution of some of the first neutrino multipoles F^i as a function of the 
scale factor, for the mode k = 0.1/iMpc^^ which enters into the Hubble radius 
close to the end of radiation domination, when a/ag = 2.5 x 10~^. We imposed an 
initial condition k^/'^ij) = — 10~^. We see that after Hubble crossing, the multipoles 
are populated one after each other, reflecting the fact that an observer would see 
locally a superposition of a growing number of neutrino flows, coming from all 
directions. 

It is useful for pedagogical purpose to compare the evolution equation of neu- 
trinos with that of photons before recombination. For this purpose, we combine 
the continuity and Euler equations into a second-order equation 

K = -\kH, - + 4V^ - h^a, . (Ill) 

In absence of shear, this equation would be similar to the photon equation of 
evolution before recombination and in the limit of small baryon-to-photon ra- 
tio Rh = pb/p-y- The terms are easy to identify: —\k'^5y is the force arising from 
relativistic pressure, —^k'^cj) the gravitational force, and Axjj the dilation effect: 
a local increase of ijj is equivalent to a local decrease of the scale factor and of 
the photon/neutrino wavelength, and hence to an increase of the blackbody 
temperature and density. In the case of photons, in the tight-coupling limit, 
the microscopic interactions guarantee that vanishes (even if <C 1), and 
the equivalent of the above equation is that of an harmonic oscillator driven 
by metric fluctuations. In the case of neutrinos, the shear becomes important 
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inside the Hubble length, and couples the above equation to the whole multi- 
pole hierarchy; qualitatively, the shear acts like viscosity in a fluid; instead of 
experiencing driven oscillations, the neutrino overdensity 5y is damped, while 
the power is transfered to higher multipoles. This behavior is illustrated in 
Figs. 8 and 9 which derive from exact numerical simulations (using the code 
CMBFAST [49]). In particular, one can sec explicitly on Fig. 8 that in the re- 
gion where the photon density contrast experiences acoustic oscillations, the 
neutrino density contrasts is damped. 

It is possible to gain even more intuition by going back to the relativistic 
Vlasov equation (101), for which the full solution can be decomposed into those 
of the homogeneous and inhomogeneous differential equations. The general 
solution of the homogeneous equation is actually trivial 

n, r) = - ^(^, n) ^^'^^ ^ . (112) 



It corresponds to the free propagation of plane waves of temperature pertur- 
bations, i.e. to free-streaming. In the Fi,i basis, this yields F^i{k,T) oc ji{kT), 
where the ji^s are the spherical Bessel functions of the first kind, and the con- 
tribution to 5i, is equal to [sin(A;r)/(A;r)] times a constant of integration. As 
expected, this contribution tends to zero after Hubble crossing (because of 
the averaging between various flows), while higher multipoles are populated 
one after each other, with a peak around r ~ ///c (the numerical solution of 
Fig. 9 looks slightly more comphcated than plain Bessel functions: this is due 
to the small time- variation of and ip caused by neutrino shear). On top of 
these terms, the full solution receives a contribution from the solution of the 
inhomogeneous equation. During radiation domination, tjj and (p decay, and 
this solution is complicated. However, it is clear that it should also decay. 
Instead, during matter domination, ip and (p are frozen, and the solution is 
obvious: = In the F^i basis, this term only contributes to F^^q — 5^, 

for which the full solution reads 

5^ = -A(f)(k) + a{k) sm(kT)/{kT) , r > Teq , (113) 



where the constant of integration a{k) could be found my matching with 
the radiation-dominated solution. The flrst term corresponds to the station- 
ary equilibrium conflguration for which, in any potential well, the relativistic 
pressure compensates exactly the gravitational force. This behavior can be 
clearly seen on Fig. 8 (note that after decoupling, the photons free-stream like 
the relativistic neutrinos, so 5-y also tends towards —4 0). 

In conclusion, during the relativistic regime, the combined effects of free- 
streaming and relativistic pressure result in a suppression of the neutrino 
density contrast inside the Hubble radius: during radiation domination, 5i, is 
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driven to zero instead of oscillating like 5^ and ^b, and during matter domi- 
nation, it remains constant instead of growing like (5b and (5cdm- 



4.5.5 Neutrino perturbations during the non-relativistic regime 

After the non-relativistic transition, the g-dependence of the phase-space dis- 
tribution becomes non-trivial. The geodesic equation 

^^q^^^q^^a'm^f/'h,d,4> (114) 
dr 



shows that the amount of blueshift or redshift dq/q experienced by a neutrino 
when it travels over metric fluctuations depends on how non-relativistic it is. 
In a given point and direction (a;*,nj), there is a superposition of neutrinos 
of different momenta which are more or less deep inside the non-relativistic 
regime. Lowest and highest energy neutrinos experience different redshifts, and 
the thermal spectrum of Eqs. (102) and (103) is inevitably distorted at first 
order in perturbation (but not at zero order, since the distortion is induced 
by metric perturbations and not by the background geometry). Therefore, it 
is impossible to eliminate the g-dependence from the problem, as we did in 
the rclativistic regime. However, it is still possible to get rid partially of the 
direction dependence by expanding ^ in multipole space 

00 

*(^, g, n, r) = E(-0'(2/ + 1)*K^, r)Pi{k ■ h) . (115) 



Then, the Vlasov equation becomes 



= — *i - V'^i , 

e am g 

qk , , ek dXn fn , ^ 

^1 = ^ ^2 - 2^0 + ^ -KT^ , 116 
3e 3q omq 

= ^^j^[(^ + 1)^^+1 - , ^>2. 

The perturbed energy density, pressure, energy flux divergence and shear stress 
in k-space are given by 



Sp„^4:TTa ^ j q^dqefo{q)'^o, 

4:71 f g^ 
Spu = Y"'~^ J ^^^^ — /o(?)*o , 
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Fig. 10. Evolution of the neutrino isotropic pressure perturbation 6pi, and 
anisotropic stress ai, in units of the density perturbation (for the mode 
k = 0.1 /iMpc~^ which enters into the Hubble radius when a/uQ = 2.5 x 10~^) com- 
pared to the evolution of the background equation of state parameter Pu/pu- The 
thin curves show the case of massless neutrinos for which Pu/pu = ^Pu/^Pu = 1/3 
and (Jylby decays slowly after Hubble crossing. The thick lines correspond to neutri- 
nos with rriy = 0.1 eV, which become non-relativistic around a/ao = 5 x 10"'^. After 
the transition, pressure and shear perturbations become negligible with respect to 
density perturbations, just like for cold dark matter. 

{p„ + pu)du = 4:rrka~^ J g^(igg/o(g)^i , 

(p. + P.)a. = -y J q'dq ^^^2 • (117) 

When for a given family the neutrinos are well inside the non-relativistic 
regime, i.e. when most of the momenta for which fo{q) is non-negligible satisfy 
q e am, the above equations show that both 6pi, and {pu + py)ay are 
suppressed with respect to 5py (as can be checked in Fig. (10)). As expected, 
the continuity and Euler equations (39) and (40) become gradually identical 
to those of ordinary non-relativistic matter 



5. = ^. + 3^ 



a 



(118) 
(119) 
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Like for CDM (and baryons after recombination), this system leads to the 
second-order equation of evolution 



6, + -6, = -k'(i) + 3{ij + -ilj). (120) 
a a 



We are interested in neutrinos becoming non-rclativistic during matter dom- 
ination, when (p and ip are approximately constant and a oa r^. In this case, 
the solution of (120) reads 



(5, = AlnT + S-^^0 (121) 

-^'-^^^-\(^)^^ ^^22) 

where A, B (or A, B) are some constants of integration. The last term (the 
solution of the inhomogeneous equation) grows like r^, i.e. like the scale factor. 
This solution has a simple interpretation: 

• first, for k > /Cnr, i-e. for modes which are inside the Hubble radius at 
the time of the non-relativistic transition, the density contrast 5^, is much 
smaller than 5cdm at that time (as a consequence of free streaming), and 

then grows faster than 5cdm oc a because of the term (A In a). At some point 
the first and third terms in Eq. (122) will be equal, and then asymptotically 
Si, will be dominated by 




which corresponds exactly to the solution of the Poisson equation in a 
matter-dominated Universe (see Eq. (76) with Qra=l). Indeed, this asymp- 
totic value is simply the one imposed to all non-relativistic species by New- 
tonian gravitation, and after an infinite amount of time 6^ would be equal to 
6cdm on all scales. The same is also true for baryons after photon decoupling, 
but since baryons did not experience free-streaming, the equality between 
(5b and Scdm is reached quickly after r = Tdec- Instead, for neutrinos with 
masses of order of 1 eV or smaller, the asymptote is far from being reached 
today, except for scales very close to A; = km- On smaller scales, Si^ remains 
much smaller than ^cdm — ^b- 



• second, for k < /Cnr, i-e. for modes which are outside the Hubble radius at 

the time of the non-relativistic transition, the density contrast is time- 
independent and of the same order as (0, 5cdm) on super-Hubble scales. After 
horizon crossing, the solution is dominated again by (123) and becomes 
rapidly equal to 5cdm- 
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The conclusion of this section is that if we could picture the three density 
contrasts at a fixed time near the end of matter domination, we would expect 
= ^cdm = for GoHq < k < k^j-, while for k > fcnr the ratio 5y/5cdm = ^i//^b 
would be smaller than one, and tending towards zero as k goes to infinity. 
The recent stage of A domination does not change this conclusion: for qoHq < 
/c < /cnr, the three density contrasts are damped by the same amount, and on 
smaller scales the ratio 6iy/6cdm does not change significantly. All these features 
can be checked in the numerical simulations presented in Sec. 4.5.7 (see Fig. 
11). Therefore, the matter power spectrum today (or at any time after the 
non-relativistic transition of the heaviest neutrinos) is given by 



P{k)^ 



' 5 Pcdm + 5Ph + 5 P. 
, Pcdm + Ph + Pu 

^cdm '^cdm ~l~ ^b '^b 




(^cdm) 

[1 - {si^J 




for k < k-c, 
for k :$> kr. 



(124) 



with Qm = flcdm + fih + ^u- If the neutrinos did not induce a gravitational 
backreaction effect on the evolution of the metric and of other perturbations, 
then (^cdm) would be the same with and without neutrinos and the role of 
the neutrino masses would be simply to cut the power spectrum by a factor 
[1 — fliy/flm]"^ for k S> knr- However, we will see in the next subsection that the 
presence of neutrinos actually modifies the evolution of the CDM and baryon 
density contrasts in such way that the suppression factor is greatly enhanced, 
more or less by a factor four. 



4-5.6 Modified evolution of matter perturbations 
Neutrino backreaction effects: generalities. 

In principle, the impact of neutrinos on the evolution of other perturbations 
(baryons, cold dark matter, photons) can be divided in two categories: 

(1) effects from the homogeneous density and pressure of neutrinos, which 
affect the Friedmann law. Since the equations of evolution of all pertur- 
bations contain the expansion rate, the solutions can be affected. 

(2) direct gravitational back-reaction effect at the level of perturbations. The 

neutrino density, pressure and velocity perturbations provide a signifi- 
cant contribution to the total STj^, essentially on scales larger than the 
free-streaming scale. So, they can change the evolution of metric pertur- 
bations, in a similar way as an increased photon density. In addition. 
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the neutrino shear induces a non-zero difference {(j) — ip), which is con- 
stant on super-horizon scale during radiation domination and becomes 
progressively negligible on all scales during matter domination. 

Some of these effects are independent of the neutrino mass, and their combi- 
nation results in a very significant difference between the CMB and LSS power 
spectra in absence of neutrinos, and those in presence of three families of ultra- 
relativistic neutrinos (see Fig. 14). These effects are investigated in details in 
[91], and are extremely relevant from the observational point of view. For in- 
stance, they give an opportunity to measure the effective neutrino number A^eff 
(see Sec. 3.2). Also, it was shown in Ref. [92] that the specific gravitational 
backreaction effect induced by neutrino shear is actually confirmed by current 
observations of the CMB temperature spectrum by WMAP. 

However, in this review, we are more interested in the effect of the neutrino 
mass, in realistic scenarios for which the non-relativist ic transition of at least 
two species takes place during matter domination. Our goal is to understand 
how the matter power spectrum of Eq. (124) depends on the neutrino mass 
in the two limits k < /Cnr and k > /Cnr- Therefore, we should study how the 
evolution of 5cdm = is affected by non-relativistic neutrinos deep in the 
matter-dominated regime. 

How neutrinos slow down the growth of matter perturbations. 

We are arriving to the most important point of this section, since we are going 
to derive the dominant observable effect of neutrino masses. Before getting into 
the equations, let us try to describe the effect qualitatively. For this purpose, 
we first need to come back to the issue of linear clustering of CDM and baryons 
during matter domination in absence of neutrinos. In Sec. 4.4, we stated that 
(f) — ip remains constant during matter domination, basing ourselves on Eq. 
(58) with a zero sound speed (or equivalently on the perturbed Einstein Eqs. 
(37) and (38) with zero pressure perturbation and shear on the right-hand 
side). Then, we used the Poisson equation 

- — ^ = iirGSp , (125) 

to show that inside the Hubble radius, the density contrast Scdm = grows 
like [a^(pcdm + Ph)]~^, i-e. like a linear function of the scale factor. The fact 
that linear structure formation corresponds to a constant gravitational poten- 
tial sounds counter- intuitive at first sight. Actually, it results from an exact 
compensation between clustering and the stretching of spacetime. During grav- 
itational collapse, the non-relativistic matter density perturbations Sp do not 
grow with time, as one would naively expect from a common intuition based 
on a static Universe. Instead, 5p decreases, but only like and not as fast 
as the dilution factor a"^, so that the density contrasts Sp/p actually grows 
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like a. The behavior 5p oc ar^ corresponds precisely to a static potential, since 
the gradient A-^ — —{k/a)'^ip also decreases like due to the stretching of 
space-time. 

Let us restore now the neutrino background. During matter domination and on 
scales smaller than the free-streaming scale, the neutrino perturbations 6iy do 
not contribute to gravitational clustering. Indeed, we have seen in Sees. 4.5.4 
and 4.5.5 that free-streaming leads to 61, <^ ^cdm, and since pi, < Pcdm we see 
that 5pi, <^ ^PcAm- so, neutrinos can be simply omitted from the Poisson equa- 
tion. On the other hand, they do contribute to the homogeneous expansion 
through Friedmann equation. Therefore the exact compensation between clus- 
tering and expansion described in the previous paragraph is slightly shifted: 
the balance is displaced in favour of the expansion effect, and we expect (p = ip 
to decay slowly, while 5cdm = should grow not as fast as the scale factor. 
This mechanism, first described and quantified in Ref . [70] , belongs to the first 
of the two categories defined at the beginning of this subsection: physically, it 
is a pure background effect, which leads to a modified evolution of metric and 
matter perturbations. 

Let us estimate the magnitude of this effect in the same fashion as in the 
pioneering paper [70]. At any time we can combine the continuity and Euler 
equations of CDM perturbations (39) and (40) into 

^cdm + - 4dm = -/cV + 3(V; + -^) , (126) 

a a 

but well-inside the Hubble radius the source term is dominated by the co- 
moving gradient —k'^(p, which is approximately equal to —k'^ip as soon as the 
neutrino shear can be neglected. This gradient is given by the Poisson equa- 
tion, and we obtain 

4dm + - 4dm = 4710? a" Sp , (127) 

where 5p is the total density perturbation. In absence of neutrinos and deep 
inside the matter-dominated regime, we would have Scdm — Sh, a. total density 
perturbation Sp = (pcdm + Pb) ^cdm, an expansion rate given by 3(d/a)^ = 
87rGa^(pcdm + Ph) oc a~^, a scale factor a oc r^, and the equation would read 

2 ■ 6 

<^cdm H Scdm o ^cdm — , (128) 

with two solutions Scdm oc and Scdm oc r~^. Neglecting the decaying mode 
we recover the standard result 5cdm oc a. Now, let us consider instead the case 
with massive neutrinos, still deep inside the matter-dominated regime and 
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on scales k ^ knr, so that Spi^ does not contribntc to the Poisson cqnation: 
5p = (pcdm + Pb) ^cdiii! while the neutrino background density does contribute 
to the expansion rate: 3 {a/af = 87rGa^(pcdm + Pb + Pi^)- Let us assume that 
Pi, is dominated by non-relativistic neutrinos, so that it decays approximately 
like a~^, and the number 

/. ^ 7 T—T-^ = ^ (129) 

{Pcdm + Ph + Pu) i'm 



remains approximately constant. Then, the scale factor still evolves like 
and the equation of evolution reads 

2 • 6 

Scdm + - Scdm ^(1 - U) <^cdm = . (130) 

T T 

Looking for solutions in 5cdm oc r^^, we find two roots 



-1± ^1 + 24(1-/.) 
P± = A > (131) 



and we conclude that the growing solution for the CDM density contrast reads 

(5edm « aP+ ~ a'-lf-^ , (132) 



where in the last step we assumed /j, ^ 1. As expected, the growth of 6cdm is 
reduced due to the fact that one of the component in the Universe contributes 
to the homogeneous expansion rate but not to the gravitational clustering. 
The Poisson equation gives 

- fcV oc aP+-^ ~ a-t^-- , (133) 



showing that for the same reason the gravitational potential slowly decays 
during matter domination. 

At the end of matter domination and during A domination, we have already 
seen that in absence of neutrinos (p decays like g{a) and ^cdm grows like ag{a) 
(we recall that the damping factor g{a) is normalized to = 1 for a <C ca). 
The combined effect of A and of neutrinos on the growth of Scdm can be well 
approximated by [93] 

5cdm OC [ag{a)f+ ~ \a gia^-l^'^ . (134) 



Matter power spectrum for massive versus massless neutrinos. 
Let us try to predict analytically the difference between the power spectrum 
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of two cosmological models with respectively massive and massless neutrinos 
(AMDM versus ACDM) assuming that they share the same values of (c^m, cjb, 
fl\, A, Us, t). The difference between the two models lies in the values of cUi, 
and tUcdm = f-i^m — (^b ~ '^i/, and can be conveniently parametrized by f^. 

(1) for very large scales k < knr, neutrino perturbations were never affected 
by free streaming, and in the non-relativistic regime they become indis- 
tinguishable from CDM perturbations. In particular, inside the Hubble 
radius, neutrinos contribute both to the Poisson equation with 61, = 5cdm 
and to the background evolution, so the density contrast grows like a. So, 
for fixed cUm, this branch of the spectrum is completely insensitive to the 
value of fi, and of the neutrino mass. 

(2) for very small scales with k ^ k^v and k ^ /Cgq, the two spectra will differ 
for various reasons. First, equality does not take place for the same value 
of the scale factor in the two models, since one always has (aeq/ao) = 
C(Ji/(^b + ^cdm) = (1 ~ fu)~^i^i/i^m, whcrc uj^ iucludcs the density of 
photons and of three effectively massless neutrinos. Since the two models 
share the same value of Uj. and cjm, the ratio between the two values of 
the scale factor at equahty is Oeq/'^^q^^ = ~ fi)~^- At equahty, and 
more generally at any time before the non-relativistic transition, the two 
models are rigorously equivalent, modulo a shift in the scale factor by the 
above factor. So, the identity 

^fdmM=^fd:°[(l -/.)«] (135) 

is valid as long as a < Onr- After a = Onr and in the massive model, 
the mechanism described in Eqs. (130)-(133) will take place. So, in the 
approximation of a time-independent fraction /j, for a > a^r (which is not 
strictly valid in the case of several massive non-degenerate neutrinos), the 
growth of CDM perturbations is given by Eq. (134), i.e. by a constant 
logarithmic slope (1 — ^f^) for the scale factor (or for the product [ag{a)] 
during dark energy domination), 

«c'3.M=f^)"''"<Sc'3„Kl. (136) 

In the massless case, everything is identical to the massive model until 
the time at which a = (1 — /i^) Onr (i-e. until the value of the scale factor 
for which neutrinos become non-relativistic in the massive case "shifted 
back" to the massless case according to Eq.(135)). Therefore, at this time, 
the logarithmic slope of 5fdm° must also be equal to (1 — |/i/), but then 
it will increase progressively until unity, so that Sl^~^ becomes a linear 
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function of a or [a g'( a)]. In a crude approximation, we can write 



but this tends to overestimate the growth of perturbations in the massless 
case: it assumes that right after a = the logarithmic slope is equal to 
one, which is not true immediately. Indeed, a comparison with numerical 
results shows that the total growth factor is a bit smaller. 

Using this semi-analytic result, we find that the ratio between the present 
value of Scdm in the two models reads 

|2^ = (l-«"=(^)"''. (139) 

<^cdm [«0] V «nr J 

According to Eq. (124) this means that the total matter power spectrum 
is reduced by 

(HO) 



P{k) 

Finally, we can replace (ao/onr) by 2000mj,/(l eV) and, assuming that 
the mass rriu is shared by a number iVj, of families, we can use ttIi, ~ 
[ujy/Ny) 93.2 eV. We obtain an expression that depends only on (/j,, N^,, 



{l-Uf 1.9xlO'g{ao)uj^U/K ■ (Ml) 



We show in Fig. 12 that this semi-analytic expression is a very good 
approximation of the exact numerical result, and also that for plausible 
values of (a;in, N^,, ft\) and for f^, < 0.07, it can be approximated by the 
well-known linear expression [94] 

^^^^^"^ - -8/.. (142) 



In conclusion -and as confirmed by the numerical simulations in the next 
subsection- the ratio P{k)^" / Plk)^"^^ computed for fixed parameters {ui^, 
Qa) smoothly interpolates between one for k < k-^r and a plateau given by Eq. 
(141) for k > ^eq- 
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Fig. 11. Evolution of the metric and density perturbations as a function of the 
scale factor (normalized to ao = 1 today), in the longitudinal gauge, for modes 
10~^/iMpc~^ < k < 1/iMpc^^ (from top to bottom), and for two cosmological 
models: ACDM (left) and AMDM (right), both with uj^ = 0.147 and S^a = 0.7. The 
integration has been performed with the code cmbfast starting from the initial 
condition k^^'^cp = —10^^. The AMDM model has three degenerate neutrinos with 
nil, = 0.46 eV, corresponding to = 0.1. 
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^.5.7 Numerical results 

In Fig. 11, we tried to summarize all the features described in this section by 
showing the numerical solutions for the metric and density perturbations {(f), 
"0) 4dm, (^b, (^7, Sj,) at different wavenumbers, for a ACDM (left) and AMDM 
(right) a model. The two models share common values of (o^m, i^b, ^A, t) 
and the initial amplitudes are the same. In the ACDM model, the three neu- 
trino families are massless, while in the AMDM model they share a common 
value of the mass corresponding to fi, — 0.1. At the bottom of each plot, we 
show the value of the scale factor at the time of Hubble entry (k — aH), 
radiation/matter equality (oeq) and matter/A equality (ca). Decoupling takes 
place soon after equality, around a = 9 x 10~"^. These solutions have been com- 
puted with the public code cmbfast [49] in the synchronous gauge, and then 
translated into the longitudinal gauge. Let us briefly summarize the behavior 
of each quantity. 

For a ACDM model (left plots in Fig. 11), one has that 

• the metric perturbations are essentially constant on large scales k < keq dur- 
ing radiation and matter domination. Instead, on small scales k > keq, they 
decay between Hubble entry and decoupling, damped by the propagation of 
acoustic waves. On all scales, they are also damped during A domination. 
Initially the difference — is constant and non-zero, due to neutrino shear. 
After equality, the difference becomes negligible. 

• the photon perturbations undergo acoustic oscillations between Hubble en- 
try and decoupling. Then, during matter domination, they free-stream, 
maintaining and average density contrast = —40, corresponding (as we 
have seen) to gravity vs. pressure equilibrium. 

• the relativistic neutrino perturbations are damped during radiation domina- 
tion after Hubble crossing, because of free-streaming. During matter dom- 
ination free-streaming goes on but Si^ maintains itself around —40 like for 
photons. 

• the CDM perturbations always grow inside the Hubble radius. During radi- 
ation domination this growth starts abruptly, but then becomes only loga- 
rithmic; during matter domination it is linear; finally during A domination it 
slows down. In the bottom left plot for /c = 1 hMpc~^, the solid red line is a 
linear function of a normalized in order to flt 5cdm during matter domination. 

• the baryon perturbations follow closely photon perturbations until decou- 
pling, experiencing the same acoustic oscillations. After decoupling, they 
grow very quickly towards (5cdm, and after a ~ 0.01 the baryon perturba- 
tions remain equal to the CDM ones. 
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Fig. 12. Ratio of the matter power spectrum including Ni, = 1 or 3 degenerate 
massive neutrinos to that with three massless neutrinos, computed at the scale 
/c* = 5/iMpc~^ for fixed parameters (wmj ^a) = (0.147,0.70), and plotted as a 
function of the density fraction f^. The numerical result is compared with the 
semi-analytical approximation of Eq. (141) and with the linear approximation —8/^. 



For a AMDM model (right plots in Fig. 11), the main differences with respect 
to the previous case are 

(1) the fact that equality takes place a bit later, by a fraction (1 — fu)~^ = l-l- 

(2) the behavior of S^, after the time of the non-relativistic transition, which 
is shown at the bottom of the plots by the label "flnr"- For a > Onr, the neu- 
trino density perturbations track 5cdm = instead of 6^. For large scales with 
k < k-oT (upper plot), the non-relativistic transition takes place when the wave- 
length is larger than the Hubble radius, so 5cdm and 5y are not very different 
when and soon after 5y reaches (5cdm- For slightly smaller scales like 
k = 10~^/;,Mpc~^ (right plot in second row), the neutrino perturbations are 
much smaller than the CDM ones when a = am-, but then they have time 
to grow and reach ^cdm before the present time. For k > 10~^/iMpc~^ (lower 
plots), they do not have time to do so and remain subdominant today. 

(3) the growth of CDM perturbations is suppressed in the AMDM case for two 
reasons: first because equality takes place a bit later, and more importantly 
because 6cdm does not grow as fast as the scale factor for a > a^^. This can be 
seen very clearly on the lower right plot, where we show exactly the same linear 
function as in the lower left one. The change in the time of equality shifts the 
(5cdm normalization down, while the effect of non-relativistic neutrinos lowers 
the value of the slope as predicted analytically in the last subsection. 
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Fig. 13. Ratio of the matter power spectrum including three degenerate massive 
neutrinos with density fraction fu to that with three massless neutrinos. The pa- 
rameters {ujra, ^a) = (0.147, 0.70) are kept fixed, and from top to bottom the curves 
correspond to fu = 0.01, 0.02, 0.03, . . . , 0.10. The individual masses niu range from 
0.046 eV to 0.46 eV, and the scale from 2.1 x lO^^/iMpc"^ to 6.7 x 10^^/iMpc"^ 
as shown on the top of the figure, keq is approximately equal to 1.5 x 10^^/iMpc~^. 

Looking now at all wavenumbers, we plot in Fig. 12 the ratio of the matter 
power spectrum for AMDM over that of ACDM, for different values of f^, 
but for fixed parameters (cUm, ^a)- Here again, the AMDM model has three 
degenerate massive neutrinos. As expected from the analytical results, this 
ratio is a step-like function, equal to one for k < k^r and to a constant for 
k ^ keq. The value of the small-scale suppression factor is plotted in Fig. 13 
as a function of f^, and of the number N^, of degenerate massive neutrinos, still 
for fixed {oo^, ^a)- The numerical result is found to be in excellent agreement 
with the analytical prediction of Eq. (141). For simplicity, the growth factor 
g{aQ) ~ 0.8 can even be replaced by one in Eq. (141) without changing the 
result significantly. The well-known formula P{k)^'' / P{k)^''^^ ^ —8fu is a 
reasonable first-order approximation for < < 0.07. 

4-6 Summary of the neutrino mass effects 

4-6.1 Effects on CMB and LSS power spectra for fixed {u^, VL^) and degen- 
erate masses 

In Fig. 14, we show Cf and P{k) for two models: ACDM with fu = and 
AMDM with = 3 massive neutrinos and a total density fraction = 0.1. 
We also display for comparison the neutrinoless model of Sec. 4.4.6. In all 
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Fig. 14. CMB temperature anisotropy spectrum Cf and matter power spectrum 
P{k) for three models: the neutrinoless ACDM model of section 4.4.6, a more re- 
alistic ACDM model with three massless neutrinos {f^ ~ 0), and finally a AMDM 
model with three massive degenerate neutrinos and a total density fraction fj^ = 0.1. 
In all models, the values of (wb, ^m, ^A; ^s, n, r) have been kept fixed. 
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Fig. 15. CMB temperature anisotropy spectrum Cf and matter power spectrum 
P{k) for three models: the same ACDM model as in the previous figure, with three 
massless neutrinos {fy ~ 0); and two models with three massive degenerate neutri- 
nos and a total density fraction f^, = 0.1, sharing the same value of and oJcdm as 
the massless model, which implies a shift either in h (green dashed) or in Q,\ (blue 
dotted). 



models, the values of (wb, tUm, f^A, ^s, r) have been kept fixed, with the 
increase in uj^ being compensated by a decrease in cUcdm- There is a clear 
difference between the neutrinoless and massless neutrino cases, caused by a 
large change in the time of equality and by the role of the neutrino energy- 
momentum fluctuations in the perturbed Einstein equation [91]. However our 
purpose is to focus on the impact of the mass, i.e. on the difference between 
the solid (red) and thick dashed (green) curves in Fig. 14. 

Impact on the CMB temperature spectrum. For fy < 0.1, the three 
neutrino species are still relativistic at the time of decoupling, and the di- 
rect effect of free-streaming neutrinos on the evolution of the baryon-photon 
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acoustic oscillations is the same in the ACDM and AMDM cases. Therefore, 
the effect of the mass is indirect, appearing only at the level of the background 
evolution: the fact that the neutrinos account today for a fraction fl^, of the 
critical density implies some change either in the present value of the spatial 
curvature, or in the relative density of other species. In this section, we choose 
to maintain a flat Universe with fixed (cjm, t^b, ^^a) -other options would be 
possible, as discussed in the next section. Thus, while and Oa are con- 
stant, flcdm is constrained to decrease as fl^, increases. The main effect on the 
CMB anisotropy spectrum results from a change in the time of equality. Since 
neutrinos are still relativistic at decoupling, they should be counted as radia- 
tion instead of matter around the time of equality, which is found by solving 
Ph + Pcdm ^ Pj + Pu- This gives 

aeq = , (143) 



where co,. stands for the radiation density extrapolated until today assuming 
that all neutrinos would remain massless, given by Eq. (11) with A^^cfr = 3.04. 
So, when f^y increases, Ocq increases proportionally to [1 — equality is 

postponed. We already saw in Sec. 4.4.7 that this produces an enhancement of 
small-scale perturbations, especially near the first acoustic peak. Also, post- 
poning the time of equality increases slightly the size of the sound horizon at 
recombination. These two features explain why in Fig. 14 the acoustic peaks 
are slightly enhanced and shifted to the left in the AMDM case. 

If neutrinos were heavier than a few eV, they would already be non-relativistic 
at decoupling. This case would have more complicated consequences for the 
CMB, as described in [85]. However, we will see in section 5 that this situation 
is disfavoured by current upper bounds on the neutrino mass. 

Impact on the matter power spectrum. This issue has been discussed 
in details in the previous subsection: the combined effect of the shift in the 
time of equality and of the reduced CDM fluctuation growth during matter 
domination produces an attenuation of small-scale perturbations for k > /Cnr, 
obeying to Eq. (141) in the large k limit. Various studies, including Refs. 
[93,95,96,97,98,99], have derived some analytical approximation to the full 
MDM or AMDM matter power spectrum, valid for arbitrary scales and red- 
shifts. However, the A 7^ case is only considered in Refs. [93,98]. 



4.6.2 Degeneracy with other cosmological parameters 

Possible parameter degeneracy in the CMB temperature spectrum. 

Since the effect of the neutrino mass on CMB fluctuations is indirect and ap- 
pears only at the background level, one could think that by changing the value 
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of other cosmological parameters it would be possible to cancel exactly this ef- 
fect. For the simplest AMDM model with only seven cosmological parameters, 
this is impossible: one cannot vary the neutrino mass while keeping fixed Oeq 
and all other quantities governing the CMB spectrum, that we already listed 
in Sec. 4.4.7. Let us illustrate this with a particular example. Since the time 
of equality is related to {uji, + o^cdm) rather than cUm, we could try to maintain 
this sum fixed. In this case, it is now uj„^ that will increase like [1 — fu]^^- Wc 
still have some freedom concerning other parameters: we could choose either 
to maintain a constant f^A, in which case h will increase as [1 — or 
to maintain a constant h, in which case (1 — Qa) will increase as [1 — fi]'^. 
The result of these two possibilities is displayed on Fig. 15. Since the value 
of the scale factor at equality is now more or less fixed, the change in peak 
amplitude is smaller than in Fig. 14. The counterpart is a larger shift in peak 
location: since we now have variations of either /i or JIa, the angular scale of 
the sound horizon at decoupling cannot be kept constant. 

We conclude that for a seven parameter AMDM model, there is no parameter 
degeneracy in the CMB spectrum. Therefore, it is possible to constrain the 
neutrino mass using CMB experiments alone (although neutrinos are still rel- 
ativistic at decoupling). We will see in the next sections how this applies to 
current and future experiments. 

This conclusion can be altered in more complicated models with extra cos- 
mological parameters. For instance, allowing for an open Universe, it would 
be possible to increase fi, with fixed (o^cdm, <^b; ^a); and to tune the spatial 
curvature in order to shift back the peaks to their original location. One could 
also vary the number of relativistic degrees of freedom in order to play with the 
time of equality (see e.g. [100]). In such extended models the CMB alone is not 
sufficient for constraining the mass, but fortunately the LSS power spectrum 
can lift the degeneracy. 

Possible parameter degeneracy in the matter power spectrum. The 

possibility of a parameter degeneracy in the matter power spectrum depends 
on the interval [/Cmin, ^max] in which it can be accurately measured. Ideally, 
if we could have k^in < 10~^/iMpc~^ and fcmax > l/iMpc~^, the effect of 
the neutrino mass would be non-degenerate, because of its very characteristic 
step-like effect. In contrast, parameters like the scalar tilt or the tilt running 
change the spectrum slope on all scales, while the cosmological parameters 
which determine the shape of the matter power spectrum below the scale of 
equality change the slope for all wavenumbers k > k^q- 

The problem is that usually, the matter power spectrum can only be accu- 
rately measured in the intermediate region where the mass effect is neither 
null nor maximal: in other words, many experiments only have access to the 
transition region in the step-like transfer function. In this region, the neutrino 
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mass affects tlie slope of tlie matter power spectrum in a way wfiicfi can be 
easily confused with the effect of the scalar tilt, tilt running, or many other 
cosmological parameters. Because of these parameter degeneracies, the LSS 
data alone cannot provide significant constraints on the neutrino mass, and 
it is necessary to combine them with CMB data. The latter usually lift all 
degeneracies. For instance, the CMB probes primordial fluctuations on a very 
broad range of scales and makes a very good determination of the tilt and 
tilt running: so, degeneracies with Ug or dug/dk are usually not a problem. We 
will see however in Sec. 5.6 that when exotic models with e.g. extra relativistic 
degrees of freedom, a constant equation-of-state parameter of the dark energy 
different from —1 or a non-power- law primordial spectrum are introduced, the 
neutrino mass bound can become significantly weaker. 



4-6.3 Effects of non- degenerate neutrino masses 

It is often said that cosmological perturbations only probe the total neutrino 
mass J2i rui ~ mi + m2 + (from now on denoted as M^). Let us scrutinize 
this point. 

We saw in the previous subsection that small neutrino masses only have an 
indirect background effect on the CMB anisotropy spectrum, related to the 
current density cjj,, and therefore to the total mass (see Eq. 12). So, in 
very good approximation, CMB anisotropics probe the total mass with no 
information on the splitting between different mass eigenstates. 

As far as LSS is concerned, the situation is a bit more subtle. In Sec. 4.5, 
we always assumed for simplicity that there was one, two or three massive 
neutrino species with a common mass. In this case, we have seen that the 
position of the break in the power spectrum P{k) is related to /Cnr, i-e. to 
the time of the non-relativistic transition which depends on the individual 
mass rriu. Instead, the amplitude of the break is related both to /j, and Onr, 
which depend respectively on and my. However, the results in Fig. 12 
present a much stronger dependence on fy than on Ny-. thus like for the CMB, 
the dominant effect is that of the total mass, but the individual masses are 
not completely irrelevant. In conclusion, we expect various scenarios with the 
same total mass My and different numbers of massive species Ny — 1, 2, 3 to 
share the same CMB spectrum, but not the same matter power spectrum, 
with small variations in the scale and amplitude of the damping caused by 
neutrino free-streaming. For the mass schemes shown in Fig. 1, the normal 
hierarchy (NH) and the inverted hierarchy (IH), the differences are expected 
to be even more subtle, because there are three distinct values of m,^, and 
h 

This problem has been studied numerically in Ref. [101], using a modified 



64 




Fig. 16. (Left) Evolution of the total neutrino energy density as a function of the 
scale factor of the Universe for models where the same total mass Mj, = 0.12 eV 
is distributed among the three species in various ways. Each line corresponds to 
the energy density in one of four different cases (only one massive eigenstates (1), 
two degenerate massive states (2), normal hierarchy (NH), inverted hierarchy (IH)) 
normalized to the case with three degenerate massive states. (Right) Comparison 
of the matter power spectrum obtained for the same models, divided each time by 
that with three degenerate massive states. Differences in the mass splitting affect 
the position and amplitude of the break in the power spectrum. 



version of CMBFAST with three independent neutrino sectors^ . Fig. 16 is taken 
from this reference. First, for a fixed total neutrino mass Mi, = 0.12 eV, we 
see on the left plot the ratio of the total neutrino density Pu{a) in the NH, 
IH, A^j, = 1 or degenerate — 2 scenario over that of the degenerate N^y — 3 
scenario. The bump, more or less pronounced depending on the case, is related 
to the different epochs of transition to the non-relativistic regime. In the right 
plot, we can see the ratio between the matter power spectrum P{k) in the 
same cases over that of the degenerate N,^ = 3 scenario. The bump around 
k — 10~^ /iMpc~^ is caused by the shift in the break scale, and the large-A; 
asymptotic values differ from one because of the different break amplitudes. 
Generally speaking, the NH scenario interpolates between the = 1 and the 
degenerate A^^ = 3 cases, depending on the value of the lightest neutrino; the 
IH case interpolates between the degenerate Ny = 2 and Ny = 3 cases. But the 
difference between the NH and IH power spectra always remain very small: 
for instance, in the case shown in Fig. 16 with = 0.12 eV, the two P{k) 
differ at most by 0.1%. It is very unlikely that this effect can be measured one 
day (see however section 6). 



^ See also the recent discussion in [102], as well as the study in [103], where CAMB 
was modified in order to perform calculations with two different neutrino masses. 
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Fig. 17. Shape distortion of the matter power spectrum P{k, z) at different redshifts. 
Upper solid curves: the power spectra are divided by [ag(a)]^ = [g{z)/{l + z)]"^, 
which accounts for the evolution on large scales. From top to bottom, the curves 
correspond to z = 10,8,6,5,4,3,2,1,0. Lower dashed curves: the power spectra 
are divided instead by [a g{a)]'^~^^^ ■^'^ = [g{z)/{l + 2)]^^^/^-^", which approximates 
the power spectrum evolution on small scales. From top to bottom, the curves 
correspond to z = 0, 1, 2, 3, 4, 5, 6, 8, 10. All cosmological parameters are kept fixed 
and the neutrino density fraction is fi, = 0.1. 

4.6.4 Massive neutrinos and redshift dependence of the matter power spec- 
trum 

In absence of neutrinos, we have already seen that the shape of the power 
spectrum during matter and A domination is redshift-independent for modes 
well inside the Hubble radius, while the global power spectrum amplitude 
grows like [0(7(0)]^, where g{a) is a damping factor normalized to g = 1 for 
a a\ (here a\ is the scale factor at matter/A equality). In terms of redshift 
evolution, this gives 

Pik, zy'^=' = ( (Y^l^) ^(^' 0)^^=° ' ^ » «^ • (144) 

The reason for this shape invariance is simply that the matter fluctuations 
obey to a scale-independent equation of evolution (see for instance Eq. (128), 
which can easily be generalized to the case of A domination). 

In presence of massive neutrinos, the shape of the power spectrum is not in- 
variant with redshift. Indeed, we are effectively dealing with a two-component 
system (baryons and CDM on the one hand, neutrinos on the other) in which 
the largest scales satisfy 6^ = Scdm [ag{a)], while for the smallest scales 
^ 5cdm oc [a g{a)Y~^^^ , as explained in Sec. 4.5.6. Therefore, the power 
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spectrum has a different redshift dependence on large and small scales, obeying 



P{k,z) 



I V(i + %(o). 



P{k,0) 



for aH < k < kjir , 



2-6/5/^ 



(145) 



P{k, 0) for 



k:^ k„ 



We illustrate this behavior in Fig. 17. Had we plotted the ratio P{k, z) / P{k, 0) , 
we would have seen both the shape and amplitude evolution of the matter 
power spectrum. Instead, we want to visualize only the shape distortion. So, 
we divide all power spectra either by [g(z)/(l + z)]'^ or [g(z)/{l + z)]'^~^^^ ■I''' . 
As expected, in the first case, the largest scales are corrected for their time 
evolution, and the ratio remains equal to one for k < 0.01 /iMpc^; in the 
second case, the smallest scales are corrected for their time-evolution and the 
ratio is equal to one for k > 2h Mpc~^. Between these two limits, the neutrino 
perturbations are just about to reach the Si, — > Scdm asymptotic behavior 
around the present time (see Sec. 4.5.7). For fi, = 0.1, the maximal shape 
distortion of the matter power spectrum is as large as 6% for z — 1 and 15% 
for z — 3. 



5 Current observations and bounds 



In this section we review how the available cosmological data can be used 
to get information on the absolute scale of neutrino masses, complementary 
to laboratory experiments. After a brief introduction to the statistical meth- 
ods used in cosmology (in particular, Bayesian inference), we describe the 
importance of each set or combination of cosmological data with respect to 
neutrino masses, and explain how the bounds depend on the data used and 
the underlying cosmological model. We summarize the current limits in Tables 
1-3, where the upper bounds on neutrino masses are given at 95% CL. after 
marginalization over all free cosmological parameters. We briefly discuss how 
the current bounds on neutrino masses can be relaxed or even avoided when 
neutrinos possess non-standard properties. 



5. 1 Statistical methods for parameter inference from cosmological data 



The theory of cosmological perturbations is a stochastic theory. For a given 
model, corresponding to some parameter values 6i, the theory predicts the 
probability of any particular realization, i.e. of any Universe with given per- 
turbations. Therefore, it can also predict the probability of obtaining a given 
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data set. The likelihood function jC{xn\Oi) is defined as the probability of ob- 
serving the data a;„ given the model parameters 9i. Such a function can be 
built for any data which has been properly modeled, i.e. for which the instru- 
mental noise properties have been fully understood, and the level of correlation 
between different data points has been estimated. 

The likelihood function can be used as a tool for solving the inverse problem, 
i.e. measuring the goodness-of-fit of a model given the data, or deriving con- 
fidence hmits on the model parameters 9i. This highly non-trivial task can be 
performed in different ways; most particle physicists follow the so-called fre- 
quentist approach, while until now, cosmologists usually prefer the Bayesian 
approach (with very few exceptions, see e.g. [104,105,106]). None of these 
methods is better in absolute; they are not directly comparable, because they 
address different questions. The selection of one method rather than the other 
should depend of the particular goal of the analysis, on how much the theory 
is constrained and understood, and on how much we trust the data modeling 
(for a detailed comparison of the two methods, see [107] or [108] and references 
therein) . 

The neutrino mass bounds derived in the next subsections are all based on the 
Bayesian inference method, which has presumably no secrets for readers with 
some experience in cosmology. On the other hand, many particle physicists 
are so familiar with the frequentist philosophy that cosmological bounds leave 
them perplex. In the rest of this subsection, we briefiy summarize the method 
used for deriving Bayesian confidence limits. Then, we emphasize some prop- 
erties of Bayesian inference which are radically different from their frequentist 
counterpart. 

In any Bayesian parameter extraction, the first step is to assume a theoretical 
model, with a set of N free parameters 9i, and N probability distributions 
n(^j) representing our knowledge on each parameter a priori, before using 
the data. These priors can be simply fiat distributions over the interval in 
which the parameter is meaningful. For instance, in the case of the neutrino 
fraction f^, it seems reasonable to define n(/,^) = 1 for < /j, < 1 and 
n(/,^) = otherwise: such a prior is called a top-hat. Priors can also contain 
some information from external data sets. For instance, if one of the free 
parameters is the current expansion rate Hq, one can choose to derive its 
prior from the results of the Hubble Space Telescope (HST) Key Project, and 
define II{Hq) as a Gaussian probability with central value 72 kms~^Mpc~^ 
and standard deviation u = 8 kms~^ Mpc~^ [109]. 

Given the likelihood C{xn\6i) and the priors Il{6i), one can compute the prob- 
ability distribution of any parameter 6j, or the joint probability for any subset 
of parameters. Such a distribution is called a posterior, since it represents the 
parameter probability a posteriori, after including the full data. Posterior dis- 
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Fig. 18. Situations which are possible in Bayesian inference, while they would be 
paradoxical in a frequentist analysis. (Left) Adding more data does not always lead 
to stronger bounds. We consider here a model with one free parameter 6 with flat 
prior. A first data set has a Gaussian likelihood Li{6) which predicts 6 = 0.3 it 0.1 
at the 68% confidence level (this interval corresponds to the solid vertical lines). A 
second data set has a non-Gaussian likelihood L2{6). The product Li{6) x L2{0) 
leads to a wider allowed range, corresponding to the dotted vertical lines. (Right) 
Adding more free parameters does not always lead to weaker bounds. We study here 
a two-dimensional model with likelihood L{6i,62) and flat priors. If there are good 
physical motivations to fix ^2 = 1, one can also study the one-dimensional model 
with likelihood L{6i, 1) (the thick line on the plot). We choose the likelihood in such 
way that J d9iL{6i, 1) = 1: the thick line can be directly interpreted as the posterior 
distribution for 9i in the one-dimensional model. For the two-dimensional model, the 
dashed line shows the posterior P{9i) obtained by marginalizing over 62, and nor- 
malizing to one. Given the very specific form of the likelihood, the two-dimensional 
model leads to a narrower posterior and to stronger bounds on 9i. 

tributions are obtained by marginalizing over the remaining free parameters. 
In practice, this amounts in integrating the product of the likelihood times 
the priors over these parameters. The result is a function of one or several 
parameters 9j, which can be normalized to one. One finally obtains a proba- 
bility distribution V{6j), from which it is straightforward to derive confidence 
limits: the n% confidence interval / is defined as the region in which 

(1) the cumulative probability is equal to n/100: Jjd9jV{9j) = n/100, and 

(2) the probability is maximal: the value of V{9j) for any point inside / is 
larger than that for any point outside /. 

The Bayesian inference method is well-defined, efficient, and has several ad- 
vantages. However, Bayesian results cannot be interpreted in the same way as 
frequentist ones, and there are significant differences that should be kept in 
mind: 

• Bayesian inference always produce some confidence intervals, even if the 
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assumed model is unable to provide good fits to the data. The results of a 
Bayesian analysis contain no statements on the goodness-of-fit of the model. 
Therefore, most authors provide some complementary information, like the 
maximum value of the likelihood £max (or equivalently, the minimum value 
of the effective x^, which is defined - even for a non-Gaussian likehhood - 
as XeS = — 21n£). Then, it is possible to estimate in the frequentist way 
how good is the best fit ^ . 

• for each parameter, the best-fit value 9'^ generally differs from the most likely 
value 91- Indeed, 9f corresponds to the best-fit model, while 9l maximizes 
the marginalization integral. They can be quite different, for instance if for 
9i fixed to 9^ there are many bad fits and very few excellent fits, while for 
9i fixed to 9l there are many good ones. In a future analysis, it would be 
conceivable that the best- fit value of the total neutrino mass be zero, while 
the most likely value be strictly positive. This is a natural outcome of the 
Bayesian approach, which is based on probability distributions in parame- 
ter space, and therefore favours regions in which models are good in average. 

• the exact meaning of a Bayesian confidence interval is: "given a particular 
underlying model and given the data, it is likely to be in this region, and un- 
likely to be outside" . One should be careful not to over-interpret a Bayesian 
confidence interval. It is risky to say that a parameter value sitting outside 
of the interval is cxchidcd. First, it is possible that with such a value one 
can find a very good fit to the data (in the frequentist sense): the Bayesian 
analysis could penalize this point just because the well-fitting region has a 
tiny volume in the subspace of parameters to be marginalized over. Second, 
one could possibly assume another model (with extra free parameters) in 
which the same parameter value would be comfortably allowed. In summary, 
Bayesian confidence intervals are well-defined, but they are attached to a 
particular underlying model. It makes no sense to quote a Bayesian upper 
limit, e.g. on the neutrino masses, without explicitly specifying the model 
assumed. 

• in the frequentist approach, adding more data always lead to stronger 
bounds on each parameter, while adding more free parameters can only 
enlarge the allowed region. In the Bayesian approach, this is true in general 
but not always: it is easy to build counter-examples, as we did in Fig. 18. 

Again, Bayesian inference offers many advantages for cosmological parameter 
extraction, and the previous items should not be regarded as criticisms, but 



However, comparing the merits of two models is a non-trivial task. A simple com- 
parison of the two maximum likelihood values is insufficient. Following the Bayesian 
philosophy, the most rigorous approach consists in computing the so-called Bayesian 
Evidence [110,111]. 
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rather as a list of warnings for correctly interpreting the results quoted in this 
section. 

5.2 CMB anisotropies 

The experimental situation of the measurement of the CMB anisotropies 
changed radically with the first year release of WMAP data (WMAPl), that 
includes the temperature x temperature (TT) [112] and temperature x E- 
polarization (TE) [113] correlation functions. The recent third-year release 
(WMAP3) improves the TT and TE power spectra [114], and adds a detection 
of the E-polarization self-correlation spectrum (EE) [115]. The WMAP exper- 
iment measured with high precision the TT spectrum on large angular scales 
(/ < 900), while at similar or smaller scales we have results from experiments 
that are either ground-based such as the Arcminute Cosmology Bolometer Ar- 
ray Receiver (ACBAR) [116], the Very Small Array (VSA) [117], the Cosmic 
Background Imager (CBI) [118] and the Degree Angular Scale Interferome- 
ter (DASI) [119], or balloon-borne such as ARCHEOPS [120], BOOMERANG 
(TT [121], TE [122] and EE [123]) and MAXIMA [124]. When using data from 
different CMB experiments, one should take into account that they overlap in 
some multipole region, and not all data are uncorrelated. 

Small neutrino masses only have a small influence on the power spectrum 
of CMB anisotropies, since for rrii, smaller than approximately 0.5 eV the 
neutrinos are still relativistic at the time of recombination, so that their direct 
effect on the CMB perturbations is identical to that of massless neutrinos. 
However, the signature of a mass smaller than 0.5 eV does not vanish: we 
recalled in Sec. 4.6 that there is always a background effect, proportional to 
fl^, which changes some characteristic times and scales in the evolution of the 
Universe, and affects mainly the amplitude of the first acoustic peak as well 
as the location of all the peaks, as shown in Figs. 14 and 15. Therefore, it is 
possible to constrain neutrino masses using CMB experiments only, down to 
the level at which this background effect is masked by instrumental noise, or by 
cosmic variance, or by parameter degeneracies in the case of some cosmological 
models beyond the minimal A Mixed Dark matter (AMDM) framework. 

We list in Table 1 the bounds on the sum of neutrino masses that were obtained 
in analyses of CMB data after the first WMAP results appeared. In these 
works, a minimal AMDM scenario was assumed: the total neutrino mass was 
the only additional parameter with respect to a flat ACDM cosmological model 
characterized by 6 parameters'^' . However, there were some differences in the 
analysis procedure. 

^ This will be the case for the bounds reviewed in this section, unless specified 
otherwise. 
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Table 1 

Upper bounds on Mj^ (eV, 95% CL) from recent analyses of CMB data only. Note 
that the results of Refs. [127,129,130,131] are consistent with each other (see the 
text for details). 



Ref. 


Bound 


Data used 


[125] 


10.6 


WMAPl 


[127] 


2.0 


WMAPl 


[65] 


2.0 


WMAP3 


[129] 


2.1 


WMAPl, VSA, ACBAR & CBI 


[130] 


1.6 


WMAPl, VSA, ACBAR & CBI 


[131] 


3.1 


WMAPl, BOOMERANG03, VSA, ACBAR, 



CBI. DASI k MAXIMA 



The analysis in Ref. [125] used a Monte Carlo Markov Chain (MCMC) method 
for the minimization, concluding that WMAP data are consistent with neutri- 
nos making 100% of the dark matter component (see also [126]), thus giving 
the loose bound on the total neutrino mass M,y < 10.6 cV. This would mean 
that eV masses are basically indistinguishable from CDM at the epoch of last 
scattering, and therefore have little impact on the CMB scales measured by 
WMAP. 

However, a very different result was discussed in ref. [127] (see also [128]), 
where an upper limit of 2 eV was found from the analysis of WMAPl only 
using a grid-based method. The authors of this last analysis ascribe the dif- 
ference with respect to [125] to a better marginalization procedure. However, 
using the WMAPl data and an MCMC method, we find the same result as in 
[127]. Finally, the 2 eV bound is now confirmed by the analysis of the WMAP3 
data in Ref. [65] . 

In Table 1, we also show the results from two recent works [129,131] based 
on the combination of WMAPl with other recent CMB experiments. We also 
show the result of a very recent analysis [130] in a cosmological model with 11 
parameters (similar to that in Table 4 of Ref. [132], but in a flat Universe). The 
2cr bound obtained by [129,130] is still close to 2 eV. The result of [131], which 
include more data from BOOMERANG, DASI and MAXIMA, is larger by 
50%. We saw in Sec. 5.1 that in a Bayesian analysis, it is possible to add more 
data and obtain weaker bounds. In the present careful examination 

of Fig. 8 in [131] shows that their analysis based on WMAPl only is in very 
good agreement with [127,129]; however, at least one of the DASI or MAXIMA 
data set tends to push the bound upwards, because it is better fitted with a 
non-zero neutrino mass. 
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In conclusion, many analyses support the conclusion that a sensible bound on 
neutrino masses exists, of order of 2 — 3 eV at 2a. This is an important result, 
since it does not depend on the uncertainties from LSS data (see Sec. 5.3). 
The only way to relax this bound would be to complicate the cosmological 
scenario assumed for parameter extraction, but this issue will be discussed in 
Sees. 5.6 and 5.7. 

However, these results do not take advantage of the most significant observable 
signature of massive neutrinos: the suppression of small-scale matter fluctua- 
tions induced by neutrino free-streaming. In order to get improved bounds on 
neutrino masses it is crucial to combine CMB data with independent large- 
scale structure observations. 

5.3 Galaxy redshift surveys 

We have seen in Sec. 4.6 that free-streaming of massive neutrinos produces a 
direct effect on the formation of cosmological structures. As shown in Fig. 13, 
the presence of neutrino masses leads to an attenuation of the linear matter 
power spectrum on small scales. 

In a seminal paper, Hu and collaborators [94] showed that an efficient way 
to probe neutrino masses of order eV was to use data from large redshift 
surveys, which measure the distance to a large number of galaxies, giving us 
a three-dimensional picture of the universe (modulo the effect of unknown 
peculiar velocities on top of the Hubble flow, which can be corrected only 
statistically). At present, we have data from two large projects: the 2 degree 
Field (2dF) galaxy redshift survey [133,134], whose final results were obtained 
from more than 220,000 galaxy redshifts, and the Sloan Digital Sky Survey 
(SDSS) [135,89], which will be completed soon with data from one million 
galaxies ^ . 

One of the main goals of galaxy redshift surveys is to reconstruct the power 
spectrum of matter fluctuations on very large scales, whose cosmological evo- 
lution is described entirely by linear perturbation theory. However, the linear 
power spectrum must be reconstructed from individual galaxies, i.e. from ob- 
jects which underwent a strongly non-linear evolution. A simple analytic model 
of structure formation suggests that on large scales, the galaxy-galaxy corre- 
lation function should be, not equal, but proportional to the linear matter 
density power spectrum, up to a constant factor that is called the light-to- 
mass bias. This parameter, denoted by b, can be obtained from independent 
methods (e.g. high-order statistics beyond the two-point correlation function). 

^ An extension of this experiment, SDSS-II was recently funded (see the SDSS 
webpage http : / / www . sdss . org ) 
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Table 2 

Upper bounds on My (eV, 95% CL) from recent analyses combining CMB, galaxy 
redshift surveys (shape of the matter power spectrum only) and other data. 



Ref. 


Bound 


Data used fin addition to WMAPl excent for Eef f65l^ 


[561 


1.2 


other CMB (pre-WMAP), 2dF-gal 




1.0 


Drcvious + HST SNIa 


[1251 


1.74 


SDSS-g-al 


flSQl 
L±o»j 


n 75 


other CMB fnrp-WMAPl 2dF-e-al SDSS-e-al HST 


flOfll 

[lUUJ 


1.0 


ACBAR 2dF-ffal SDSS-eal 




0.6 


Drevious + HST SNIa 




Q6 


VSA 2dF-B-al 

V ^Li-L ftCt-L 


fl37l 

[XO ( J 


1.54 


SDSS-eal SNIa 


[140] 


1.4 


Other CMB , 2dF-gal, HST, SNIa 


[131] 


1.2 


other CMB , 2dF-gal, SDSS-gal 


[129] 


1.27 


other CMB , SDSS-gal 




1.16 


other CMB , 2dF-gal 


[65] 


0.87 


\\'MAP3. 2dF-gal 



Currently, these methods tend to confirm that the linear biasing assumption 
is correct, at least in first approximation. 

A conservative way to use the measurements of galaxy-galaxy correlations in 
an analysis of cosmological data is to take the bias as a free parameter, i.e. to 
consider only the shape of the matter power spectrum at the corresponding 
scales and not its amplitude. In this case, we will refer to the corresponding 
data of the 2dF and SDSS surveys as 2dF-gal and SDSS-gal, respectively. 
Alternatively, one can add a measurement of b, as recently done in three anal- 
yses. The value b = 1.04 ± 0.11(1(7) [2dF-bias] [136] from the 2dF bispectrum 
analysis was used in [55], while [137] adopted the estimate b — 0.99±0.11(lcr) 
[SDSS-bias] from galaxy-galaxy lensing from SDSS [138]. Finally, in Ref. [131] 
the intermediate estimate 6 = 1.0 ± 0.10(1(7) was used. 

We summarize in Table 2 the bounds on neutrino masses that were obtained 
in combined analyses of CMB and LSS data, using galaxy clustering data in 
the range of wavenumbers (0.015 - 0.02) h/Mpc < k < {0.15 - 0.20) h/Mpc 
(including non-linear corrections for the smallest wavelengths [141]). The re- 
sults in this Table were obtained leaving the bias as a free parameter. One 
concludes from Table 2 that the addition of LSS data improves the bounds 
with respect to the numbers in Table 1, with an upper limit on between 
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0.9 and 1.7 eV for SDSS-gal and/or 2dF-gal added to CMB data. In general, 
one obtains weaker bounds on neutrino masses using preliminary SDSS results 
instead of 2dF data, but this conclusion is expected to change after the next 
SDSS releases. The improvement with respect to previous results is evident: 
in 2001, ref. [142] found < 4.2 eV (95% CL) from the combination of the 
PSCz survey [143] with pre-WMAP CMB data, while the analysis in [144] 
used early 2dF data to obtain Mi, < 2.2 eV with priors on the values of Jim 
and h. 

There exist other measurements of cosmological parameters that are often 
added to CMB and LSS data, as shown in Table 2. For instance, the present 
value of the Hubble parameter was measured by the Key Project of the Hub- 
ble Space Telescope (HST), giving h = 0.72 ± 0.08 (la) [109], which excludes 
low values of h and leads to a stronger upper bound on the total neutrino 
mass. In addition, one can include the constraints on the current density of 
the dark energy component deduced from the redshift dependence of type la 
supernovae (SNIa) luminosity, which measures the late evolution of the ex- 
pansion rate of the Universe. For a flat Universe with a cosmological constant, 
these constraints can be translated into bounds for the matter density fl^, 
which range from the conservative value fim = 0.28 ± 0.14 (lo") from [145] to 
the more recent fl^ = 0.29l[j:[j| (la) [146] from the analysis of SNIa with the 
HST (other recent works include [147,148]). 

As shown in Table 3, the bounds on neutrino masses are more stringent when 
the amplitude of the matter power spectrum is fixed with a measurement of 
the bias [55,131,137], instead of leaving it as a free parameter. The upper 
limits on Mi, are reduced to values of order 0.5 — 0.9 eV (95% CL), although 
some analyses also add Lyman-a data (see next subsection). 

Finally, a galaxy redshift survey performed in a large volume can also be 
sensitive to the imprint created by the baryon acoustic oscillations (BAO) at 
large scales on the power spectrum of non-relativistic matter, that we discussed 
in Sec. 4.4 (see Fig. 7). Since baryons are only a subdominant component of 
the non-relativistic matter, the BAO feature is manifested as a small single 
peak in the galaxy correlation function in real space that was recently detected 
from the analysis of the SDSS luminous red galaxy (LRG) sample [89] at a 
separation of 100 h^^ Mpc. The observed position of this baryon oscillation 
peak provides a way to measure the angular diameter distance out to the 
typical LRG redshift of 2; = 0.35, which in turn can be used to constrain the 
parameters of the underlying cosmological model. The SDSS measurement of 
the angular diameter distance at z = 0.35 [SDSS-BAO] was included in the 
analysis of Ref. [152] to get a bound of 0.44 eV (95% CL) on the total neutrino 
mass M„, as shown in Table 3. 
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Table 3 

Upper bounds on M^, (eV, 95% CL) from recent analyses combining CMB, LSS 
(including bias and/or Lya) and other cosmological data. 



Ref. 




Dfltfl used (\r\ additinn tn WMAPI p-xcmt for TRef [fi,5l^ 


[551 


0.63 


ACBAR, CBI, 2dF-gal, 2dF-bias 




0.68 


"niT'vimm -1- T/vrv 1 4Q 1 P^D 




0.54 


SDSS-g-al SNIa SDSS-Lvo' 




0.42 


previous + SDSS-bias 


fl40l 


0.47 


other CMB 2dF-eal HST SNIa SDSS-Lva 


fl511 


0.65 


SDSS-ffal HST Lva fl49 1501 


[131] 


0.48 


Other CMB, 2dF-gal, SDSS-gal, bias 6 = 1.0 ± 0.10 


[152] 


0.44 


other CMB, 2dF-gal, SDSS-gal, HST, SNIa, SDSS-BAO 




0.30 


previous + SDSS-Lya 


[65] 


0.91 


WMAP3, SDSS-gal, SDSS-bias 




0.68 


previous -|- other CMB, 2dF-gal, SNIa 



5.4 Lyman-a forest 

The matter power spectrum on small scales can also be inferred from data on 
the so-called Lyman-a forest. This corresponds to the Lyman-a absorption of 
photons traveling from distant quasars (2; ~ 2 — 3) by the neutral hydrogen 
in the intergalactic medium. As an effect of the Universe expansion, photons 
arc continuously red-shifted along the line of sight, and can be absorbed when 
they reach a wavelength of 1216 A in the rest-frame of the intervening medium. 
Therefore, the quasar spectrum contains a series of absorption lines, whose 
amplitude as a function of wavelength traces back the density and temperature 
fluctuations of neutral hydrogen along the line of sight. It is then possible to 
infer the matter density fluctuations in the linear or quasi-linear regime (see 
e.g. [153,154,155,156]). 

In order to use the Lyman-a forest data, one needs to recover the matter power 
spectrum from the spectrum of the transmitted flux, a task that requires the 
use of hydro-dynamical simulations for the corresponding cosmological model. 
This matter power spectrum is again sensitive to the suppression of growth 
of mass fluctuations caused by massive neutrinos. In 1999, Ref. [157] obtained 
an upper bound of rui, < 5.5 eV (for one massive neutrino species) from the 
combined analysis of Lyman-a data and other cosmological measurements. 

Given the various systematics involved in the analysis pipeline, the robustness 
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of Lyman-o; forest data is still a subject of intense discTission between ex- 
perts. Among these systcmatics, we can cite the uncertainty existing on: the 
mean absorption by the intergalactic medium (the effective optical depth), the 
temperature-density relation within the absorbing gas, the effect of galactic 
winds, or the proper use and accuracy of hydro-dynamical simulations in the 
mildly non-linear regime. Different authors tend to describe or to deal with 
these uncertainties in slightly different ways [149,150,158,159], an indication 
that the field might not be fully mature at the moment. In addition, the way in 
which neutrinos are implemented in hydro-dynamical simulations is generally 
rather simplistic. 

Some of the results shown in Table 3 were obtained adding Lyman-a data 
to CMB and other LSS data. For a free bias, one finds that Lyman-a data 
help to reduce the upper bounds on the total neutrino mass: using the SDSS 
Lyman-a forest [159,160], Ref. [137] finds < 0.54 eV (95% CL), while 
Ref. [151] quotes 0.65 eV using Lyman-a data from [149,150]. At least two 
recent analyses combined results on the bias and the Lyman-a forest: in Ref. 
[137] the addition of SDSS-bias allowed to reduced the upper bound on M^, 
to 0.42 eV, while the first year analysis of WMAP collaboration [55] found a 
slightly weaker bound of 0.68 eV when adding 2dF-bias and Lyman-a data 
from [149,150]. Finally, Ref. [152] found the upper bound < 0.30 eV (95% 
CL) adding simultaneously SDSS-Lya and SDSS-BAO data. 

Thus, at present, the impact of Lyman-o; data on the neutrino mass bounds 
is comparable to that of galaxy bias measurements, if not weaker. We can 
conclude that both the Lyman-a data and the light-to-mass bias determination 
are powerful tools for getting information on neutrino masses, but due to 
the many involved systematics, the current bounds obtained with them are 
probably not as robust as those in Tables 1 and 2. 



5.5 Summary of current bounds 

The results shown in Tables 1-3 are representative of an important fact: a 
single cosmological bound on neutrino masses does not exist. A graphical 
summary is presented in Fig. 19, where the cosmological bounds correspond 
to the vertical bands, which were grouped according to the ranges quoted in 
Tables 1, 2 and 3, respectively. The thickness of these bands roughly describe 
the spread of values obtained from similar cosmological data: 2 — 3 eV for 
CMB only, 0.9 - 1.7 eV for CMB and 2dF/SDSS-gal or 0.3 - 0.9 eV with the 
inclusion of a measurement of the bias and/or Lyman-a forest data and/or 
the SDSS measurement of the baryon oscillation peak. 

One can see from Fig. 19 that current cosmological data probe the region of 
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Present bounds 




lightest rriy (eV) 

Fig. 19. Current upper bounds (95%CL) from cosmological data on the sum of 
neutrino masses, compared to the values in agreement at a 3a level with neutrino 
oscillation data in Eq. (2). 

neutrino masses where the 3 neutrino states are degenerate, with a mass M^/S. 
This mass region is conservatively bounded to values below approximately 1 eV 
from CMB results combined only with galaxy clustering data from 2dF and/or 
SDSS (i.e. the shape of the matter power spectrum for the relevant scales). The 
addition of further data leads to an improvement of the bounds, which reach 
the lowest values when data from Lyman-a and/or the SDSS measurement of 
the baryon oscillation peak are included or the bias is fixed. In such cases the 
contribution of a total neutrino mass of the order 0.3-0.6 eV seems already 
disfavoured, in particular covering the recent positive 0i>2f3 signal from [31] 
(of course with the corresponding caveats, see e.g. the discussion in [14]). 

There exist other measurements of cosmological parameters that were not 
included in the analyses we have reviewed so far, although they could play 
a role for neutrino mass determination. In particular, we have seen that any 
measurement that fixes the amplitude of the matter power spectrum, such 
as the bias, would be very important. For instance, the parameter ag -which 
gives the linear rms mass fluctuations in spheres of radius 8 h~^Mpc- can be 
measured, for instance, from the determination of cluster masses through X- 
ray observations (see e.g. [161]) or from weak lensing (see e.g. [162]). If various 
methods happened to converge towards an accurate determination of ag, the 
determination of neutrino masses would improve significantly, as discussed in 
Ref. [163]. Since erg probes the amplitude of matter fluctuations on scales of 
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the order a few Mpc which are deeply affected by neutrino free-streaming, 
a large value of would push down the upper bound on M^, while a low 
value would potentially bring evidence for non-zero masses. However, at the 
moment, the ensemble of erg determinations is rather spread (see Table 5 in 
[125]), suggesting that not all systematics errors are properly modelled. 

One of the first analysis of neutrino mass bounds from cosmology used data 
on (Ts from cluster abundance to obtain an upper bound of order a few eV 
[164]. Another possibility was recently explored in [165], using data on the 
redshift-evolution of the number density of massive galaxy clusters to find the 
bound < 2.4 eV (95% C.L.). Ref. [166] used another probe of the small- 
scale power spectrum amplitude from X-ray cluster data, and found some 
marginal evidence for a non-zero neutrino total mass, Mj, = 0.56lo;26 eV at 
68% C.L. This result is related to the fact that the value of as preferred by 
these X-ray data is as small as to contradict most other determinations of 
this parameter, for instance from weak Icnsing. This example shows that the 
sensitivity of cosmological observations to neutrino masses is a powerful tool, 
but its implications should not be extracted without care. 



5. 6 Eodra parameters 

Assuming that the relic neutrinos are standard, we have seen that the limits 
mainly depend on the cosmological data used in the analysis, although minor 
differences on the quoted bounds exist even when using similar data. These 
differences arise because not all the analyses in the Tables used the same set of 
cosmological parameters and priors (such as the assumption of a flat Universe). 
It is thus important to test whether the impressive cosmological bounds on 
neutrino masses change much if additional cosmological parameters, beyond 
those included in the minimal ACDM, are allowed. This could be the case 
whenever a new parameter degeneracy with the neutrino masses arises. 

An interesting case of degeneracy between cosmological parameters is that be- 
tween neutrino masses and the radiation content of the universe (parametrized 
via the effective number of neutrinos A^eflf)- The extra radiation partially com- 
pensates the effect of neutrino masses, provided that other cosmological pa- 
rameters such as fijn and h arc varied, leading to a less stringent bound on 
Mj, [56,100,167,168]. This parameter degeneracy is shown in Fig. 20, taken 
from Ref. [100]. For instance, when one considers a model with four instead 
of three species of massive neutrinos, the upper bound found from CMB and 
2dF/SDSS-gal is relaxed from 0.8 to 1.2 eV when one of the states is much 
heavier than the others. These results are interesting for the 4-neutrino mass 
schemes that also incorporate the results of the LSND experiment [19], that 
we discussed in Sec. 2. At present, the LSND regions in the space of oscillation 
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Fig. 20. Two-dimensional likelihood in (A'efT, M^) space, marginalized over the other 
cosmological parameters of the model. We plot the la (green / dark) and 2a (yellow 
/ hght) allowed regions. Here we used CMB (WMAPl & ACBAR) and LSS (2dF 
& SDSS) data, adding extra priors on h (HST) and il^a (SN99 or SN03) as follows: 
(a) no priors, (b) HST, (c) HST+SN99, (d) HST+SN03. For details, see ref. [100]. 

parameters are not yet completely disfavoured by cosmological data. 

Recently a degeneracy between neutrino masses and the parameter w, that 
characterizes the equation of state of the dark energy component X {px = 
wpx), was studied in [170]. For constant values w < —1, in the so-called 
phantom energy regime (which could describe particular scenarios based on 
string theory or modified gravity), it was shown that the neutrino mass bound 
can be relaxed to < 1.48 eV (95% C.L.), compared with < 0.65 eV 
for fixed w = —1. Similar conclusions were obtained by the WMAP team (see 
Fig. 18 in [65]). Ref. [170] studied the effect of this parameter degeneracy on 
the analysis of future cosmological data (see also [86] for models with a time- 
varying equation of state of dark energy). Finally, Ref. [152] showed that data 

^ However, a very recent analysis [169] showed that in a 4-neutrino scenario, the 
heavier neutrino mass is bounded to values smaller than 0.55 eV (95% CL) when 
data from CMB, LSS and Lyman-a are included. 
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on baryon acoustic oscillations (see Sec. 5.3) helps in breaking the rrii, — w 
degeneracy. 

Finally, the neutrino mass bounds could certainly be relaxed by assuming 
that the primordial power spectrum has a non-trivial shape or Broken Scale 
Invariance (BSI). This eventuality cannot be excluded, for instance if a phase 
transition takes place during inflation. Such an assumption opens so many 
degrees of freedom in the parametrization of cosmological models that a sys- 
tematic study is almost impossible. However, some authors [171] used this 
possibility in combination with a non-zero neutrino mass for building a cos- 
mological scenario without any form of dark energy, that could fit most of the 
data available at the time of publication, except for SNIa data. This model 
is now ruled out by the combination of WMAP with the SDSS galaxy-galaxy 
correlation function, but it illustrates the fact that if one is ready to add more 
degrees of freedom in cosmological scenarios (in the case of Ref. [171], they 
were actually not so many), then there is most probably a way to evade the 
bounds presented in this review. 



5. 7 Non-standard relic neutrinos 

Finally, let us remind the reader that the cosmological implications of neutrino 
masses could be very different if the spectrum or evolution of the cosmic neu- 
trino background was non-standard. Our aim here is not to give an exhaustive 
list of all possibilities (many results on non-standard neutrinos were reviewed 
in [6]), but only to briefly summarize some recent works. 

For instance, the bounds on neutrino masses would be modified if relic neutri- 
nos presented a momentum spectrum significantly different from the equilib- 
rium Fermi-Dirac distribution of Eq. (6). For instance, neutrinos could violate 
the Pauli exclusion principle and obey Bose-Einstcin statistics with important 
cosmological and astrophysical implications [172,173]. Non-equilibrium neu- 
trino spectra can also be produced in low-reheating scenarios [174,175,176,177,178,179], 
from the decay of massive neutrinos into relativistic products (see e.g. [180,181]) 
or as the result of active-sterile neutrino oscillations after decouphng (see 
[182,183] and references therein). A recent analysis of cosmological data in- 
cluding large non-thermal corrections in the neutrino spectra was done in Ref. 
[184]. 

We conclude giving two examples of models where the cosmological bounds 

on neutrino masses that we have discussed can be evaded. In the first case, 
massive neutrinos are strongly coupled to a light scalar field, in such a way 
that they annihilate as soon as they become non-relativistic and do not affect 
the formation of LSS [185]. The viability of this kind of models depends on 
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the cosmological data included in the analysis and on the number of neutrino 
species strongly interacting (see the discussions in Refs. [186,187,188]). 

The second possibility, described in [189,190] and many other recent works (see 
e.g. the discussion in [191]), concerns neutrinos coupled to the dark energy so 
that the dark energy density is a function of the neutrino mass (Mass- Varying 
Neutrinos or MaVaNs). Then, neutrino masses are only relevant for the very 
recent stages of the Universe, so that no bounds can be placed from high- 
redshift cosmological data. Note, however, that this is not true for all MaVaNs 
scenarios, see for instance [192]. 



6 Future sensitivities and new experimental techniques 

If the characteristics of future experiments are known with some precision, 
it is possible to assume a "fiducial model", i.e. a cosmological model that 
would yield the best fit to future data, and to estimate the error bar on a 
particular parameter that will be obtained after marginalizing the hypothet- 
ical likelihood distribution over all the other free parameters (see Sec. 5.1). 
Technically, the simplest way to forecast this error is to compute a Fisher 
matrix, as explained in the next subsection. This technique has been widely 
used in the literature, for many different models and hypothetical datasets. 
Here we will focus on the results for a{M^), the forecast 68% CL error on 
the total neutrino mass, assuming various combinations of future observations 
(CMB anisotropics, galaxy redshift surveys, weak lensing surveys, . . . ) and 
various cosmological models with fiducial parameter values. Many papers de- 
rived some predictions for a{M^); we will summarize here the main results 
and see whether the expected sensitivities will be sufficient for detecting the 
total neutrino mass, even if it is close to the minimum values guaranteed by 
present data on flavour neutrino oscillations. 



6.1 Fisher matrix 



In a Bayesian analysis, the probability distribution in the space of cosmological 
parameters is inferred from the likelihood jC{xn\Oi) of the data Xn given a 
theoretical model described by a set of parameters 6i (see Sec. 5.1). If we 
assume that the maximum likelihood £max is obtained for a particular model 
9^ (the "flducial model"), and if we approximate the likelihood function as a 
multivariate Gaussian in the vicinity of the maximum, the symmetrical error 
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on each parameter can be inferred from the Fisher matrix, defined as 



F 



d9id9j 



(146) 
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The la (or 68% confidence) error bar on a particular parameter 6j assuming 
that all other parameters are fixed to the fiducial values 9^ would be given 
by {Fjj)^^/'^. However, the really interesting quantity is the error on 9j as- 
suming that the other parameters arc unknown (and simply marginalized out 
by integrating over the likelihood). In that case, the error can be shown to 
be [193] 



Geometrically, and in the multivariate Gaussian approximation, the Fisher 
matrix gives the coefficients of the ellipsoids Fij9i9j = f{c) corresponding to 
given confidence levels c in parameter space. Therefore, the study of its eigen- 
vectors and eigenvalues gives an idea of the linear parameter combinations 
with the best and poorest determination (the latter correspond to parameter 
degeneracies) . 

In many cases one refers to data points Xi {i = l,...,n) with a Gaussian 
probability distribution centered in zero 

C {x\9) = [(27r)"|C|]-V2e-i^"^^^ , (148) 



where the "data covariance matrix" C contains the variance expected for the 
data, C = (,x,x*), after a detailed modelling of all possible contributions. Typ- 
ically, the data covariance matrix will be the sum of the theoretical model 
prediction (for parameters 9) and of the noise variance estimated from the 
instrumental characteristics. In general, the data covariance matrix is non- 
diagonal and the data points are correlated with each other. A purely alge- 
braic calculation shows that the Gaussian likelihood of Eq. (148) inserted in 
Eq. (146) leads to the following expression for the Fisher matrix [193] 



Fij = - Trace 



d9i 09, 



(149) 



6.2 Future CMB experiments 



The quality of the data from the WMAP satellite [55,65], complemented by 
the results of other experiments at smaller angular scales such as those hsted 
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in section 5.2, has shown the importance of CMB data as a probe of cosmo- 
logical parameters. There are several other projects in operational or design 
stage, aimed mainly at improving the sensitivity towards small scales and/or 
polarization anisotropics. The neutrino mass is not the primary target of these 
experiments, because as we have already seen in Sec. 4.6, small masses only 
have a modest signature on the CMB. But still, we know that there is a back- 
ground effect proportional to uo^ which renders the primary CMB anisotropics 
sensitive -at least in principle- to very small values of the total neutrino 
mass. Instead, typical targets are the scalar tilt and the tensor-to-scalar ratio, 
in order to better constrain the predictions of inflationary theories. However, 
in order to measure with precision the neutrino masses, it will be useful to 
combine some large scale structure observations with the best possible CMB 
dataset in order to reduce as much as possible parameter degeneracies. In this 
section we review the Fisher matrix computation and expected sensitivity of 
future CMB experiments alone; in the next sections the CMB contribution to 
the Fisher matrix will be combined with that from other hypothetical obser- 
vations. 

The CMB temperature maps can be decomposed in multipoles 
AT 

— (^,0)=^aLllU^», (150) 

l,m 



and similarly the polarization map arc described by multipoles af^ and af^. 
The data covariance matrix C = {af^aff^^,) receives contributions from: 

• primary anisotropics, for which the expected variance is predicted by the 
cosmological pcrtiirbation theory: {a,f„fl^^i) = 5w5mm'C^-^ . For parity rea- 
sons the non- vanishing power spectra are Cf^, C^^, C^^^ and Cj^ . 

• secondary anisotropics imprinted at redshifts z -C 1100, like the late Inte- 
grated Sachs- Wolfe (ISW) effect and the weak lensing induced by nearby 
clusters, as well as various foreground contaminations. The late ISW effect 
is usually included in the theoretical power spectrum computation, while 
the lensing effect is small enough to give a neghgible contribution to C 
(excepted for the small-scale B-mode variance). Here we will neglect the 
lensing distortions, which will be the object of Sec. 6.4. Most of the fore- 
ground contamination is expected to be accurately removed from the raw 
data using a multi-frequency analysis (since most foreground signals do not 
have a blackbody distribution). 

• experimental noise, which is usually parametrized by two quantities: the 
angular resolution and the sensitivity per pixel (here, a pixel represents the 
smallest area in a map on which an independent measurement is performed) . 
Intuitively, a given angular resolution implies that the error-bar will grow 
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dramatically below a given scale, and the finite sensitivity per pixel implies 
that even in absence of a primary signal, the covariance matrix would get 
some non-zero contribution at a given level. Technically, the angular resolu- 
tion is expressed as an angle Of,, the Full- Width at Half-Maximum (FWHM) 
of the instrument beam. The sensitivity (or noise squared variance) per pixel 
(for X — T, E, B) is usTially given in units of a squared temperature, 
and the dimensionless sensitivity per pixel is obtained dividing by the 
square of Tq ~ 2.726 K. The solid angle of a pixel can be approximated by 
9l, so that the dimensionless sensitivity per steradians reads (6'^A^). For 
one instrument, the noise variance can be shown to be [194] 

{aLaL') = 5u'5^m'5^^' Nr , ^ {elK\)B^' , (151) 

where the Gaussian beam window function Bi = exp[— /(/ -|- l)dl / (16 In 2)] 
has the effect of cutting the sensitivity for / ^ 6^^. Note that for X 7^ X' 
the variance vanishes, since the noise in one map (for mode T, E or B) is 
expected to be uncorrelated with the noise in another map. Finally, most 
experiments consist of various detectors, for which a global sensitivity can 
be built by summing the inverse variances 

Detectors are generally grouped in channels^ i.e. ensembles of detectors 
working at a given frequency. An efficient foreground extraction requires 
a large number of channels, at the expense of decreasing the number of 
detectors per channel, and hence increasing the noise per channel. 

The theoretical and noise variance are usually assumed to be uncorrelated, 
which allows to express the covariance matrix in a rather simple way. If we 
represent the data as a vector 

^)im, where the indexes 

{1,171) run over all multipoles, the full data covariance matrix can be decom- 
posed in diagonal blocks representing the covariance (x^^x*^) for a given 

{l,m) 



Clm — Ci 



fjTE fjEE ^ ^EE 

Cf ^ + iVf ^ 



(153) 



The index m does not appear in this matrix, since all multipoles (/,m) with 
the same / share the same variance. If we consider the data as Gaussian (which 
is a very good approximation except for small-scale 5- mode multipoles), we 
can use Eq. (149). The first step of the calculation is to estimate the number 
of multipoles (/,m) with a fixed I. For full-sky experiments, m runs from —I 
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to / and this number is {21 + 1). For experiments probing a fraction /sky of the 
sky, the number of independent multipole measurement at a given / is close 
to {21 + l)/sky This leads to 

^CMB^_^(2/ + l)/3kyTrace 

^ 1=2 



09, 



89, 



(154) 



which can be further reduced to an expression of the type [195] 

'max Pin^^ 
PS"" - Y.Y.^iCov,)-l,^ , (165) 

1=2 a,l3 '^^i '^"j 



where a and (3 run over TT, EE, TE, BB and Cov^ is the "spectra covariance 
matrix" . For space reason we do not give here the full expression of Cov; (see 
e.g. Ref. [196]). As an example, the first coefficient of the matrix reads 

{Covi)tt,tt = T^T-^rw— + • (156) 

[ZL + ij/sky 



The factor {21 + l)/sky which appears in the denominator of Cov; accounts for 
the effect of cosmic variance: for large / values, an experiment makes many 
independent measurements (corresponding to several m's) and the variance 
in absence of noise is much smaller than the power spectrum C^^ of a single 
multipole a^. For small /, there are only few independent measurement and 
the variance is of the same order of magnitude as C^^ . 

Let us apply this machinery to concrete cases. In Table 4, we present the 
expected sensitivity of a non-exhaustive list of future CMB experiments. The 
design of projected experiments is subject to constant evolution, and this com- 
pilation should be regarded as indicative only. We did not take the numbers 
from some official documentation edited by each experiment, but from the 
papers cited below, which worked out some neutrino mass error forecasts on 
the basis of the numbers presented here. We refer to these papers for more 
details and references concerning the sensitivity parameters. 

After WMAP, we can expect some very interesting results from a large number 
of ground-based experiments, mapping the anisotropics only in a small regions 
of the sky but usually with excellent sensitivity and resolution. Most of these 
experiments are optimized for measuring E-polarization or even B-polarization 
anisotropics, which will still be poorly constrained after the completion of 
WMAP, and even after Planck in the case of B-polarization. Among these 
collaborations, we can cite the South Pole Telescope* (SPT, in construction) 

* http : //spt . uchicago . edu7| 
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Table 4 

Experimental parameters of a few CMB projects, used in the neutrino mass forecast 
papers mentioned in this section. Here /sky is the observed fraction of the sky, v is 
the channel frequency in GHz, measures the FWHM of the beam in arc-minutes 
and At- p is the square root of the sensitivity per pixel in /xK for temperature and 
polarization. 



Experiment 


j sky 


V 


Ob 


At 


Ap 


SPT [197] 


0.1 


217 


0.9' 


12 


17 


BICEP [198] 


0.03 


100 


60' 


0.33 


0.47 






150 


42' 


0.35 


0.49 


OUaD fl98l 


0.025 


100 


6.3' 


3.5 


5.0 






150 


4.2' 


4.6 


6.6 


BRAIN [198] 


0.03 


100 


50' 


0.23 


0.33 






150 


50' 


0.27 


0.38 






220 


50' 


0.40 


0.56 


ClOVER [1981 


0.018 


100 


15' 


0.19 


0.30 






143 


15' 


0.25 


0.35 






217 


lo 


0.55 


0.7G 


Planck [198] 


0.65 


30 


33' 


4.4 


6.2 






44 


23' 


6.5 


9.2 






70 


14' 


9.8 


13.9 






100 


9.5' 


6.8 


10.9 






143 


7.1' 


6.0 


11.4 






217 


5.0' 


13.1 


26.7 






353 


5.0' 


40.1 


81.2 






545 


5.0' 


401 


oo 






857 


5.0' 


18300 


oo 


SAMPAN [198] 


0.65 


100 


42' 


0.13 


0.18 






143 


30' 


0.16 


0.22 






217 


20' 


0.26 


0.37 


CMBpol [197] 


0.65 


217 


3.0' 


1 


1.4 


Inflation Probe 


0.65 


70 


6.0' 


0.29 


0.41 


(hypothetical) [198] 




100 


4.2' 


0.42 


0.59 






150 


2.8' 


0.63 


0.88 






220 


1.9' 


0.92 


1.30 
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[199], the Atacama Cosmology Telescope^ (ACT, funded in January 2004), 
and two complementary experiments: Background Imaging of Cosmic Extra- 
galactic Polarization (BICEP) [200], designed for large angular scale, and 
QUest at DASI (QUaD) [201], designed for small angular scales. The second 
experiment, which is already collecting data, is composed of the Q and U Ex- 
tragalactic Sub-mm Telescope (QUEST) instrument mounted on the structure 
of the DASI experiment. A second set of experiments is scheduled in Antarc- 
tica at the French-Italian Concordia station and in the Atacama plateau in 
Chile, for unprecedented precision measurements of the i?-mode for / < 1000 
(which is particularly useful for probing the primordial gravitational waves 
from inflation): the B-modes Radiation measurement from Antarctica with a 
bolometric INterferometer (BRAIN) [202] instrument for measuring large 
scales, and the CI ObserVER^^ (ClOVER) [203] instrument for intermediate 
scales. ClOVER was approved for funding by PPARC in late 2004 and could 
be operational by 2008. 

At that time, a CMB satellite (the third one after COBE and WMAP) should 
be already collecting data: the Planck satellite [204] has already been built 
and should be launched in late 2007 or early 2008 by the European Space 
Agency (ESA). The temperature sensitivity of Planck will be so good that 
it can be thought as the "ultimate" observation in the sector of temperature 
anisotropics. However, it will still be possible to improve the measurement of 
E-polarization on small angular scales, and B-polarization will be poorly con- 
strained by Planck. On intermediate scales, the ground-based experiments 
should be quite efficient, but progress will still be needed on very large scales 
(requiring full-sky coverage) and very small scales (requiring both high reso- 
lution and excellent sensitivity). Beyond Planck, at least two space projects 
are under investigation. The mini-satellite SAMPAN (SAtellite to Measure 
the Polarized ANisotropies) [205] is targeted for large scales, and aims at 
improving the ClOVER measurement of the i?-mode for I < 1000, thanks 
to its full-sky coverage and slightly better sensitivity. The NASA wishes to 
launch a more ambitious satellite project in order to make the "ultimate" 
measurement of E-polarization, like Planck for temperature, and a very good 
measurement of B-polarization on all scales, with for the first time an instru- 
mental noise smaller than the B-mode cosmic variance up to / ~ 1000 or even 
maybe / ~ 1500. The generic name of the NASA call for projects is Infla- 
tion Probe . Its design is still very uncertain. A very preliminary project 



^ http : //www . hep . upenn . edu/ ~angelica/act/act . html 
http : //www. astro . caltech. edu/-^lgg/bicepjEront .html 
http: //www. astro . cf . ac .uk/groups/instrumentation/proj ects/quad/l 
http : //apc-p7 . org/APC _CS/Experiences/Brain/ index . phtm^ 



http : //www-astro . physics . ox . ac . uk/-^act/clover . html 
http : //sci . esa. int/science-e/www/area/ index . cfm?f areaid=17 
http: //universe .gsfc .nasa.gov/prograiii/inf lation. html. 



was called CMBpol; this name is mentioned in some works presented there- 
after [197,206,101], with sketchy characteristics mentioned in Table 4. There 
are now more advanced projects under investigation. Nobody knows which 
project and which design will be eventually approved. We mention in the 
Table some hypothetical characteristics for Inflation Probe suggested in 
Ref. [198], assuming a satellite project with a bolometer array and a passively 
cooled telescope of 3-4 m of aperture. This is one of the various possibilities 
which are being discussed, and it should be regarded as purely indicative. 

Some predictions concerning the sensitivity of these experiments (not com- 
bined with other observations) to the neutrino mass are published in Refs. 
[196,207,101,198] for Planck and in Ref. [198] for the combinations BICEP 
+QUaD, BRAIN+CIOVER, as weU as SAMPAN alone, Planck+SAMPAN 
and the hypothetical version of Inflation probe (see also [208,209] for fore- 
casts including a large neutrino asymmetry, now disfavoured as we saw in 
Sec. 3.1). Like for current bounds, forecast errors are not straightforward to 
compare from paper to paper, because different authors use 

• different fiducial models with different fiducial values for the cosmological 
parameters. In Ref. [196] , the fiducial model which is closer to the current 
concordance ACDM model is the one in Table II of that paper; here, we 
will only quote results from this table. There are eleven free parameters, 
which are the usual six of the ACDM model, the mass of a single massive 
neutrino, the tensor-to-scale ratio T/S, the scalar tilt running a, the cur- 
vature fraction ftk and the ionized helium fraction Ine- The fiducial value 
of Ml, is very close to zero (so that the fiducial CMB spectra are indistin- 
guishable from those with massless neutrinos). In Ref. [207], the number 
of free parameters is reduced to eight, by removing {T/S, a, Ihc) which 
are fixed respectively to (0, 0, 0.24). There is a single massive neutrino with 
fiducial mass 0.07 eV. Next, in Table III of Ref. [101], there are also eight 
free parameters, but the last one is Yne instead of VL^. There are three mas- 
sive neutrinos with total mass My = 0.3 eV distributed according to the 
normal hierarchy scheme. Finally, in Table II of Ref. [198], there authors 
marginalize either over eight parameters (the usual six of ACDM model, the 
mass of a single massive neutrino and Ine) or over eleven parameters (the 
same ones, plus the scalar tilt running a, the parameter w of the equation 
of state of dark energy, and the effective number N^s of extra relativistic 
degrees of freedom). The fiducial value of the neutrino mass is M^ = 0.1 eV. 

• different experimental sensitivities. For instance, the modelling of Planck 
in Refs. [196], [101] and [198] are roughly equivalent, although in [101,198] it 

was updated to the current instrument design. Instead, in Ref. [207], the au- 
thor replaces the expression of the noise variance A^";"^"^, Nf'^ [see Eq. (151)] 
by two step functions equal to zero below some critical multipoles {l^^, 
^max) = (2500,1500) and to infinity above. This is equivalent to assuming a 
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Table 5 

Forecast error a(Mi,) in eV using Planck alone: compared results from the liter- 
ature. The third column shows the free parameters over which the final result is 
marginalized ("7" stands for six usual parameters of the minimal ACDM model, 
plus total neutrino mass M,^). The fourth column gives the assumed fiducial value 
of the mass in eV. In the last two lines, the two numbers correspond to optimistic or 
pessimistic assumptions concerning the foregrounds contamination of the primary 
C MB signal. 



Ref. 




parameters 


fiducial M^, 


Planck sensitivity 


[196] 


0.3 


7+{a, T/S, Ok, YHe} 





slightly optimistic 


[207] 


0.07 


7+{J7k} 


0.07 


very optimistic 


[101] 


0.3 


7+{rHe} 


0.3 


up-to-date 


[198] 


0.45-0.49 


7+{yHe} 


0.1 


up-to-date 


[198] 


0.51-0.56 


7+{a, w, iVeff, Yne} 


0.1 


up-to-date 



slightly over-optimistic version of Planck. 

• different parameter basis and different finite steps for the computation of 
the derivatives dC^^' /d9i. The likelihood C (x\9j can be close to a multi- 
variate Gaussian in one basis, and more complicated in another basis - in 
which case Eq. (155) may give only a poor approximation of the Fisher 
matrix. The possible non-gaussianity also implies that the numerical value 
of the derivative dC^^ /dOi depends on the step A^j used for its calculation. 
However, authors usually optimize their choice of basis and step size [see 
Ref. [196])] in order to obtain robust results. 

In order to illustrate how these different assumptions are translated into differ- 
ent forecasts from paper to paper, we summarize in Table 5 the predictions for 
Planck published by the four references [196,207,101,198]. In the case of Ref. 
[198] we quote two numbers, which represent some optimistic or pessimistic 
assumptions concerning the subtraction of astrophysical foregrounds from the 
observed CMB maps: these two values, which are not very different from each 
other, are supposed to bracket the true realistic Icr error. The difference be- 
tween a(M^) = 0.3 eV in [196] and a{M^) = 0.07 eV in [207] (taken from 
Fig. 3 of this reference) probably comes from the sketchy and very optimistic 
modelling of the noise in the second reference. The results of [196] and [101] 
seem to be consistent with each other, since in the latter reference the number 
of parameters is smaller but the experimental characteristics are more real- 
istic. The difference between a(M^) = 0.3 eV in [101] and a{M^) ^ 0.45 eV 
in [198], for identical cosmological models and experimental characteristics, 
comes from variations in the fiducial value of the total neutrino mass, which 
decreases from 0.3 to 0.1 eV. Indeed, Ref. [101] explores the dependence of 
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Table 6 

Forecast error a{My) in eV using various CMB experiments or combinations of 
them. These numbers are all taken from Ref. [198], which assumes either an eight 
parameter model (the usual six ones of ACDM plus and ine), or an eleven 
parameter model (adding a, w and iVcfr). The fiducial value of the total mass is taken 
to be Mv = 0.1 eV, and in each case the two numbers correspond to optimistic or 
pessimistic assumptions concerning the foregrounds contamination of the primary 
CMB signal. 



Experiment 8 parameters 11 parameters 



BICEP+QUaD 


1.3 ■ 


■ 1.6 


1.5 ■ 


■ 1.9 


BRAIN+CIOVER 


1.5 ■ 


- 1.8 


1.7 - 


■ 2.0 


Planck 


0.45 - 


- 0.49 


0.51 - 


■ 0.56 


SAMPAN 


0.34 - 


- 0.40 


0.37 - 


- 0.44 


Planck+SAMPAN 


0.32 - 


- 0.36 


0.34 - 


■ 0.40 


Inflation probe 


0.14 - 


■ 0.16 


0.25 ■ 


■ 0.26 



a{My) on the fiducial value My (see Fig. 22), and find that when My decreases 
from approximately 0.5 to 0.1 eV, a{My) is multiphed roughly by a factor two. 
Finally, the comparison of the two results obtained in Rcf. [198] shows that 
the determination of the neutrino mass by Planck alone does not suffer sig- 
nificantly from parameter degeneracies which appear when a, w or A^^cfr arc 
included as extra free parameters. Similarly, the comparison of [196] and [101] 
suggests that this also applies for the possible inclusion of T/S and Qk- 

For other CMB experiments than Planck, not combined with other type of 
observations, the only available reference is [198], whose results are summa- 
rized in Table 6. We see in this table that ground based experiments alone are 
not helpful for constraining the neutrino mass: this is due to that fact that 
with a small sky coverage, and in spite of their excellent sensitivity, they have 
relatively large error-bars (with respect to Planck) for small and intermedi- 
ate values of /. We will see in Sec. 6.4 how this conclusion changes when one 
assumes that the results of these experiments are analyzed with the "lensing 
extraction" technique. 

In summary, the satellites Planck, SAMPAN or the combination of both 
would reach a la error in the range a{My) = 0.3 — 0.4 eV, depending on the 
case. This is of the same order of magnitude as the current neutrino mass 
bounds presented in section 5, when all available CMB and LSS data are 
combined with each other. Finally, if the Inflation probe mission can actually 
reach the sensitivity assumed in Table 4, it will provide a sensitivity ranging 
from a{My) = 0.15 eV to 0.25 eV depending on the number of marginalized 
parameters. 
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6.3 Future galaxy redshift surveys 

As already discussed in Sec. 5.3, Galaxy Redshift Surveys (GRS) allow for 
a reconstruction of the matter power spectrum modulo a light-to-mass bias 
factor. The analytic models for non-linear structure formation show that the 
relation between the reconstructed power spectrum and the total matter lin- 
ear power spectrum should be scale-independent, Pobs(^) = b^P{k), at least 
up to a wavenumber /cmax where non-linear corrections start to induce scale- 
dependent biasing. At redshift zero and for the concordance ACDM model, 
this wavenumber is expected to be around A^ax = 0.15/iMpc~^ The usual 
approach is to discard any information for k > k^s_y^ in order to avoid many 
complications related to biasing, non-Gaussian statistics and uncertainties in 
the predictions of the non-linear matter power spectrum. 

A Fisher matrix can be derived for future surveys following the lines of Sec. 
6.1. It was shown by Tegmark [210] that the Fisher matrix takes a very simple 
form under a couple of assumptions. 

• For k < kuiax, one assumes that the power spectrum reconstructed from 
the data, Pobs{k), is really a tracer of the total matter power spectrum, 
not significantly affected by non-linear corrections or by redshift-space dis- 
tortions. Wc recall that what is measured is the redshift of each galaxy, 
not its physical distance; therefore, strictly speaking, the power spectrum is 
not reconstructed in three-dimensional real space but in redshift space; the 
redshift-space and real-space power spectra do not coincide exactly, they are 
related through a function which depends on the cosmological parameters, 
and this complication is ignored here in first approximation. 

• One assumes that the survey is volume-limited rather than brightness- 
limited, i.e. the spectrum is reconstructed from data within a volume V 
in which the selection function n(r), which represents the number den- 
sity of observed galaxies, is approximately constant (instead of being sup- 
pressed near the edges due to the fact that only the brightest galaxies would 
be seen). Moreover, one assumes that the density is large enough so that 
l/n{r) < P{k) for all relevant wavenumbers fcmin < k < /cmax, which means 
that the galaxy correlation function really probes cosmological information 
rather than shot noise due to insufficient sampling. 

Under these assumptions, the Fisher matrix takes a form which is analogous 
to its CMB counterpart of Eq. (155): 




(157) 
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where the weight function w{k) of the survey is defined as 



w{k) 



V V 



(158) 



A3 (27r/fc)3 ■ 



In Eq. (157), the factor A7[w{k)dlnk accounts for cosmic (or samphng) vari- 
ance: it gives the number of independent measurements with wavenumber 
in the range [k,k + dk] within the survey volume. This factor is neghgible 
for wavenumbers smaller than fcmin, the smallest wavenumber probed by the 
survey: so, for practical purposes, the integral can be performed starting from 
zero. The result of Ref. [210] is actually more general than above, and includes 
the case where shot noise is not negligible and where the selection function is 
not constant inside the volume V. In these cases, one must replace V in Eqs. 
(157) and (158) by an effective volume Ves{k) defined as 



However, for the future surveys discussed in this section, we can systematically 
assume that the power spectrum reconstruction will be performed inside a 
volume V chosen in such way that throughout this volume l/n{r) is smaller 
than Pohs{k) for all fcmin < k < A;max, so that Vcs{k) can be replaced by V. 

For instance, for the completed SDSS Bright Red Galaxy (BRG) survey, this 
volume will include galaxies up to a distance A ~ 1 h^^Gpc with a sky cov- 
erage /sky = 0.25. So, one can take V ~ /sky(4/3)7rA3 = 1 {Gpc/hf [196]. 
Beyond SDSS, there are various plans for larger surveys, like for instance the 
SDSS-II LEGACY project, the Advanced Large, Homogeneous Area Medium 
Band Redshift Astronomical survey (ALHAMBRA) [211], the Dark Energy 
Survey ^"^ (DES) [212,213,214], or the Kilo- Aperture Optical Spectrograph^^ 
(KAOS). Some of these surveys will go to such high redshift that it will be 
possible, first, to compare the observed power spectrum with the linear one 
up to wave-numbers significantly larger than the usual fcmax ~ O.lbh Mpc~^ 
(since the non-linear evolution affects larger and larger modes as time passes 
by); and second, to do some tomography, i.e. to reconstruct P{k) in various 
redshift bins, corresponding to different times in the evolution of the Uni- 
verse, in order to measure the variations of the linear growth factor and to 
get some good handle on the dark energy variables. Tomography is also ideal 
for neutrino mass extraction, since massive neutrinos induce a very pecu- 
liar redshift-dependence on the matter power spectrum, as we have seen in 

http : //alhambra .iaa.es: 8080/' 
http : //www . darkenergysurv ey . org/| 
http : //www . noao . edu/kaos^ 




(159) 
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Sec. 4.6.4. For instance, KAOS could build two catalogs centered around red- 
shifts z = 1 and z = 3, corresponding roughly to k^a.^ ~ 0.2 h Mpc~^ and 
^max ~ 0.48 h Mpc~^ respectively, instead of k^ax ~ 0.1 /i Mpc~^ for current 
surveys (conservative values). In both catalogs, the number density would be 
such that l/n ~ P{kraa.x)i and the effective volume of the two samples close 
to 1^ ~ 0.5 (Gpc//i)3 and V ~ 0.6 {G^c/ hf respectively . 

In order to forecast the sensitivity of a galaxy survey combined with a CMB 
experiment, one can simply add the two Fisher matrices and invert their sum. 

The sensitivity of large redshift surveys to the neutrino mass was fist studied 
in Ref . [94] ; however this pioneering paper was rather crude in the way to com- 
bine GRS with CMB experiments. This work was updated and generalized in 
[196], and later in [207,101,102]. In Table 7, we quote from these papers the 
forecast error a{My) obtained assuming that the completed SDSS survey is 
combined either with the WMAP or Planck data. These results are based on 
Table II in [196], Figs. 4 and 5 in [207], and Table III in [101], for a common 
conservative value fcmax = 0.1/iMpc~^. With respect to the results quoted in 
the previous subsection, the number of free parameters has increased by one 
unit, since there is now an unknown free bias parameter h = [Pobs(^)/-P(^)]^^^ 
which is marginalized out. We see that SDSS alone gives essentially no informa- 
tion on M,^ given the various parameter degeneracies. For the SDSS+WMAP 
combination, there is an excellent agreement between [196] and [207], and 
the forecast error a{My) = 0.3 eV is three times smaller than the current 
error based on incomplete WMAP+SDSS results a{Mu) = 0.9 eV [135]. For 
SDSS+Planck, it seems again that the Planck noise modelling in [207] is 
much too optimistic, while Refs. [196] and [101] agree on a{M^) = 0.2 eV. If 
we compare this with the result obtained in [101] under the same assumptions 
for Planck alone, a{My) = 0.3 eV, we see that SDSS allows for a significant 
improvement. 

Concerning experiments beyond SDSS, Ref. [207] and later Ref. [101] present 
forecast errors in a very general way, by simply showing a {My) as a function 
of two free parameters V and fcmax- Instead Ref. [102] assumes some detailed 
characteristics for three high-redshift surveys, each of them operating in three 
redshift bins. Here we will summarize the results of Ref. [101], assuming that 
an hypothetical future redshift survey with given (V, fcmax) is combined with 
a CMB experiment that could be SPTpol, Planck, CMBpol (see Table 4 
for the value of assumed sensitivity parameters), or an "ideal CMB experi- 
ment" that would be limited only by cosmic variance up to / = 2500 (both 
for temperature and polarization). The main difficulty for reaching this goal 
would be to subtract accurately small-scale foregrounds (point-like sources. 



^ The characteristics of KAOS are taken from the "Purple Book" available on-line 
at http : //www . noao . edu/kaos 
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Table 7 



Forecast error a{Mi,) in eV for the completed SDSS survey, with volume V = 
l(Gpc//i)^, eventually combined with WMAP or Planck data: compared results 
from the literature, assuming the conservative value femax = 0-1^ Mpc"^. 



Ref. 


SDSS 


+WMAP 


+Planck 


assumptions 


[196] 


9 


0.3 


0.2 


12 free parameters, 
optimistic Planck characteristics 


[207] 




0.3 


0.06 


9 free parameters, 
very optimistic WMAP/Planck charact. 


[101] 


7 




0.2 


9 (other) free parameters, 
up-to-date Planck characteristics 



dusty galaxies, etc.) but even with current technology such an ideal experi- 
ment is not unconceivable. For these four cases. Fig. 21 shows the predicted 
2a error on the total neutrino mass M^. In these plots the fiducial value of 
the total mass was fixed to = 0.11 eV, and distributed according to the 
NH scheme. For SDSS (or for any survey with 2; < 1) we expect the relevant 
value of /cmax to be around 0.15 /i Mpc~^. However, depending on the over- 
all amplitude of the matter power spectrum (often parametrized by cTg, and 
still poorly constrained) and on future improvements in our understanding of 
non-linear corrections, this value might appear to be cither too optimistic or 
too pessimistic: this is the reason why it is interesting to leave it as a free pa- 
rameter. The figure shows the importance of employing high-volume surveys, 
which have the potential to improved the forecast errors even for the best 
CMB experiments. Keeping fcmax = 0.15 /i Mpc~^, we see that a two-sigma de- 
tection of the minimal mass in the IH scenario would require V = 40 (Gpc//i)^ 
after Planck, V = 12 (Gpc//i)2 after CMBpol, while the "ideal CMB exper- 
iment" alone would suffice. A two-sigma detection of the minimal mass in the 
NH scenario seems to be unreachable for any plausible survey volume. 

As in the previous subsection, we should mention that the results for a{Mi,) 
depend significantly on the assumed fiducial value M^,. This is illustrated in 
Fig. 22, taken from Ref. [101]. For the combination of either Planck with 
SDSS, or of the "ideal CMB experiment" with a larger hypothetical survey, 
we can see that the 2a error is typically twice larger for M^, ~ 0.1 eV than for 
- 0.5 eV. 

It is legitimate to wonder whether the most precise combination of experiments 
would be sensitive to different neutrino mass splittings between the three 
families, for a fixed value of the total mass. Indeed, we have seen in Sec. 4.6 
that unlike CMB observables, the power spectrum of large scale structure 
is in principle sensitive both to the total and the individual values of the 
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Fig. 21. Results from Ref. [101] on the neutrino mass sensitivity of future redshift 
surveys, combined with future CMB experiment. The level of gray gives the 2a" ex- 
pected errors on the total mass M^, (eV) for a fiducial value of 0.11 eV, as a function 
of the largest wavenumber fcmax used in the analysis [in units of hMpc~^] and of 
the survey volume V [in units of Vsdss = 1 (Gpc/h)^]. The vertical lines indicate 
the cut-off wavenumber fcmax for the linear matter power spectrum at the conser- 
vative (optimistic) value 0.15(0.2) h Mpc^^. The thin contours are (from bottom to 
top) for 0.3 and 0.2 eV, while the thick contours correspond to the minimum values 
of Ml, in the IH (lower lines) and NH (upper lines) schemes, assuming a squared 
mass difference Amatm = 2-6 x 10~^ eV^ (thick solid lines). The thick dashed lines 
correspond to the same limits for a value Am^^jj^ = 3.7 x 10~^ eV^ which is now 
essentially ruled out by current data on neutrino oscillations. 

three masses. The analysis of Ref. [101] demonstrates however that the tiny 
difference between the power spectrum obtained in the NH and IH cases for 
a fixed value of the total mass (close to 0.1 eV) is far too small for being 
detected even with the combination of a huge redshift survey of volume V — 
40 {Gpc/hf with the "ideal CMB experiment" . 

The issue of degeneracies between the neutrino mass and other parameters, 
such as those describing the dark energy evolution, probably needs further 
investigation. Assuming that dark energy has a constant equation-of-state pa- 
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fiducial M (eV) fiduciai IVI (eV) 

Fig. 22. Predicted 2a error on the total neutrino mass as a function of the as- 
sumed fiducial value of M^,, using: (left) Planck+SDSS, (right) an ideal CMB 
experiment (limited only by cosmic variance up to / = 2500, both for temperature 
and polarization) and a redshift survey covering 75% of the sky up to 2 ~ 0.8 
[VcS = 40 (Gpc/h)^]. Both plots assume a limit fcmax = 0.15 /i Mpc~^ and a squared 
mass difference Am^^^ = 2.6 x 10"'^ eV^. In each case, we show the results assuming 
either NH or IH. Figure from Ref. [101]. 

rameter w, Ref. [170] showed that future analyses of combined CMB and LSS 
data could be affected by a significant degeneracy between Mi, and w, a result 
which was not found in Ref. [101]. The authors of Ref. [86] focused on this 
degeneracy, including a possible dependence of w on redshift. They point out 
that in order to help resolving the degeneracy, it is useful to take into ac- 
count the fact that the CMB and LSS data are not statistically independent, 
as assumed for simplicity in all the works cited in this section. Instead, it is 
well-known that due to the late Integrated Sachs- Wolfe (ISW) effect (see Sec. 
4.4.7), the power spectrum of temperature anisotropics encodes some infor- 
mation on neighboring structures (like galaxy clusters). So, there is a non-zero 
cross-correlation between temperature and galaxy maps, which has actually 
been already measured [215,216], still with a rather low level of significance. In 
the future, precise measurements of this cross-correlation could be very useful 
for neutrino mass extraction. 



6.4 CMB weak lensing 
6.4-1 Principle 

It was recently realized that the CMB anisotropies encode even more cos- 
mological information than expected, because it should be possible in a near 
future to measure the deflection field caused by the weak lensing of CMB pho- 
tons by the large scale structure of the neighboring universe up to a redshift 
2; ~ 3 [217,218,219] (see [220] for a recent review on CMB lensing). The power 
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spectrum of the deflection field encodes some information concerning struc- 
ture formation mainly in the linear or quasi-linear regime, and is therefore 
extremely useful for measuring parameters like the total neutrino mass (or 
the dark energy equation-of-state) which mildly affect the primary anisotropy 
[197]. So, the next generation of CMB experiments could output for free a LSS 
power spectrum, without suffering from the systematics induced by mass-to- 
light bias and by strong non-linear corrections on small scales at ^ < 0.2, as 
is the case for galaxy redshift surveys. 

More precisely, weak lensing produces a deflection of light rays parametrized 
by a two-dimensional deflection field d(n) such that an object located in a 
direction n will appear in a direction n + d(n). At leading order [221,222], the 
deflection field can be written as the gradient of a lensing potential, d = V0. 
The lensing potential contains information on the gravitational potential along 
the photon trajectory, convolved by a window function which depends on the 
distance to the source 

(f){n;rs) = J dr W{r;rs) ^{rh,r] — —r) , (160) 





where r is the comoving coordinate along the line of sight, Vg is the source 
coordinate, $ the gravitational potential, and rj the conformal time defined 
by dt — —a dr] and 77 = today (with this definition, the photons emitted by 
the source and reaching us today travel according to r{rj) = —rj). Using the 
geodesic equation, it is possible to show that the window function W{r; r^) is 
equal to 

W(r;rs) = 2^^^^^ . (161) 



In the present case, the source is the last scattering surface and stands for 
its comoving radius. For a given cosmological model, the primordial spectrum 
(|$|^) and its time evolution can be calculated; the same numerical code which 
computes the CMB anisotropy spectrum C{^^ and the matter power spectrum 
P{k) can also easily compute the deflection fleld power spectrum Cf'^, which 
encodes some information on LSS at redshifts of order ^ ~ 3 or so. Hence, if 
we could map the CMB deflection flcld d, we could measure Cf'^ and compare 
it with theoretical predictions in order to constrain cosmological parameters. 

There are several methods on the market for extracting the deflection map 
[221,223,224,225], all based on the non-Caussianity induced by lensing [217]. 
These methods start from the assumption that both the primary anisotropics 
and the deflection fleld are Gaussian; they also assume that the noise present 
in the temperature and polarization maps is Gaussian and uncorrelated with 
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the signal. None of these assTimptions is exactly true. Ref. [226] estimated to 
which extent the lensing extraction will be biased, first, by the non-Gaussianity 
of the lensing potential caused by the non-linear growth of matter perturba- 
tions on small scales, and second, by the imperfect cleaning of the CMB maps 
from the kinetic Sunyaev-Zel'dovich effect, which also has a blackbody spec- 
trum, induces non-Gaussianity, and features spatial correlations with many 
of the structures responsible for the lensing. Both effects were found to be 
relevant (i.e. to induce a significant bias in the estimators). However, they are 
small enough to preserve the vahdity of the method. Ref. [227] estimated the 
impact on lensing extraction of a small level of non-Gaussianity produced in 
the Early Universe by inflation, reheating or some other mechanisms (for a 
review, see [228]), and concluded that given the current bounds on primor- 
dial non-gaussianity inferred from CMB maps, this non-gaussianity cannot be 
large enough to threaten lensing extraction. 

6.4-2 Quadratic estimator method 

This field is rather technical, but we can stress the main ideas without entering 
into many details. For simplicity, we will write the main formulas in the flat- 
sky approximation, which amounts in replacing the expansion of CMB maps 
in harmonic space by an expansion in two-dimensional Fourier space, and to 
use two-dimensional vectors 1 instead of pairs of integers (/, m). In the flat-sky 
approximation the Fourier modes are fully described by the power spectra (7"^ 
where a and b belong to the basis {T, E, B} of temperature, T-mode polariza- 
tion and £^-mode polarization. Weak lensing correlates the lensed multipoles 
[229,217] according to 

(a(l)6(l'))cMB = {2nyS{l + l')Cf + /"^(l, l')0(l + 1') (162) 

where the average holds over different realizations (or different Hubble patches) 
of a given cosmological model with fixed primordial spectrum and background 
evolution (i.e. fixed cosmological parameters). In this average, the lensing po- 
tential (/) is also kept fixed by convention, which makes sense because the CMB 
anisotropics and the LSS that we observe in our past light-cone are statisti- 
cally independent, at least as long as we neglect the integrated Sachs- Wolfe 
effect. The above functions are defined in [224]. For instance, for ab — TT 
one has 

/^^(l, V)C = C^i^ + 1') ■ 1 + Cj^i^ + 1') • 1' • (163) 

Let us introduce the basic principle of the quadratic estimator method of Hu 
& Okamoto [223,224,225] (which is equivalent in terms of precision to the 
alternative iterative estimator method of Hirata & Seljak [221] as long as 
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CMB experiments will make noise-dominated measurements of the B-mode, 
i.e. at least for the next decade). By inverting Eq. (162), one builds a quadratic 
combination of the temperature and polarization observed Fourier modes 

d^'iL) = Jj^AWWih, h) (164) 

where I2 = L — li, and in which the normalization condition 

-1 

(165) 



ensures that d"*" is an unbiased estimator of the lensing potential 

(d»''(L))cMB = ^L0(L) = d(L) . (166) 

Note that, so far, the coefficients 5f"^(li,l2) are still arbitrary. From the ob- 
served temperature and polarization maps, one could compute each mode of 
d"^ and obtain various estimates of the deflection modes, precise up to cosmic 
variance and experimental errors. In order to quantify the total error, it is 
necessary to compute the power spectra of the quadratic estimators 

(d"'*(L)d''''(L)) = (27r)25(L - L')Cf ^"''^ (167) 



/Tr(li,l2)y"'(li,l2) 



where the average is now taken over both CMB and LSS realizations, since 
0(L) is also a stochastic quantity. In this definition, the power spectra are writ- 
ten with a superscript dd{ab) in order to be distinguished from the actual power 
spectrum of the true deflection fleld. These spectra feature the four-point cor- 
relation function of the observed (lensed) Fourier modes (0(11)6(12)0(13)6(14)), 
which should be expanded at order two in 0(L) in order to catch the leading 
non-Gaussian contribution. 

The four-point correlation functions are composed as usual of a connected and 

an unconnected piece. The connected piece is by definition a function of the 
power spectra Cf** in which we now include all sources of variance: cosmic 
variance, lensing contribution and experimental noise. The unconnected piece 
is a function of the same spectra plus the deflection spectrum Cf^, and as usual 
it can be decomposed in three terms corresponding to the different pairings 
of the four indexes [232]: (/i,/2), (^3,^4) or (/i,/3), {l2,U) or (/i,//i), {h,h)- 
The first term leads to considerable simplifications when it is plugged into 
the expression of the quadratic estimator power spectrum, and the result is 
simply Cf'^, as one would expect naively from squaring Eq. (166). The other 
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terms lead to more complicated expressions that we will write as a noise term 

which represent respectively the contribution from the connected piece and 
from the two non-trivial terms of the unconnected piece [230,232]. In order 
to get an efficient estimator, one should adopt the set of coefficients ^"^(li, I2) 
which minimize the noise terms. It is actually much easier to minimize the 
connected term only, which leads to the simple results 



for a = b (for a ^ b see [224]). With such a choice, the unconnected piece 
contribution iV£^ can be shown to be smaller than Af^, but not completely 
neghgible [230]. 

The various estimators d"'' can be constructed for each pair of modes, except 
for the pair BB, because the spectrum C^^ is dominated by lensing at least on 
small scales, which invalidates the present method. Therefore, the quadratic 
estimator technique would not be optimal for long-term CMB experiments 
with cosmic- variance-dominated measurement of the B mode [221,231]. For an 
experiment of given sensitivity, the five other estimators can be combined into 
a final minimum variance estimator, which gives the best possible estimate of 
the deflection fleld by weighing each estimator accordingly to its noise level. 
The sensitivity of the Planck satellite is slightly above the threshold for 
successful lensing extraction, but only at intermediate angular scales, and 
with essentially all the signal coming from the d^^ estimator. The following 
generation of experiments - such as the CMBpol or Inflation probe project - 
should obtain the lowest noise level from the d^^ estimator [224]. 



6.4-3 Neutrino mass from CMB lensing extraction. 

The authors of Refs. [197] and [198] studied the limits that could be obtained 
on neutrino mass by using future CMB experiments only, assuming that the 
lensing power spectrum could be extracted using the quadratic estimator tech- 
nique. In this way, it is possible to combine information on CMB acoustic 
oscillations and on the surrounding LSS of the Universe at redshifts of order 
2; < 3, avoiding many of the complicated features of galaxy redshift surveys 
related to the mass-to-light bias and to the strongly non-linear evolution of 
small-scale perturbations at small redshift. 

The principle of Refs. [197,198] is to treat the power spectrum of primordial 
anisotropics (for temperature and polarization) and of the deflection field as 
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two independent information on the cosmological parameters. The expected 
error-bars on the temperature and polarization multipoles are estimated from 
the characteristics of a given CMB experiment, while those on on the lens- 
ing multipoles are taken to be the dominant term [N^f' in Eq. (168), which 
can be computed from the same experimental characteristics using essentially 
Eqs. (169) and (165). Then, the authors assume a fiducial model and per- 
form a Fisher matrix analysis in the same fashion as described in the previous 
subsection. 

The fiducial model of Ref. [197] includes ten free parameters: the usual six 
parameters of the minimal ACDM model, the neutrino mass M^, in the case 
of three degenerate neutrinos, the equation of state parameter w, the running 
of the tilt a and the primordial Helium abundance Ine- The authors conclude 
that the Planck satellite will output a la error-bar a{My) ~ 0.15 eV, which 
shows that the use of the lensing extraction technique should improve the 
global sensitivity in a rather spectacular way at least for this parameter. In 
the case of the CMBpol project an impressive error of a{Mi^ — 0.044 eV is 
predicted. 

This study was repeated in Ref. [206], with one extra parameter representing 
the variation of the dark energy equation-of-state with respect to the scale 
factor. Wc saw before that this parameter could potentially introduce extra 
degeneracies involving the neutrino mass, but the forecast errors are found 
to be approximately the same: oihly) = 0.16 eV and a{M^) = 0.046 eV for 
Planck and CMBpol, respectively. 

Finally, the authors of [198] assumed two possible fiducial models, one with 
eighth free parameters (the same as in Ref. [197] with the exception of a 
and w) and one with eleven free parameters (same as in Ref. [197] plus the 
effective number of extra relativistic degrees of freedom A^eff)- They consider six 
experiments or combinations of experiments as listed in Table 8, and for each 
case they quote two numbers, corresponding cither to optimistic or pessimistic 
assumptions concerning the subtraction of astrophysical foregrounds from the 
observed CMB maps. 

It can be seen in Table 8 that the combination QUaD-l-BICEP benefits a lot 
from lensing extraction, since the error a{M^) decreases from approximately 
1.5 eV to at least 0.4 cV. Thus, with QUaD+BICEP it should be possible 
to reach in a near future -using CMB only- the same precision that we have 
today combining many observations of different types (galaxy-galaxy correla- 
tion function, Lyman-a forest) which are affected by various systematics. The 
situation is almost the same for BRAIN+CIOVER, which should also achieve 
a{M^) ~ 0.4 eV using lensing extraction. Planck should make a decisive 
improvement, lowering the error to a{Mi,) ~ 0.15 eV, in excellent agreement 
with the results of Ref. [197]. Note that without lensing extraction the error 
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Table 8 

Forecast error a(Mi,) in eV using various CMB experiments or combinations of 
them, and assuming that the lensing map is extracted following the quadratic esti- 
mator method. These numbers are all taken from Ref. [198], which assumes either 
an eight parameter model (the usual six of ACDM plus Mi, and Ihc); or an eleven 
parameter model (the same ones plus a, w, Nes). The fiducial value of the total 
mass is taken to be Mi, = 0.1 eV, and in each case the two numbers correspond to 
optimistic or pessimistic assumptions concerning the foregrounds contamination of 
the primary CMB signal. 



Experiment 8 parameters 11 parameters 



BICEP+QUaD 


0.31 - 


0.36 


0.36 - 


0.40 


BRAIN+CIOVER 


0.34 - 


0.43 


0.42 - 


0.51 


Planck 


0.13 - 


0.14 


0.15 - 


0.15 


SAMPAN 


0.10 - 


0.17 


0.12 - 


0.18 


Planck+SAMPAN 


0.08 - 


0.10 


0.10 - 


0.12 


Inflation Probe 


0.032 - 


0.036 


0.035 - 


0.039 



would be multiplied by three (by four in the case with extra free parameters). 
No significant difference between the forecast errors in the eight and eleven 
parameter models was found in Ref. [198]. SAMPAN alone is slightly more 
efficient than Planck, and the combination Planck+SAMPAN is the first 
one to reach a{Mj,) ~ 0.1 eV, even in the pessimistic case of large foreground 
residuals and extra free parameters. Finally, the version of the Inflation 
Probe satellite considered in [198], which is slightly more ambitious than 
the CMBpol assumption of Ref. [197], is able to reach a{My) ~ 0.035 eV 
both in the eight and eleven parameter cases. In any case, the results seem 
to be relatively robust against pessimistic assumptions concerning foreground 
contamination. 



6.5 Galaxy weak lensing (cosmic shear) 
6. 5. 1 Principle 

Weak lensing changes the apparent shape of galaxies in various ways. Among 
others, it changes slightly their apparent ellipticity: i.e. galaxies which would 
be apparently spherical without weak lensing can look in fact elliptical, stret- 
ched in one direction and squeezed in the orthogonal direction. This effect, 
which is coherent over the angular size of the lensing gravitational field, i.e. 
potentially over many sources, is called cosmic shear. It can be detected pro- 
vided one has a very dense sample of galaxies with enough resolution in order 
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to measure each individual shape. In order to be able to reconstruct the lens- 
ing gravitational field, it is also necessary to know cither the redshift of each 
source galaxy, or at least their coarse-grained number density g{z,n) in red- 
shift space and in direction n. Let us assume for simplicity that we have an 
isotropic sample with a spherical number density g{z) (or g{r), where r is the 
comoving radius). We can divide the sample in small solid angles or pixels. 
Since intrinsic ellipticities are randomly distributed, the average ellipticity in 
each pixel will give a measurement of the average cosmic shear, up to a shot 
noise term proportional to N^^J , where Npi^ is the number of galaxies per 
pixel. Therefore, the size of the pixels must be chosen optimally in a trade-off 
between high angular resolution and high signal-to-noise. 

The cosmic shear field reconstructed by this method derives from the same 
lensing potential (p as the deflection field discussed in the previous subsection. 
The relation between the lensing potential and the gravitational potential 
$ along the line of sight is the same as in Eqs. (160) and (161), except that 
one must now average over all individual sources of comoving distance G 
j-g^^maxj ^^iQ coutinuous limit and for an ensemble of sources described by 
the density function g{r), this gives 

„max 
' s 

j drs g{rs) <j)in;r,) 

(t>{n;g) = — — 

j drs g{rs) 





j,max 
' s 

j dr W{r;g) ^{rh,7] = -r) , (170) 





with 



■ s 

j drs g{rs) W{r;rs) 
W{r- g) = ' . (171) 

drs g{rs) 







Like for CMB anisotropics, once we have the lensing potential map 0(n; g) 
we can compute the power spectrum Cf'^'^ , which gives us an information 
on the surrounding large scale structure projected along the line of sight. 
Therefore, cosmic shear surveys offer another opportunity to measure the two- 
point correlation function of total matter fluctuations, complementary to the 
three-dimensional Fourier power spectrum P{k) obtained from galaxy redshift 
surveys (probing lower redshifts) and to the CMB lensing power spectrum 
(probing higher redshifts). 

If a shear survey contains a large enough number of galaxies with known 



104 



redshifts, it is even possible to refine this technique by sphtting the sources 
into redshift (or distance) bins, from < < ri to r„_i < < r„. It is 
then possible to measure n lensing potential maps (l){n^gi), such that the bin 
distribution gi{rs) vanishes outside of the range rj_i < < rj. From these 
maps one can obtain n power spectra cf'^'^" and n{n — l)/2 cross-correlation 
spectra cf'^'^'^^ (cross-correlations do not vanish because the same structure 
can contribute to the lensing of sources in different bins). Since the highest 
redshift bin probes structures at larger redshift, this method allows to some 
extent for a three-dimensional reconstruction of the gravitational fluctuations. 
This technique is called tomography, and brings some hope of tracking the 
evolution of P{k, z) as a function of redshift caused by dark energy [233] and 
by neutrino masses (see Sec. 4.6.4). 



6.5.2 Neutrino mass from cosmic shear surveys 

There are various ongoing and planned cosmic shear surveys, like for instance 
the Canada- France-Hawaii Telescope Legacy Survey (CFHTLS) [162]; the 
Dark Energy Survey (DES) [212,213,214]; the SuperNova Acceleration Probe 
(SNAP); the Panoramic Survey Telescope And Rapid Response System 
(Pan-STARRS); or the Large Synoptic Survey Telescope (LSST) [234,235]. 
The sensitivity of this type of observation to neutrino masses has been ad- 
dressed first in [236], and later in [237,206,238]. Here we will review the re- 
sults of [206], where four hypothetical experiments are considered. The first 
two, called S300 and SIOOO, stand for future satellite surveys with an average 
source density of 100 galaxies per squared arc- minute and a sky coverage of 
300 square degrees (/sky = 0.0073) and 1000 squared degrees (/sky = 0.024), 
respectively. The SNAP lensing survey could be close to one of these two as- 
sumptions. The last two, called G27r and G47r, stand for future ground-based 
surveys with an average source density of 65 galaxies per square arc-minute 
and a sky coverage of respectively /sky = 0.5 (as expected for LSST) and 
/sky = 1 (as expected for Pan-STARRS). The authors make some realistic 
assumption concerning the redshift distribution of the sources (given experi- 
mental limitations) and perform a Fisher matrix analysis for the combination 
of the cosmic shear surveys with either Planck or CMBpol (see Table 4), 
including in each case information on the CMB lensing, extracted with the 
quadratic estimator method described in Sec. 6.4. They assume a tomographic 
analysis of the galaxy-ellipticity data with eight redshifts bins. Thanks to the 
quadratic estimator method, the CMB experiments also provide a measure- 
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Table 9 

Forecast error a{Mj^) in eV from reference [206], using some hypothetical cosmic 
shear surveys (S300, SIOOO, G27r and G47r) described in the text, combined with 
either Planck or CMBpol (including, in each case, information on the CMB lensing, 
extracted with the quadratic estimator method). 





S300 


SIOOO 


G27r 


G47r 


Planck 


0.12 


0.090 


0.052 


0.045 


CMBpol 


0.047 


0.045 


0.031 


0.027 



ment of the cosmic shear power spectrum, sensitive to higher redshifts: here 
the information from CMB lensing plays the role of a ninth redshift bin. 

There arc eleven free cosmological parameters in the analysis of Ref. [206]: 
the usual six ACDM parameters, plus the total neutrino mass (attributed 
to a single family), the primordial helium fraction, the running of the scalar 
tilt, and two parameters describing the dark energy equation of state: w{a) = 
Wq + (ao — a)wa. The results arc summarized in Tabic 9, from which one 
concludes that cosmic shear surveys constitute an extremely powerful tool 
for constraining the total neutrino mass. The sensitivities obtained with this 
method are the absolute best forecasts derived so far. This is not surprising, 
given the ideal properties of weak lensing observations: precise, high redshift 
measurement of the matter power spectrum over a wide range of scales, and 
tomography. 

Following Ref. [206] and Table 9, the most ambitious hypothetical experiment 
GAti combined with the Planck data would reach (7{M^) ~ 0.045 eV, i.e. 
the same sensitivity as CMBpol alone (see Sec. 6.4.3); while C47r combined 
with CMBpol would reach a{M^) ~ 0.027 eV, which means that a 2a detec- 
tion of the neutrino mass would occur even for the smallest M^, in the NH 
scheme, of order (Am^tm)^^^ — 0-05 eV. The results of Table 9 are relatively 
conservative in the sense that in the Fisher matrix computation, the cosmic 
shear power spectrum was limited to the angular scales 6 > 0.18 deg (in a 
multipole expansion, / < 1000) below which non-linearities render the power 
spectrum reconstruction more difficult. However, the authors of Ref. [206] re- 
fer to a method for extracting information on smaller scales proposed by Refs. 
[239,240]. Ref. [206] speculates that using this method, one could gain a fac- 
tor three (two) with respect to the results of Table 9 for the combination 
of Planck (CMBpol) with G27r or G47r. This would imply a real precision 
measurement of the total neutrino mass with cosmological observables. 

Finally, the recent analysis in Ref. [238] considers the combination of Planck 
(no lensing extraction) with two lensing surveys inspired by the SNAP and 
LSST projects, using up to five tomography bins. The results are consistent 
with those of Ref. [206]. 
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6. 6 Galaxy cluster surveys 



The last method that we will discuss is expected to bring some extra inde- 
pendent information, and could further improve the bounds discussed in the 
previous subsection. 

In the past years, the study of the cluster abundance evolution (the function 
A'"(^), where N is the number density of galaxy clusters and z the redshift) has 
lead to interesting constraints on the normalization of the matter power spec- 
trum at small scales (defined in Sec. 5.5), and on the matter density fraction 
fim- For a given cosmological model, the function N{z) can be predicted from 
non-linear simulations, and compared with experimental data from galaxy 
cluster surveys. The non-linear simulations generally introduce systematic er- 
rors, but since cluster scales are not too deep inside the non-linear regime even 
today, the method can be kept under control. 

The function N{z) -or better, its derivative dN/dz- is affected by the evolution 
of the linear growth factor on cluster scales. Therefore, as emphasized recently 
in [241], future precise measurements of cluster abundances will be useful 
probes of the neutrino free-streaming effect. 

On the observational side, the prospects for galaxy cluster surveys are very 
promising, since catalogs will be derived from various techniques, by mapping 
the X-ray emission of the hot cluster gas, or the scattering of CMB photons in 
this gas (Sunyaev-Zel'dovich effect, SZE), or the cosmic shear measured from 
galaxy weak lensing (as discussed in the previous subsection). We have already 
seen that ambitious projects arc being planned for galaxy weak lensing; this is 
also the case for X-ray and SZE observations. The full data will then provide: 
first, a measurement of the matter power spectrum P{k, z) (modulo some 
bias and redshift-space distortion effects), which can be measured in various 
redshift bins, like for weak lensing tomography; and second, an estimate of 
the abundance evolution dN/dz, also at various redshifts. 

An expression of the Fisher matrix describing future P{k,z) and dN/dz{z) 
observations has been proposed in [242]; this forecast method was applied to 
neutrino mass extraction in [241], using the expected instrumental sensitivity 
of the following projects; for CMB temperature and polarization anisotropics, 
Planck (no lensing extraction); for X-ray observations, the Dark Universe 
Observatory (DUO) survey (see [243] for the description of a similar, previ- 
ously proposed survey), which should see ~ 11,500 clusters in 6,150 deg^; 
for SZE observations, the South Pole Telescope (SPT) survey [199], yielding 
~ 20, 000 clusters in 4,000 deg^; and for galaxy weak lensing, the LSST project 
(already introduced in Sec. 6.5.2), assuming ~ 200, 000 clusters in 18,000 deg^. 
Among the three data sets, Ref. [241] finds that the most powerful should be 
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Fig. 23. Forecast la sensitivities to the total neutrino mass from future cosmological 
experiments, as described in Sec. 6, compared to the values in agreement with 
present neutrino oscillation data in Eq. (2) (assuming a future determination at 
the 5% level). Left: sensitivities expected for future CMB experiments (without 
lensing extraction), alone and combined with the completed SDSS galaxy redshift 
survey. Right: sensitivities expected for future CMB experiments including lensing 
information, alone and combined with future cosmic shear surveys. Here CMBpol 
refers to a hypothetical CMB experiment roughly corresponding to the Inflation 
Probe mission. 

the LSST lensing survey: the forecast error a(M^ is typically 50% larger for 
SPT, and 60 to 100% larger for DUO. For the combination Planck+LSST, 
it is interesting to note that the inclusion of dN/dz improves the neutrino 
mass determination by a factor two: using only information from the CMB 
and from P(/c, z), the authors find a{My) ~ 0.09 eV, while adding dN/dz they 
get a{My) ^ 0.04 eV. Including the DUO and SPT data marginally improves 
this bound. 



7 Conclusions 



Cosmology has played an important role in constraining neutrino properties, 
providing information on these elusive particles that complements the efforts of 
laboratory experiments, and neutrino cosmology will remain an active research 
field in the next years. 

In this review, we were interested in the connection between neutrino masses 
and cosmology. We have described in detail how the evolution of cosmological 
perturbations is modified by massive neutrinos, and how their effects leave an 
imprint in the different cosmological observables. We focused on the standard 
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case of three active neutrino species, known to be massive thanks to the ex- 
perimental results on flavour neutrino oscillations which, however, do not fix 
the absolute neutrino mass scale. 

We saw how the analysis of cosmological observables can provide either a pos- 
itive indication or an upper bound on a total neutrino mass in the range from 
0.05 to a few eV indicated by oscillation data and other laboratory results. 
We discussed why a unique cosmological bound does not exist: the values 
depend on the assumed cosmological model and the considered cosmological 
data. Within a minimal AMDM model, the present cosmological bounds at 
95% CL are given in Tables 1-3, ranging from 2-3 cV with CMB data only to 
more stringent values when LSS data are added: 0.9-1.7 cV considering only 
the shape of the matter power spectrum from 2dF or SDSS data, or 0.3-0.9 
eV fixing the light-to- mass bias and/or adding Lyman-a forest data and/or 
the SDSS measurement of the baryon oscillation peak. These values hmit the 
neutrino masses in the so-called degenerate region, as presented in Fig. 19. We 
also discussed how these cosmological bounds could be relaxed or even dis- 
appear when adding some extra cosmological parameters or in scenarios with 
non-standard relic neutrinos. It is interesting to note that these mass bounds 
also apply for other thermal relic particles, such as axions with masses in the 
eV range as shown in [244] . 

Finally, we have discussed how the sensitivity to neutrino masses will improve 
with new cosmological data, such as new CMB experiments or larger galaxy 
redshift surveys. In particular we have reviewed the excellent prospects of 
future weak lensing measurements, either from CMB data or with cosmic shear 
surveys. We saw how these future data will be sensitive to a total neutrino 
mass well below 1 eV, even when considering extra cosmological parameters. 

We give a graphical summary of the forecast sensitivities to neutrino masses 
of different cosmological data in Fig. 23, compared to the allowed values of 
neutrino masses in the two possible 3-neutrino schemes. One can see from this 
figure that there arc very good prospects for testing neutrino masses in the 
degenerate and quasi-degenerate mass regions above 0.2 eV or so. A detection 
at a significant level of the minimal value of the total neutrino mass in the 
inverted hierarchy scheme will demand the combination of future data from 
CMB lensing and cosmic shear surveys, whose more ambitious projects will 
provide a 2a sensitivity to the minimal value in the case of normal hierarchy 
(of order 0.05 eV). The combination of CMB observations with future galaxy 
cluster surveys (derived from the same weak lensing observations, as well as 
X-ray and Sunyaev-Zel'dovich surveys) should yield a similar sensitivity. 

The information on neutrino masses from analyses of cosmological data is 
complementary and can not replace the efforts in terrestrial projects such as 
tritium beta decay and neutrinoless double beta decay experiments (and vice 
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versa). In particular, a positive result from one of these three possibilities, 
which as reviewed in Sec. 2 do not measure the same combination of neutrino 
masses and mixing parameters, should be cross-checked by the others. Of 
course any information on the absolute neutrino mass scale will be a very 
important input for theoretical models of particle physics beyond the Standard 
Model. 

Note that light massive neutrinos could also play a role in the generation of 
the baryon asymmetry of the Universe from a previously created lepton asym- 
metry. In these leptogenesis scenarios, one can also obtain quite restrictive 
bounds on light neutrino masses, which are however model-dependent (see 
[245] for a recent review) . 

This review was entirely devoted to the case of light neutrinos with masses 
of order eV or smaller. However, we would like to mention the case of a ster- 
ile neutrino with a mass of the order of a few keV's and a small mixing 
with the flavour neutrinos. It has been shown that these keV neutrinos could 
also have interesting cosmological consequences: they would play the role of 
dark matter and replace the usual CDM component (see [246,247,248,249] 
for an incomplete list of references) . Because of their large thermal velocity 
(slightly smaller than that of active neutrinos), they would behave as Warm 
Dark Matter and erase small-scale cosmological structures. Their mass can 
be bounded from below using Lyman-a forest data from quasar spectra, and 
from above using X-ray observations [250,251,252,253,254]. This scenario, al- 
ready disfavoured according to the results in [254], is currently under careful 
examination. 



Note added 



A few new analyses of cosmological data including non-zero neutrino masses 
appeared after this review was accepted for publication. In Rcf. [255] the 
authors used data on the Lyman-a forest from SDSS (improved thanks to 
the inclusion of the mean flux constraints derived from a principal component 
analysis of quasar spectra), SDSS and 2dF galaxy clustering, SNIa and CMB 
power spectra (WMAP3 and others) to get an upper bound of M^, < 0.17 eV 
(at 95% CL) on the total mass of the three active neutrinos {Mi^ < 0.26 eV 
for four thermalized species). Instead, in Rcf. [256] the limit M^y < 2.0 eV for 
three neutrinos was found using WMAP3 data only, in agreement with the 
result of the WMAP team [65]. 
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